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Global  Tomography  with  Irregular  Parameterisatioii 
(Towards  regiorialised  travel-time  corrections 
for  body  waves) 
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and 

Malcolm  Sambridge 


Introduction 

Global  models  of  velocity  heterogeneity  have  emerged  in  the  past  decade  which  show 
significant  deviatioits  from  spherical  symmetry  in  Earth’s  structure.  As  the  quality  of 
data  and  their  coverage  increases  more  detail  is  included  in  the  models  Global 
models  currently  claim  to  resolve  structures  as  small  as  1000  km  in  lateral  extent,  but 
regional  models  resolve  structures  to  considerably  smaller  scales,  i.e.  a  few  hundred 
kilometres,  'fhe  strength  of  heterogeneity  in  the  models,  i.e.  the  peak-to-peak  variation 
of  velocity,  increases  as  more  detail  is  included.  The  smength  of  heterogeneity  in  the 
P-wave  velocity  field  reaches  ±5%  at  the  top  of  the  upper  mande  at  the  scales  which 
are  resolved  regionally.  It  is  clear  that  such  strong  heterogeneity  will  significantly 
affect  the  travel  times  of  the  body-wave  phases  most  commonly  used  in  earthquake 
location.  Accounting  for  this  effect  can  therefore  significantly  improve  the  quality  of 
event  location  which  is  of  obvious  importance  for  nuclear  monitoring. 

It  is  clear  from  models  derived  from  body  waves  and  surface  waves  that  the  level 
of  velocity  heterogeneity  in  the  Eai'th  vaires  dramatically  with  depth  and  is  sti'ongly 
concentrated  in  the  top  300-400  km  (e.g.  Gudmundsson  et  al.  [1990]).  Most  of  the 
effect  of  structure  on  mavel  times  is  therefore  due  to  structures  at  relatively  shallow 
depth  through  which  teleseismic  body  waves  have  a  relatively  simple  ray  geometry 
(steep).  The  image  of  upper  mantle  heterogeneity  that  has  emerged  is  one  that  clearly 
correlates  with  the  overall  tectonic  processes  on  the  surface.  Young  oceanic  regions 
have  a  low  velocity  as  do  regions  beneath  tectonically  active  continents,  while  old  oce¬ 
anic  regions  and  old  stable  continents  are  underlain  by  a  mantle  of  high  velocity. 

Accounting  for  the  effects  of  velocity  heterogeneity  on  travel  times  and  thus 
event  location  must  include  the  main  features  of  the  velocity  structure  and  remain  a 
simple  procedure  in  order  to  be  implementable.  We  attempt  to  meet  these  prere¬ 
quisites  with  the  following  strategy: 

We  use  the  high  level  of  correlation  of  global  velocity  models  with  surface  tec¬ 
tonics  to  define  the  main  features  of  Earth  structure.  A  detailed  regionalisation 
defines  irregularly  shaped  bodies  on  the  Earth’s  surface  which  in  turn  define 
discrete  model  parameters.  With  irregular  sampling  of  the  regionalisation  we  can 
represent  irregular  bodies  and  variable  scales.  The  sampling  is  dense  where 
needed  (continental  margms,  island  arcs  and  tectonic  continent),  and  sparse  where 
that  suffices  (intraregional,  oceanic). 
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We  use  Delaunay  tesselation  and  Voronoi  cells  around  each  sample  of  the  region¬ 
alization  to  define  discrete  regions  belonging  to  the  various  tectonic  types.  This 
allows  us  to  collaps  the  number  of  parametres  needed  tp  represent  the  complex 
regions  significantly.  Efficient  algorithms  are  available  for  navigation  and  book 
keeping  in  an  iiiegular  grid  (Delaunay  tetrahedra  or  Voronoi  cells)  (Sambridge  et 
al.,  1995). 

We  use  a  select  data  set  of  travel  times  from  well  located  events  and  nuclear 
explosions  in  order  to  alleviate  the  mapping  of  mislocation  into  structure.  We 
invert  these  data  for  an  irregularly  parameterised.  regionalised  upper-mantle 
model  (RUM).  Our  objective  is  not  to  find  the  best  possible  model,  but  to  find 
the  simplest  model  that  still  contains  enough  information  about  structure  to  be 
useful  for  generating  corrections  to  travel  times  used  in  event  location, 

Teleseismic  travel  times  can  then  be  corrected  for  this  model  by  a  relatively  sim¬ 
ple  look-up  table  based  on  a  tabulation  of  the  tectonic  regions  in  which  stations 
and  events  fall. 
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Method 

Figure  1  shows  our  regionalisation.  It  is  primarily  based  on  maps  published  by 
Sclater  et  al.  (1980)  and  Jordan's  (1981)  GRTi  regionalization,  but  also  incorporates 
information  from  e.g.  Mtyashixo  et  al,  (1982).  Eight  tectonic  regions  are  defined.  The 

oce.an  is  divided  into  three  age  provinces,  the  same  as  the  GRTI  i  gionalization  of  Jor-  • 

dan  [1981],  This  regionalization  has  a  resolution  of  two  degrees  compai'ed  to  five 

degrees  in  GRTI.  The  continents  ai'e  divided  into  5  regions  according  to  age,  four 

regions  of  low  volcanic  and  earthquake  activity,  in  addition  to  a  region  termed  tectonic. 

continent  where  activity  is  high.  We  have  termed  the  regions  Young  Continent, 

Intermediate-age  Continent,  Old  Continent  and  Ancient  Continent  The  chosen  age 
boundaries  are  listed  in  Figure  1  and  were  selected  by  comparing  the  age  zonation  of 
the  continents  with  tomographic  images  of  the  upper  mantle  (c.g,  Woodhuuse  and 
Dziewonski,  1984).  We  have  perhaps  unconventionally  termed  the  Iceland-F'aeroe- 
Greenland  ridge  continental  because  of  its  highly  anomalous  crustal  thickness  for  oce¬ 
anic  crust  (25-35  kin  according  to  Bott  and  Gunnarson,  1980). 

The  regionalization  is  reprerented  by  16200  paiameters.  Roughly  a  quarter  of 
those  aie  selected  to  pararneterise  the  boundaries  between  llie  regions  (4100  parame¬ 
ters).  The  boundaries  ai'e  all  parameterised  with  a  tv/o-degree  resolution  except  the 
boundaries  between  oceanic  regions.  Thus,  the  density  ot  parameteiisation  is  high 
where  we  expect  potential  sharp  contrasts,  i.e.  continental  margins  and  boundaries 
between  tectonic  and  atectonic  regions,  and  low  where  we  don’t,  i.e,  within  the  oceans. 

Figure  2a  shows  the  Voronoi  cells  around  the  selected  samples  of  the  regionalization. 

The  figure  contains  about  4100  Voronoi  cells.  We  reproduce  the  full  complexity  of 
the  regionalization  in  Figtue  1  with  only  about  a  quarter  of  the  number  of  parameters 
in  the  original  (see  Figure  2b).  Each  Voronoi  cell  is  taken  to  encompas,s  a  region  of 
uniform  velocity  at  each  depth.  The  regionalization  is  extended  from  the  surface  to  a 
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depth  of  660  km  to  span  the  whole  of  the  upper  mantle. 

Subducting  slabs  are  signihcant  structures  in  the  upper  mantle  which  are  clearly 
oblique  and  cannot  be  included  in  a  two-dimensional  regionalization  as  above.  Most 
of  the  earths  seismicity  occurs  in  or  around  subducting  slabs.  Slabs  are  thus  likely  to 
su  ongly  affect  travel  times  in  general.  Wc  have  included  oblique  slab  structures  in  our 
parameterization  of  the  upper  mantle  as  the  ninth  region.  First  we  contoured  slab- 
related  seimicity  world  wide.  Examples  of  slab  contoures  are  presented  in  Figure  3a 
for  the  slabs  of  the  northwest  Pacific  and  the  Tonga/Kermadec  slab.  The  contours  are 
drawn  at  fifty  km  depth  intervals  and  such  that  they  lie  near'  the  top  of  the  seismogenic 
region.  We  used  the  relocated  NEIC  catalog  of  Engdahl,  van  der  Hilst  and  Buland  to 
define  the  seismicity.  The  catalog  is  not  complete  and  contains  significant  gaps  which 
are  difficult  to  Interpret.  Such  gaps  are  not  interpolated  if  they  exceed  a  few  hundred 
km  in  size.  The  slab  contours  obviously  do  not  define  a  slab  where  there  is  no  seismi¬ 
city.  It  is  difficult  to  invisage  the  thermal  anomaly  of  a  subducting  slab  to  terminate 
as  abruptly  as  the  seismicity  Li  many  places.  Thus  oui'  model  of  slabs  is  likely  to  be 
incomplete.  However,  slabs  are  defined  where  we  most  need  them,  i.e.  where  the 
eai'thquakes  occur  for  which  travel-time  corrections  are  needed. 

The  slab  contours  define  a  surface  when  inteipolated  which  we  take  to  define  the 
top  of  the  slab  (or  what  used  to  be  the  top  of  the  oceanic  lithosphere  before  subduc- 
tion).  We  then  define  a  complimentai'y  surface  assuming  the  slab  is  200  km  thick. 
This  is  a  simplification  of  slab  geometry,  particularly  at  depth.  Projecting  the  original 
surface  200  km  along  its  local  noiinal  results  in  the  folding  of  the  original  surface. 
Such  folds  were  simply  edited  out  interactively  with  a  graphics  tool.  We  sample  the 
contours  discretely  at  a  spacing  of  about  10  km.  Our  samples  are  not  internal  to  the 
slabs  as  we  define  them,  but  on  their  surface.  Thus  constant  velocity  Voronoi  cells  arc 
not  suitable  to  describe  their  complex  geometry.  Instead  we  place  the  slabs  in  a  global 
coordinate  system  and  consti'uct  the  three-dimensional  Delaunay  tesselation  for  the 
23,000  slab  samples  for  all  24  individual  slab  stnictures.  Those  Delaunay  tretrahedra 
which  Join  the  phases  of  the  same  slab  sti'ucture  then  define  the  volume  of  slabs  in  the 
upper  mantle.  Two  examples  are  shown  in  Figures  3b  and  3c.  Those  are  two  of  the 
larger  and  deeper  slabs  which  display  significant  complexity  in  their  shape. 

In  summai-y  the  parameterization  of  the  upper  mantle  consists  of  two  distinct  ele¬ 
ments;  Delaunay  tetrahedra  to  di  .e  slabs  and  a  two-dimensional  grid  of  Voronoi  cells 
to  define  tectonic  regions.  The  slabs  are  overlaid  on  the  regionalization.  To  determine 
which  type  of  environment  a  given  point  (along  a  ray)  is  in  we  first  evaluate  if  it  is 
within  a  slab.  Lf  it  is  not,  its  position  in  the  regionalization  is  determined.  These  tools 
allow  us  to  represent  complex  .structures  efficiently  and  economic  algorithms  are  avail¬ 
able  to  locate  any  arbltrai'y  point  in  the  iiregular  grid. 

We  use  a  select  data  set  of  travel  times  from  well  located  events  and  nucleai' 
explosions  in  order  to  alleviate  the  mapping  of  mislocation  into  structure.  The  events 
are  the  105  even'.s  of  the  test-data  set  used  by  Kennett  and  Engdahl  (1991)  supple¬ 
mented  with  92  large  earthquakes  to  get  a  more  even  sampling  of  the  globe  and  the  9 
tectonic  regions.  The  criterion  for  the  selection  of  these  events  was  that  they  be 
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recorded  by  a  large  number  of  globally  distributed  stations.  Their  locations  are  taken 
from  the  relocated  lTn;iC  catalog  of  Engdahl,  van  der  Hilst  and  Buland  (1996)  where 
they  have  reassociated  ISC  data,  included  pP  and  pwP  phases  in  the  location,  used 
IASP91  (Kennett  and  Engdahv  1991)  and  subjected  locations  to  stringent  fitness  cri¬ 
teria.  We  use  only  teleseismic  V  and  PKP  anivals  (first  arrivals)  and  picks  reported  to 
the  ISC  as  impulsive.  We  invert  these  data  for  an  irregularly  parameterised,  regional- 
ised  Earth  model. 

Figure  4  shows  a  map  of  the  197  events  used  together  with  the  regionalization. 
The  events  are  grouped  into  three  categories.  First  there  are  21  nuclear  and  chemical 
explosions  displayed  as  white  stars  in  the  figure.  Then  there  are  39  earthquakes  which 
aie  well  located  with  local  observations  in  a  local  velocity  model.  They  are  displayed 
as  white  dots.  The  remaining  137  events  are  located  using  teleseismic  observations. 
The  locations  for  the  21  explosions  and  39  well  located  events  are  taken  fixed.  For 
the  remaining  137  events  we  use  locations  from  Engdahl  et  al.  (1996). 

We  refererence  tiavel-time  residuals  to  AX  135  (Kennett  et  al.,  1995)  and  reasso¬ 
ciate  all  impulsive  picks  reported  to  the  ISC  for  the  above  events.  We  find  a  suipris- 
ing  number  of  events  with  a  significant  mean  residual.  This  is  surprising  because 
AK135  is  almost  identical  in  the  rriantle  to  1ASP91  which  Engdahl  et  al.  (1996)  used 
for  theii'  locations.  The  explanation  must  lie  in  a  difference  in  data  selection.  We  for 
example  include  PKIKP  and  not  pP.  We  do  not  associate  pwP.  The  large  number  of 
events  with  a  significant  apparent  origin-time  shift  led  us  to  allow  origin  time  to  be  a 
free  parameter  in  our  inversion  for  regionalized  structure. 


Results  and  conclusions 

We  solve  a  linearized  inverse  problem.  A  higher  order  approach  is  not  warranted 
because  ray  bending  critically  depends  on  velocity  gradients  whicli  are  artificially 
represented  by  our  parameterization.  We  sample  depth  ten  times  at  66  km  spacing  and 
seek  9  independent  velocity  profiles,  one  for  each  of  tlic  eight  tectonic  regions  as  well 
as  slabs.  We  use  a  standai'd  dainped-least-squares  inversion  procedure.  In  regions 
with  limited  depth  distribution  of  sources  teleseismic  rays  give  limited  depth  resolution 
in  the  upper  mantle.  We  cannot  well  resolve  ten  independent  parameters  except  in 
subducting  slabs,  where  there  is  a  distribution  of  events  in  depth.  vSince  the  level  of 
heterogeneity  in  the  upper  mantle  is  strongly  concentrated  near  the  surface  (e.g.  Gud- 
mundsson  et  al.  1990),  in  the  top  300-400  km.  we  have  projected  the  depth  profiles  on 
to  a  set  of  4  Gaussian  basis  functions  which  are  not  alloyed  to  put  structure  in  the 
I  )wer  half  of  the  upper  ma  tie.  This  is  with  the  exception  of  slabs  which  are  allowed 
to  have  structure  to  the  base  of  the  upper  mantle.  The  result  we  present  are 
moderately  damped.  An  undamped  inversion  allowing  for  10  idependent  depth 
parameters  in  each  region  achieves  a  19%  variance  reduction.  Limiting  the  degrees  of 
freedom  to  four  per  region  via  the  Gaussian  basis  reduces  the  variance  reduction  to 
16%.  We  have  damped  the  solution  to  yield  a  variance  reduction  of  14%. 
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The  results  of  inversion  are  shown  in  Figure  5.  Each  profile  represents  a  given 
tectonic  region.  The  curves  have  error  bars  equivalent  to  one  standard  deviation 
assuming  that  the  eiror  variance  of  the  data  is  about  1  second  (Gudmundsson  et  al. 
1990),  The  results  are  encouraging.  Slabs  appear  as  a  continuous  fast  anomaly 
throughout  the  upper  mantle  except  at  the  top  where  we  have  tectonic  cmst  rather  than 
a  true  slab  and  where  the  result  is  probably  affected  by  the  slow  mantle  wedge.  The 
amplitude  of  the  slab  structure  is  2-3%  velocity  perturbation.  This  is  assumed  to  be 
uniform  over  a  thickness  of  200  km  and  would  thus  be  consistent  with  a  peak  pertur¬ 
bation  in  the  core  of  the  slab  of  4-5%.  Other  features  of  the  results  are  also  encourag¬ 
ing,  The  oceans  are  slow  with  old  oceans  developing  a  fast  litliospheie.  Tectonic  con¬ 
tinent  and  young  continent  are  also  slow  and  essentially  indestinguishable  while 
intermediate-age  continents  have  developed  a  fast  lithosphere.  There  is  no  resolvable 
structure  underneath  old  continents  suggesting  that  AK135  (and  1ASP91)  ai'e  represen¬ 
tative  of  continent  in  that  age  range.  Ancient  continents  aic  fast,  but  not  to  a  great 
depth,  only  about  200  km.  It  should  be  noted  that  we  have  not  corrected  for  crustal 
structure  so  the  model  should  be  inteipreted  to  include  crustal  signature.  We  refer  to 
this  model  as  the  RUM  model  (Regional  Upper  Mantle). 

Our  procedure  of  data  reduction  concludes  with  the  calculation  of  station  correc¬ 
tions  (where  sufficient  numbeis  of  picks  are  reported)  computed  after  correction  for  the 
above  model.  We  can  then  divide  the  reduction  process  into  three  steps:  relocation 
(ongin  time),  regionalized  tomography,  and  subsequent  station  corrections.  Table  1 
summarizes  the  vaiiance  and  mean  reduction  achieved  by  these  three  steps  in  the 
analysis,  Relocation  accounts  for  much  of  the  variance  (more  than  50%).  The  shuc- 
tural  model  achieves  14%  variance  reduction.  It  is  worth  noting  that  we  work  with 
individual  observations  and  have  not  subjected  the  data  to  any  averaging  procedures 
prior  to  inversion  (such  as  summary  rays).  Finally,  station  corrections  accomplish  a 
16%  variance  reduction.  This  indicates  that  our  model  leaves  significant  structural  sig¬ 
nal  in  the  data.  Histograms  of  residuals  before  and  after  analysis  arc  presented  in 
figure  6. 

Figure  7  presents  the  results  for  the  station  corrections.  Clearly  a  significant  sig¬ 
nal  is  left  in  the  data  after  inversion  for  regionalized  structure.  Note  for  example  that 
a  fast  anomaly  remains  within  the  Laurentian  craton.  On  the  other  hand  the  African 
cratons  are  left  with  slow  anomalies.  Some  tectonic  regions  ai'e  left  with  systemati¬ 
cally  fast  residuals  while  others  are  slow.  Clearly  all  regions  of  the  same  type  are  not 
the  same.  Our  parameterization  in  effect  impose.s  complex  spatial  constraints  on  this 
inversion  for  up;  er-mantle  structure  which  we  can  think  of  as  a  priori  infon,.ation,  i.e. 
the  presumption  that  the  upper  ma.ntle  Catr  be  regionalized  as  wc  have  done.  The  suc¬ 
cess  of  the  inversion  Justifies  those  constraints.  The  remaining  station  residuals  point 
to  the  fact  that  this  regionalized  description  is  not  complete,  however,  and  that  on 
order  of  half  the  structureal  signal  in  the  data  is  explained  by  it, 

Wc  have  described  upper-mantle  stiucture  with  23  degrees  of  freedom  and 
manage  to  explain  about  half  of  the  apparent  structural  variance  in  travel  time  data. 
The  data  are  individual,  not  summary  ray  data.  Preliminary  re,sults  from  S-wave 


residuals  give  results  similar  to  tire  ones  presented  here  for  P-wave  velocity  variations. 
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tectonic  regions  shown  above. 
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Black  dots  are  earthquakes  with  no  local  constraints,  but  which  have  been 
located  with  a  large  number  of  glob:".!  observations  subject  to  fitness  constraints. 


travel-time  data. 


Variance  and  mean  reduction 


mean 

[s] 

variance 

[s^l 

variance 

reduction 

Initial 

1.18 

4.30 

After  relocation 

0.69 

1.99 

54% 

After  velocity  perturbation 

0.19 

1.71 

14% 

After  station  correction 

0.04 

1.44 

16% 

Statistics  of  travel-time  residuals  at  progressive  stages  in  the 
analysis.  Much  of  the  variance  of  the  data  lies  in  earthquake 
mislocation.  A  significant  level  of  variance  reduction  is 
achieved  by  the  modeled  velocity  perturbation  although  no 
averaging  [e.g.  summary  raysl  has  been  applied  [we  use 
individual  residuals].  The  variance  reduction  achieved  by 
station  corrections  is  probably  due  to  structure  which  is 
either  on  scales  smaller  than  allowed  for,  or  which  reflects 
inconsistencies  between  regions  categorized  as  of  the  same 
tectonic  type  [see  map  of  station  corrections]. 
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SUMMARY 

The  advent  of  broad-band  seismology  has  meant  that  use  is  being  made  of  a 
vs'ide  range  of  seismic  phases  for  many  of  whicli  ellipticity  corrections  have 
not  been  readily  available.  In  particular,  when  many  seisnric  phases  are  used 
in  location  schemes,  it  is  important  that  the  systematic  effects  of  ellipticity  arc 
included  for  each  phase. 

An  efhdent  and  elieclive  procedure  for  constructing  ellipticity  corrections  is 
to  make  use  of  the  ray-based  approach  of  Dziewonksi  &  Gilbert  (1075),  as 
reformulated  by  Doornbos  (1088),  in  conjuction  with  the  rapid  evaluation  of 
tra\'el  times  and  slowness  for  a  given  range  using  the  tau-spline  procedure  of 
buland  ik  Liiaprnan  tio8,u. 

Ellipticity  coefficients  have  been  tabulated  for  a  wide  range  ol'  seismic  phases 
and  are  ai’ailable  in  electronic  form.  The  ellipticity  correction  iirocedures  ha\’e 
been  extended  to  include  an  allow.mce  for  diffraction  phenomena  e.g.  l^iin, 
.S'aifi  diffracted  along  the  core-mantle  boundary.  Corrections  for  additional 
phases  cun  be  generated  by  building  the  ellipticity  coefficients  Ironi  suitable 
coniLiinalions  of  the  coefheients  for  different  phase  segments. 

Key  words:  Ellipticity  corrections,  travel  times 
Running  Title:  Elliplicit>'  corrections  for  scisniir  phases 
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1  INTRODUCTION 

As  methods  for  the  delineation  of  the  three-dimensional  structure  from  seismic  arrivals 
have  been  developed,  a  wider  range  of  seismic  phases  have  begun  to  be  used.  It  is 
important  therefore  that  the  influence  ot  tlie  ellipticity  of  the  figure  of  the  Earth  on 
seismic  travel  times  be  taken  into  accouht  in  a  comparable  way  for  each  seismic  phase. 
The  same  need  arises  in  the  exploitation  of  later  seismic  phases  in  rapid  earthquake 
location  using  a  sparse  global  array;  the  systematic  effects  of  ellipticity  must  be  inc'uded 
for  all  phases  used  in  the  location  procedure. 

The  normal  model  for  the  calculations  of  travel  limes  for  seismic  phases  is  to  use  a 
spheri  al  earth  with  a  radial  distribution  of  velocity;  aspherical  effects  are  then  included 
via  a  perturbation  treatment.  yVn  important  component  of  such  aspbericity  arises  from  the 
ellipticity  of  the  figure  of  the  Earth.  The  resulting  intluenee  on  the  travel  time  of  a  seismic 
phase  will  then  depend  on  the  location  of  the  source  and  receh'er  points,  rather  ilian  jus; 
the  cpireniral  distance  A  between  them,  for  a  source  at  co-latitude  ^0,  and  a  path  with 
a/miuth  C  to  the  receiver  (measured  clockwise  from  North),  Bullen  (ItW7a)  introduced 
triple  entry  correction  tables  <it  ^  b((i>(,,  A,C)  to  e\aluate  the  ellipticity  correchon  lor  /' 
waN'e.s;  these  were  later  extended  to  polar  laiitudes  (cf  Engdahl  &  Gunst,  lOtiti). 

Dziew'onski  &  Gilbert  ( 1  i)7())  ha\  e  dtwcloped  an  elegant  representation  ol  the  inlliience 
of  ellipticity  on  trac'el  times  and  have  demonstrated  the  significance  ol  source  depth,  z.s, 
in  deiermininu  the  appropriate  eorreciions.  l-'or  a  seismic  phase  the  correciioii  can  be 
built  u])  I  roin  tlirei'  spherical  harmonic  components  of  order  2, 

i)t  ^  (r»i  l-tt.', ,  A)/'j.„i  (cos  f>(i)  cos  mC.  ll.U 

in  0 

file  coefficients  A)  represent  ray  path  integrals  in  terms  of  the  radial  \  eloeil\ 

disthlmiion  v(r)  as  an  iiiiegral  over  epicentral  distance.  Iiilroducliii',  as  a  scaled 

associated  Legendre  function, 

;iTj,,„(cosh),  (1.2) 

the  coeflicienls  Ai  take  tlie  form 

(r,„(c,,A)  -  [  dtb]'(/')-,--c(r),\,„(b) 

r  Jo  dr 

-v]cA,„(d)  I  1- jtA,„(())  (tj-:  -p’)''  [(].;t) 

where 

']('■)  iMii/r,  (1.4) 

p  is  the  ray  slowness,  c[r)  is  Earth's  ellipticity  of  figure,  and  the  effect  of  source  depth 
is  implicit  in  the  dependence  of  r  on  0  along  the  ra>'.  The  sutnmations  represent  the 
inlluence  of  discontinuities  in  velocity  (v)  in  transmission  (./)  and  reflection  ik).  Uoornbos 
(1988)  has  shown  liow  the  integrals  for  (;,„(Z(,A)  can  allernativi'ly  be  I'orrnulaied  in 
terms  of  integrals  over  radius  without  sacrificing  the  convenience  of  the  representation 
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(1.1).  Note  that  the  elliptictty  coefficients  were  designated  by  by  Dziewon.ski  &  Gilbert 
(197S),  we  have  changed  to  the  notation  oy,,  to  avoid  conflict  with  the  well  established 
usage  T(p)  for  the  delay  time  wliich  is  also  used  in  section  2.1. 

Dziewonski  &  Gilbert  (1975)  presented  a  tabulation  of  ellipticity  corrections  using  the 
(/„, (z^,A)  formulation  for  the  phases  P,  PcP,  l'’A.T(ab,bc.dn.  SKS  using  the  pern 

model  of  Dziewonski,  Hales  &  Lapwood  (197,5).  Corrections  were  given  for  depths  of  0, 
300  and  650  km.  llie  perturbations  in  travel  time  are  fortunately  relatively  insensili\  e  to 
the  particular  velocity  model  which  is  employed.  Kennett  (1991)  extended  the  ct,„[Z;,A) 
tabulation  using  the  iasp91  model  of  Kennett  &  Kngdahl  (1991)  to  include  a  liner 
tabuhition  in  depth  (z.,-  0,  35,  50,  10,  200,  300,  500  and  700  km)  and  additional  pihases 

,,  SeP  and  SKP. 

In  this  paper  we  demonstrate  the  convenience  of  combining  the  Biilnnd  &  Chapman 
(1983)  tau-spline  procedure  for  calculating  travel  times,  as  used  by  Kennett  &  Dngdahl 
(19i)l),  with  the  Doornbos  (1988)  calculation  of  the  a,„  ellipticity  coefficients  using 
ray  integrals.  ,\n  extensive  range  of  ellipticity  coefficients  have  been  tabulated  and  are 
available  in  electronic  form.  The  correction  procedures  have  been  extended  to  include  an 
allowance  for  diffraction  phenomena,  as  e.g.  at  the  core-mantle  boundary  in  i'aifi, 

1  urtlicr,  ue  diunonstraie  how  to  combine  the  a,,,  coefficients  for  different  phases  to 
produce  corrections  h.ir  cases  where  the  coefficients  have  not  been  tr  ilaied. 


2  TLLIl’  l  IC'TY  COKRLCT'IONS 

In  lerins  ol  the  fn,|(z,,  Al  cuelucients  the  o.xplicit  representation  of  the  time  eorrection 
to  he  added  to  the  calue  calculated  for  a  sjiherically  symmetric  harlh  modi'l  is 

hft.W-.dn.A.iC)  :|(1  I  3cos2d„)  milZv.A) 

-  -V  sill  2  9()  cos  'C  cri  ( z, ,  A ) 

.  siir’  9,1  cos 2  C(r;.(Zs,  A),  (2,1) 

where,  as  before,  f+o  is  the  epicentral  co  latitude,  A  is  the  cpicentral  disianc'',  C  is  the 
azimuth  Irom  the  epicentre  to  the  receiver,  and  Zt  is  source  depth. 

2,1  Kay  based  corrections 

Doornbos  (1988)  has  presented  an  effective  algorilhm  lor  the  eoin]niiaiion  of  the  tr,,, 
l  ilipiuiiy  coefficients  which  depends  on  a  specification  of  the  ray  paiainelei  for  any 
particular  phase.  This  approach  is  ba.sed  on  compulation  of  ray  integrals  for  fLxed 
slowness  and  in  eonsequence  needs  to  be  coupled  to  an  auxiliary  routine  if  correcdons 
are  to  be  produced  for  fixed  epicentral  distance,  A. 

liuland  &  Chapman  (1983)  have  introduced  the  tau-spline  method  for  rapid 
construction  of  travel  times  as  explicit  functions  of  range.  This  method  is  based  on  the 
pro]3erlies  of  the  delay  time  r 


r(p)  T(p}  -  pAtp). 


(2.2) 
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The  stationar>'  values  of 

0{p,x)  Tip]  +  px  =  T(p)  +  p[x  ~  A{p]]  i2.2) 

correspond  to  geometrical  arrivals  for  the  range,  x.  With  a  tau-splinc  representation  of 
Tip)  for  each  travel-time  branch  the  condition  c9[p,x)/dp  =  0  leads  to  the  solution  of 
a  quadratic  to  find  the  slowness,  Px,  and  then  a  rapid  evaluation  of  the  trawl  time  for 
range  .v  from 

Tix]  =  Tipx)  +  PxX.  i2A) 

I'he  tau-splines  can  be  established  for  many  phases  at  a  time  by  using  a  common  slowness 
discretisation  and  the  superposition  of  t  segments  for  different  ray  segments.  The 
appropriate  software  has  been  widely  distributed  with  the  release  of  the  iaspDl  model 
(Kennett  &  hngdahl,  1991)  and  can  be  readily  adapted  to  different  radial  velocity  models. 

I'he  cnitpai  from  the  I?uland&  Chapman (I98:i)  algoritlim  is  ideally  suited  for  use  with 
Ooornbos'  (1988)  ellipticity  corrections  since  the  slowTiess  is  known  e\plicitly  for  a  given 
range  and  phase  specification.  1  hus,  priwided  a  ray  description  is  available,  the  ellipticity 
coefficients,  tr,,,,  can  be  calculated  rapidly  and  a  cross  check  made  of  the  accuracy  of  the 
calculation  by  compariscm  of  the  estimated  and  desired  range.  For  the  model  akI35  of 
Kennett,  l.ngdah!  &  Huland  (199.1),  which  is  specified  by  linear  gradients  in  radius,  llu* 
convergoiiee  is  better  than  0.021  lor  all  the  major  phases. 

This  coinb'.naiion  (T  ray  algorithms  has  been  used  to  develop  an  exleusi\  e  set  of 
ellipticity  corrections  in  the  form  of  a,,,  coefficients  as  a  function  of  depth,  c,-,  and  range, 

,  fsi'i'  Tiibl”  1 ''  The'-'''  cnrt't't'i  inn'^  a>'(>  a'’ai|a'i>n  jt>  elnrirfi'-'ic  form  n't  r'  'I.  i",  .'.pyudix 

B. 

In  tile  application  of  such  t‘lh[iti(  ity  corrections  care  must  be  taken  wlu're  inierpolalioii 
is  applied  to  ensure  that  the  corrections  coresspond  to  the  same  clasi,  of  proiiagaiion 
path.  It  is  imiionant  to  separate  the  upgoiiig  /'  and  .V  waves  from  deep  evenis  il’,i,u  .S',,,, 
III  'fable  1)  from  the  main  P  and  1'  wl-iicb  have  a  separate  ray  siHicU'icalion. 

As  noted  by  l)/iewonski  &  tiilbert  ( 197())  the  elliptit  il>'  corrections  lor  /’  ami  ,S'  jihases 
having  comparahle  paths  scale  leiy  closely  with  a  constant  of  proponionaliiy  winch 
varies  Iroru  1.801  to  1.884,  i.e,  1.82  ±  0.011,  despite  the  considerable  variation  m  the 
ratios  oi  /'  to  .S'  wave  speed  witli  depth.  Tlii'  constant  of  ijroiiortionaliiy  is  e  ery  close  in 
the  value  1.8  proposed  by  Bullen  (l')27b).  The  simple  rescaling  works  well  for  iiianth' 
phases,  such  as  Td' and  i'tvS',  but  cannot  be  applied  where  the  propagaiion  paths  ha\'e 
(l.fferent  cliaiarier,  as  e.g.  core  phases  such  as  ]’KP  aiul  SKS. 

2.2  Diffracted  phases 

The  ellipticity  coeffici''nt.s  for  a  diffracted  phase  can  be  determined  from  the  values 
at  the  grazing  point  on  the  appropriate  interface.  For  a  given  source  depth,  z.^,  and 
cpiccntral  distance  for  the  grazing  ray,  A_,^,  we  cnn.struci  the  coefficients  for  the  grazing 
ra)’:  ( z, ,  ) ,  trf  ( Zy ,  Ary ) ,  a f  ( A„ ) . 

As  detailed  in  Appendix  A,  the  elliptic  ly  coefficients  for  diffracted  waves  are  specified 
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b>-  two  integrals  along  the  upward  leg  of  the  grazing  ray: 


1  -1  dv 

—  dOcos20(r/7-  — 
P  Ja^  dt- 
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(2.61 


r  ^  -I  £  sin20  I  (f]i  - 

'.I  “k  ( 

!he  change  in  arc  length  due  to  ellipicity  is  second  order  in  c  and  may  ihercl'ore  be 
Ignored. 

The  (Til,  coefficients  for  the  diffracted  phase  may  then  be  constructed  from  those  at 
grazing,  f/;,’,,  by 


ij Q  ( z  V ,  A,/  e  b ) 
(r|'(Zv,A_„  -1  b) 

(rytrs-.A.,  +  b) 


ail  ( s'f'  <5  (A'  sin  b  -r  y  cos  b ) , 

•> 

(T'l  ( Zj ,  A (, )  +  -“=  sin  b  ( y  sin  b  -  A'  cos  b ) , 
\'d 

(■ri'iZv,  Ai;)  — i=  sinb(A  sinb  +  I  cosb). 


(2.K) 

(2.91 


where  b  i.s  the  arc  distance  of  diffraction.  Only  two  numerical  integrals  need  to  be 
e\  alualed  to  generate  the  new  coefficients.  The  integrals  X  and  y  are  weakly  dependent 
on  source  depth  throttgh  the  variation  in  epieentral  distance  of  tlie  grazing  point, 

'Ihe  set  of  a  ctteificieiits  for  hjiff  and  .Sjiif  diffracted  at  the  core-mantle  boundary  are 
illustrated  in  figure  1  as  fitnclions  of  epieentral  distance,  together  with  those  for  I’  and 
.s'.  .As  we  e.vpeet  there  is  a  smooth  transition  from  the  phase  defined  by  a  geometric  ray 
to  the  diffraeted  phase,  for  short  diffracted  legs  a  simple  extrapolation  of  the  geometric 
results  pro\'ides  an  adequate  approximation. 

The  tr  coefficients  for  Pgiff  and  .S'a,tT  are  tabulated  in  Tables  2,2  in  compaiahli'  form  to 
that  employed  in  Keiuiett  (l!)i)l)  for  P  and  .V. 


2  COMllINAl  ION  Of  fLLIPl  ICHT  COlflffClIONS  fOR  COMl’OSI  1  L  I’lIA.SES 

Hie  form  ol  the  ellipticity  correction  (1.2)  is  the  same  for  all  phases  with  the  coefficients 
.t;„i  being  the  only  phase-dependent  quantity.  Ifius,  if  we  can  construct  the  w;we  path  of 
a  particular  phase  as  a  combination  of  two  or  more  elementary  phases  (e.g.  PKKI’  as  a 
combination  of  PKI^  and  PeP),  we  can  not  only  construct  the  correction  for  the  composite 
phase  from  the  corrections  for  the  elementary  pha.ses,  hut  also  construct  the  o„,  tables 
for  the  composite  phase  by  a  linear  combination  of  the  cr„|  tables  for  the  elementary 
pliases.  Consider,  for  example,  a  surface  rellected  phase  with  two  legs,  a  and  h,  where 
the  point  of  rohoction  lies  at  A,)  so  that  the  two  propagation  legs  are  (0,  A„)  and  (A„,  A) 
(sec  figure  2).  The  corrections  could  be  determined  by  calculating  the  coefficients  I'or 
the  A„  leg  and  then  iiitioducing  an  apparent  source  at  the  rellection  point  to  build  a 
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further  correction  for  the  arc  length  A|,  -  A  -  A,,.  Alternatively,  we  can  use  spherical 
trigonometry  to  develop  an  expression  for  the  composite  o,„  coefficients: 


cr(^'(z.v,A)  =  A„) 


cr,Mz_v,  A)  o  j'tz,-,  A„) 


(rUz,,A)  -  o-ytz^.A,,) 


+  jOcos-’ A„  -  l)cr^’(0.Ai,) 

-  ^sin2AafTf(0,Ai,) 

+  ^  sin-' AdCrj’(0,  A|,), 

4  y  sin2A,i(T//(0,  A/,) 

+  cos  2A„or|’(0,  Ab) 

-  4  sin2Ai,(j!’(0,  A(.). 

4  V*- sin-’ A|,y,)’(0,  A|,) 

4  1  sin2A„y{’((),  A/.) 

4  1  (1  4  COS"  A,()yj’(0.  Aj.), 


» 


(2.1) 


(2.2) 


where 


.A/i  -  A  —  A||. 


(2.4) 


The  deri\';ition  of  tliese  expressions  is  presented  in  Appendix  A. 

■is  ;i  quick  check  on  the  composition  rules  we  compared  a  direct  c.ilciilation  lor  the  * 

cfieflicients  for  PP  at  (iO  which  is  composed  of  two  I’  legs  of  SO"  each,  with  the  results 
Irotn  the  ap|)lit:ation  of  (2.  l)-(2.4),  The  resulting  coefficients  differ  by  no  more  tluiii  ().()()(i 
s,  which  rellects  the  truncation  errors  associated  with  tabulated  values. 

More  complex  cotubitiations  of  phases  can  also  be  mtidc  by  exploitipc  tlu'  linear  ► 

superposition  of  the  correctiotis  for  different  proptigalion  paths.  Consider  for  i  ample 
PKKP,  which  can  be  constructed  from  the  infurniation  from  PKP  anti  PcP.  for  a  full 
in'opagation  ptiih,  A,  w(‘  can  build  the  travel  titne  for  PKKP  by  ;i  coiiibinaiion  of  two  PKP 
legs  with  the  subtraclioii  of  a  PcP  segment:  ^ 

■/■'''''''''(Zv,  A)  -  4  -  ■/•'’‘■''(O.A,.),  (2.2) 

where 


A,  A,,  (A  Aj,).  (2.(d 

» 

A,)  and  Aj,  are  the  arc  lengths  of  the  two  PKI’  legs  which  have  to  be  louiul  b>-  matt  bin, g 
the  ray  parameter  (see  figure  5). 

The  ir,„  coefficients  for  PKKP  arc  then  to  be  fouiul  from  a  linear  combination  of  the 
loeflicients  lor  the  /•’AT  and  /'c/’ couiptinenls: 


^r('/^'’(z,.  A)  -  (7,';'"''(Z.s,  A.,)  -t 


PKP,, 


---  al'^'’{z,,Aa)  4 


^;(3cos"A„  -  l)lfr,'’^'’(0,A(. )  -  (0,  A,  )| 

V  sin 2A„((/, Aj.)  o  '’’ ''{(),  A,  )  I 
sin-  A„  1  <4’’^'’ ( 0.  A„ )  -  yf  '■  (0,  A, )  I , 
sin2A„y(',’((),  A{,) 

cos2A„[yf''^’(0,  A(,)  -  yf’^'’ (0,  A,  )  | 

I  sin  2A^|  y;f'^''’((),  A„)  y'’<  ''(0,  A,  )  |, 


» 


(2.7) 

» 


(3.8) 


» 


« 
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^  sin- A., (7^(0,  Aj.) 

4sin2ArtlorP*‘'''(0,  Ai,)  -  A,- )  ] 

5(1  -t- cos- Afl)[crj'^''(0,  A/,)  -  crf‘''’(0,  A,.)]. 


(3.9) 


4  Dl'^'CUSSlON 

■Although  the  significance  of  the  Earth’s  ellipticity  on  travel-times  has  long  been 
recognised,  applications  have  been  restricted  to  a  few  major  phases.  The  e.\tension  of 
the  corrections  for  diffracted  phases  and  for  combinations  of  seismic  phases  means  that 
it  is  now  feasible  to  generate  corrections  for  most  obsen'able  bodv-wa\'e  phases. 

W't  hope  that  the  provision  of  a  wide-range  of  ellipticity’  coefficients  in  readily  useable 
electronic  form  will  encourage  the  systematic  inclusion  of  the  effects  of  ellipticity  in  both 
regional  and  global  studies. 


APPENDIX  A:  COIERECTIONS  FOR  DIFFRACTED  AND  COMPOSITE  PHASES 
A1  Diffracted  phases 

The  geometry  of  a  phase  diffracted  around  a  spherical  boundary  within  the  Earth  is 
shown  and  defined  in  figure  .M.  Energy'  is  radiated  downward  from  the  source  at  S  at 
a  particular  angle  defining  a  ray  parameter,  p  -  pg,  such  that  when  the  wave  refracts 
toward  the  boundary  it  grazes  it  at  point  G,  an  arc  length  A„  away  from  the  source.  The 
ray  path  then  relracts  upwards  toward  point  P  on  the  surface  a  further  arc  length  Aj,  from 
G.  The  continuous  solid  curve  in  the  figiare  (SCP)  defines  the  grazi.ig  ray  I0  the  boundary. 
Some  energy  may  diffract  along  the  boundary'  frorn  G  to  11,  vvhere  it  is  refracted  upwards 
again  to  reach  the  surface  at  R,  an  arc  length  A  -=  A(,  -r  A;,  -1-  6  away  from  the  source.  At 
high  frequencies  the  diffracted  energy  is  confined  to  the  upper  edge  of  the  boundary  and 
has  a  predictable  travel  time,  which  is  constructed  by  a  path  integral  of  the  slowness  of 
the  Earth  along  the  path  (SGHR)  shown  in  figure  Af,  albeit  not  a  "lay"  path.  The  geometry 
of  this  path  and  the  slowness  encountered  along  it,  and  consequently  the  predicted  travel 
time,  will  be  perturbed  by  the  Earth’s  ellipticity  m  a  manner  much  like  a  true  "ray”  path 
as  formulated  by  Dzieiconski  and  Gilbert  (1976). 

W'e  can  divide  the  path  SGHR  (see  figure  Al)  into  tlrree  segments,  SC.  GH,  and  HR,  and 
thus  construct  the  perturbation  due  lo  ellipticity  by  combining  three  contributions.  The 
first  and  third  path  segments  are  "ray”  paths  and  may  be  treated  as  such  usfn  le  results 
of  Dziewonski  and  Gilbert  (7976),  The  middle  segment,  CH,  encounters  a  slowness  equal 
to  the  slovsuess  at  the  top  side  of  the  boundary  in  the  spherical  reference  Earth  model 
used.  It  therefore  suffers  no  effect  other  than  that  due  to  change  of  arc  length  along  the 
boundary.  This  in  turn  is  second  order  in  e,  the  Earth’s  ellipticity  of  figure,  and  may  be 
ignored.  We  are  than  left  with  the  "ray”  path  effects  along  segments  SC  and  HR. 

Before  commencing  the  derivation  from  the  formulation  of  Dziewonski  and  Gilbert 
(1976)  it  is  convenient  to  simplify'  their  result  {eq.  1.3)  by  defining: 
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f(r)  =  — 

Ps  dr 


Af,  =  -  ^  e 


/  -j  ^  \  1/2 


AF,  =  ±  je[^(r7z  -  pg-) 


depending  upon  the  nature  of  the  interaction  with  the  particular  boundary'  in  the  Earth, 
which  here  is  labelled  by  i.  We  then  have  a  more  compact  expression  for  the  coefficients 


o-,^(Zs,A)  =  d6f(r)A,„{6)  +  yGF,\,n(0,)  (A3) 

Jo 

We  can  now  follow  Dzievvonski  and  Gilbert  (1970  and  write  for  the  diffracted  phase: 
+  3}  —  f  d0/ (r) A»i  (0) -f  AF|Am(0i) 

■'o  .Tc 

+  d0/(r)A„,(0)  -r  y  Af,Am(0i)  (A4) 

where  A„  -  A.,  -  A|.  is  the  epicentral  distance  of  the  grazing  ray.  The  sums  labelled  5G 
and  HR  represent  sums  over  all  discontinuities  In  velocity  (slowness)  encountered  by  the 
SG  and  HR  paths  respectively.  The  0/  under  Zmt  refers  to  the  angular  distance  from  S 
along  the  HR  path  segment.  Note  that  the  HR  path  is  a  replica  of  the  CP  path  shifted  by 
(5  in  angular  distance  from  S.  A  simple  change  of  variable  in  the  last  two  terms  of  (A4)  to 
P'  :■  (1  _  1^  transforms  HP  to  CP  and  vv  hr'.'": 
rX, 

(j;' (-.(,  A.,  +  bj  =•■  d0/(r)A,n(0)  C  ^  AF,A,m(0i) 

,fc; 

+  [  d07(r)A,„(0'  -  b)  V  AF,A„a(V  b')  (A3) 

Now  recall  that  A„AO)  ^  -'jP2,,nicos$).  where  Pz.m  are  associated  Legendre  functions 
of  degree  3.  Simple  algebra  leads  to: 

A.-aO  ‘  bl  =  A(il flj  +  smb(sin b cos 20  +  cosb sin 20) 


A]  (()  +  b )  =  Ai  (0)  +  sin  b( sin b  sin 20  -  cos  b  cos  20) 
v3 

A^iO  +  b)  =  Aj(0)  — \=  sin btsinb  cos20  -  cos b  sin 20) 

V3 

Thus,  if  we  define; 


d0/(r)cos20+  ^AFiCos20, 

•t  .iu  j 


y=  d0/(r)sin20+  J^AF,  sjn20, 

i.GP 
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we  can  write 

a^{z,,Ag  +  5)  =  pd0/(r)Ao(0)  +  X 

I.SG 

+  pde/(r)Ao(e)  +  X  AF,Ao(d,) 

+  sin6(A’sin5  +  i'cosd) 

d6);(r)Ao(d)  +  X  AF.Ao(0,) 

I.SP 

+  sin  (5(  A'  sin  <5  4-  Y  cos  <5 ) 

=  CTy  {Zj,  A(j)  +  sindiA'sind  +  V'ros(5) 


and  similarly  for  o  f  and  a-f: 

2 

(xf  ( Zj ,  A ,,  +  (5 )  =  (jf  ( z.t ,  All )  +  sin  d  d’  sin  5  -  A'  cos  5 ) 

(r‘,*(z.,A,,  -r  A)  ~  cj1'(z,,Aj)  — L  siniSlA'sind -r  I'cosd) 

^  3 


(A8) 


(A9) 

(AlO) 


Note  that  relative  to  the  grazing  ray,  the  change  in  the  cTm  tables  is  expandable  in  terms 
of  sin  (5,  where  d  is  the  arc  length  of  diffraction.  Thus,  for  a  short  diffraction  distance 
the  o'ln  tables  for  the  diffracted  phase  are  well  represented  by  a  linear  extrapolation  of 
me  tables  tor  the  associated  geometrical  ray.  It  is  apparent  in  figure  1  that  the  ellipticity 
coefficients  are  continuous  to  first  order  across  the  transition  from  an  optical  path  to  n 
diffracted  path  (dashed  vertical  lines).  It  is  possible  to  show  analyticall>'  that  this  is  the 
case  for  surface  focus  paths  (S  at  the  surface)  and  approximately  so  in  general.  The  proof 
of  this  statement  is  beyond  the  scope  of  this  paper.  This  is  potentially  useful  because  if 
one  has  the  coefficients  sufficiently  accurately  tabulated  for  the  refracted  phase  one  can 
estimate  the  integrals  X  and  Y  from  their  derivatives. 

It  is  worth  noting  that  wtiile  figure  A1  draws  the  diffraction  boundan'  near  the 
core-rnantle  boundary  and  our  numerical  examples  are  for  diffractions  around  the 
core-mantle  boundar>',  the  arguments  and  deri\  ation  above  apply  to  diffractions  around 
arbitrar)’  spherical  boundaries  in  the  Earth.  Our  derivation  holds  true  whenever  the 
predicted  time  of  the  diffracted  phase  is  constructed  by  a  path  integral  akin  to  the  path 
depicted  in  figure  Al.  Th  derivation  applies  e.g.  to  diffractions  off  the  caustics  of  upper 
mantle  triplications. 

It  is  also  worth  noting  that  our  description  of  the  diffracted  path  involves  a 
high-frequency  assumption.  High-frequency  diffractions  decay  rapidly  with  propagation 
distance  and  the  travel  time  of  diffracted  waves  is  dispersive.  The  change  of  travel  time 
with  frequency  can  be  understood  in  terms  of  a  change  in  the  diffraction  geometry. 
We  have  not  explored  the  potential  delations  from  the  above  derivation  due  to  finite 
frequency. 


JtJ 


Jt. 


I 


» 
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A2  Composite  phases 

If  we  can  construct  a  wave  path  of  a  particular  phase  as  a  combination  of  rwo  or  more 
elementary  phases  we  can  evaluate  the  ellipticiti’  correction  for  the  composite  phase 
by  adding  the  corrections  for  the  elementary-  phases  according  to  eq.  1.2  using  the 
appropriate  epicentral  distances,  colatitudes  and  azimuths.  Because  the  form  of  the 
ellipticiry  corrections  is  the  same  for  all  phases  we  can  go  further  and  construct  the 
tables  of  ellipticity  coefficients,  o',n,  for  the  composite  phase  from  the  a,,,  coefficients  of 
the  elementary  phases.  We  start  by  introducing  the  shorthand: 

cjm (0,  =  Pe,),i(cos 0) cos  mt.  (All) 

We  can  then  uTite  the  ellipticiry  correction  for  the  composite  phase,  5t‘^.  either  in  terms 
of  the  coefhcients  for  the  composite  phase,  a^,  or  in  terms  of  the  coefficients  for  the 
elementary  phases,  all,  and  a}l,  (see  figure  A2  for  definitions  of  geometry): 

6t''  =  (fo(Os,^s)a,^(Zs,A)  +  gi(0sXs)o'i(Zs.^)  +  g2<&sXs)(r>  (Zs,  A) 

=  gn( 0, ,  ( - s .  1  ( (Zi,  A„ )  +  cy2  (0s .  Cs ) erf'  (z,- ,  A„ ) 

-  go +gi( 0,, ,  C,, ) al’ (0,  A,. )  +  ( dp .  C,. ) nj’ (0,  A;, ) ,  ( A 1 2 ) 

In  order  to  construct  the  coefficients  all,  we  need  to  tran:  scribe  the  terms  involving  O,, 
and  Cp,  wliich  arise  through  the  sliift  of  the  origin  of  the  path  for  phase  b,  in  terms  of 
0,  and  IVe  utilize  the  following  relatinnsliips  from  spherical  trigonometry  (we  apply 
the  sine  and  cosine  rules  to  the  left-hand  ;.ide  triangle  in  figure  A2): 

sitic,,  suit';,  sin(TT  -Cp)/sln0j  =  sinCp/sind, 

cos  0/,  ■-  cos  0,  cos  -1  Os  sin  A.i  co-;  Cs 

cos  0s  --  cos  Of,  t:os  A,,  i  sin  0,,  sin  A,|  cos 
from  which  we  arrive  at 
sin  Of,  sin  Cp  =  sin  d,-  sin  ^s 

sin  Of,  cos  Cp  =  sin  A„  cos  dj  -  cos  sin  0.  i:c  •  Cs 

by  simple  algebraic  manipulation.  vVe  can  vi.  '.v  rewrite  the  terms  of  eq  A I  )  which  int  oh  e 
Of,  and  Cf,,  the  terms  C7m(dp,  ^p),  as: 

go(Of,,^f,}  =  -  (1  +  3cos2A„)t7o(d5,CW  +  sin2A„e7i  (dj,  Cv)  +  —  sin"  Art£/^>( d^ ,  ), 

0i(dp,Cp)  -  sin2A„tto(di,C5) cos2A„^!((9i,^5)  +  istn2ApCt2(0j,  Ci),  (Al.l) 

^zldp.Cp)  =  —  Sim  Aa5o(ds,Cj)  -  isinZAa^tiidi.^j)  +  ^{1  +  cos"  Aa)ry.;(dj,  Cs)- 

We  find  that  they  decompose  into  combinations  of  terms  gmiOs.^s)  witii  weights  which 
depend  upon  the  angular  distance  of  the  sliift  cf  origin  of  phase  b  from  that  of  phase  c, 
A,). 
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It  is  now  straight  forward  to  match  terms  ineq  All  and  write  the  ellipticity  correction 
coefficients  for  phase  c,  cr^,  in  terms  of  the  coefficients  for  the  elementary  phases: 

afiZs.A)  =  A„)  +  ^(1  +  3cos2Ai,)cro  (O.Aj,) 

/o  'T 

-  ~  sin2A„af(0,  Ai,)  +  ^  sim  AatrJ'iO,  At,) 
cjf(Zi,A)  =  cr|'(Zj,  A.,)  +  “  sin2Aa)tro  (0,  Af,) 

+  cos2AaCrf  (O.At)  -  -  sin2A,,crj'(0,  Aj,) 

o-',  {Zs,A)  -=  o-;;‘(Zi,  A,j)  +  ^  sin- A«<Jo  (0,  Ai,) 

+  isin2Ai,(Tf(0,Aj,)  +  i(l  +  cos"  A.Jo-|'(0.  Aj,)  (Al-4) 

In  order  to  iniplenienl  these  expressions  one  needs  to  match  the  ray  parameters  of  phases 
a,b,  and  c  in  order  to  determine  A,,  and  Af,  from  A  and  then  the  o,„  tables  of  phases  a 
and  b  can  be  combined  to  calculate  the  tables  for  phase  c. 

APPENDIX  B;  EllcCTRONIC  fOP.MS  FOR  ELLIPTICITY  CORRECTIONS 

A  Fortran  subroutine  el  lip  provides  access  to  the  ellipticity  corrections  lor  ;i  given 
source  location,  epiccntral  distance  and  azimuth  using  direct  access  tables  of  ilie  o,,, 
coefficients  (as  specified  in  Table  1.  This  subroutine  and  the  necessary  direct  access 
table  el  1  i  pcor .  tbl  are  available  using  anonymous  ftp  at  the  follvving  sites: 

•  the  Research  School  of  Earth  Sciences,  Australian  National  Uni\’ersii\',  Canberra, 

Australia:  rses  .anu.cdu.au 

•  the  IRIS  Data  Management  Centre  in  Seattle,  Washington,  i;,S,.\. 

iris.waslrington.edu  I 

•  the  National  Earthquake  Information  Centre,  U.S.  Cieologicitl  Surv’ey,  tiolden, 

Colorado,  U.S.  A.:  gl  dfs  .  cr .  usgs  ,gov,  directoiy  ell  ip 

A  further  program  ttimel  provides  combines  access  to  rapid  calculation  of  seismic 
tra\  el  times,  using  the  akl35  model  of  Kennctt,  Engdahl  &  Buland  (  U)');i),  with  cdlipticity  * 

collections  for  the  full  range  of  phases  in  table  1. 
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Figure  1.  (a)  (A)  elllpticity  coefficients  for  P  and  Pdiff.  (b)  a,„  (A)  ellipticity  coefficients 

for  S  and  /idiff. 

Figure  2.  Addition  of  ray  legs  to  produce  a  composite  phase 

Figure  3.  Addition  of  ray  legs  to  construct  PKKP  from  tw'o  PKP  segments  and  PcP. 

Figure  Al.  Ellipticity  corrections  for  diffracted  phases 

Figure  A2.  Geometric  relations  on  a  sphere  for  phase  addition 
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Table  L  EUiptidty  corrections 
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Table  2.  EUipticiIy  coefbtients  /'jkj 
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Table  3,  LLlipticity  coefficients  for  S^iff 

Ucptli  of  source  |knil 


0. 

35. 

50. 

100. 

A 

c^O 

rci 

(72 

<70 

(71 

(72 

(70 

<71 

<72 

<70 

(7) 

(72 

100.0 

■1.041 

0.902 

-1.480 

•1.010 

0.896 

-1.480 

•0.997 

0,894 

-1.481 

•0.933 

0.886 

■1.481 

105.0 

-1.1C7 

1.14,1 

■1.407 

-1.137 

1.138 

•1.407 

•1.124 

1.136 

-1.407 

-1.080 

1,128 

■1.407 

110.0 

-1.329 

1.357 

•1.314 

■1.298 

1.351 

■1.314 

■1.285 

1.349 

■1.314 

-1.241 

1.341 

•1.314 

115.0 

-1.520 

1.534 

■1.204 

■1.489 

1.529 

-1.204 

-1,476 

1  526 

■1  204 

■1.432 

1.518 

■1.204 

UO.O 

-1.734 

l.ti71 

■1.080 

1,704 

1.66', 

■1.080 

■1.691 

1.662 

■1.080 

-1.647 

1.(,:.5 

■l.OSO 

125.0 

-1.900 

1.7G2 

■0,94G 

■1.936 

1.756 

■0.946 

■1.923 

1.754 

■0.946 

-1.879 

1.746 

-0.946 

KiO.O 

-2. 300 

1.805 

•0.806 

■2,178 

1.799 

■0.806 

•2.165 

1.797 

-0.806 

-2.121 

1.789 

■0.806 

15  5.0 

-2.453 

1.799 

-0,064 

■2,423 

1.794 

■0.665 

■2.410 

1.791 

■0.GG5 

■2.366 

1.783 

■0.665 

140.0 

■2.G94 

1.744 

■0.526 

•2,663 

1.739 

•0.526 

2.650 

1.736 

-0.526 

-2.606 

1.728 

•0.526 

145.0 

2.922 

1.042 

■0.394 

■2.892 

1.636 

■0.394 

■2.879 

1.634 

■0.394 

-2,835 

1.626 

-0.394 

150.0 

■3.132 

1,496 

■0.273 

3.101 

1.4  (JO 

0.273 

3.088 

1.487 

0.273 

-3.044 

1.479 

-0.273 

200. 

300. 

500. 

700. 

s 

(Jy, 

<ri 

(72 

(7(1 

<7l 

<72 

('« 

(7l 

o> 

(7(1 

(7l 

<72 

100.0 

0  870 

O.HK2 

•1.482 

■0.796 

0,860 

-1.482 

•0.652 

0.824 

■1.483 

-0.528 

0.808 

-1.484 

10  5.0 

■  1.00  1 

1  ,r24 

-1.408 

-0.<123 

1,102 

-1.40!) 

■0.778 

1.06,5 

1.409 

-0.655 

1 .050 

■1.411 

1  lO.U 

-1.105 

1,337 

■1.315 

■1.085 

1.314 

■I.:!i5 

■0.<J40 

1.278 

■1.316 

■0.816 

1.263 

■1.318 

11.5.0 

1.350 

1.515 

■1.205 

■  1 .2  75 

1  492 

■1.20.5 

1.131 

1.456 

■1.206 

-1,007 

1 .44  1 

■1,207 

120.0 

1  5  7() 

1,0:.  1 

•l.()«l 

■1.490 

1.628 

-1,081 

1.592 

•1.082 

-1,22.^ 

1.577 

■1,084 

12  5.0 

■1.803 

1,742 

■0.94  7 

■1.722 

1,71') 

-0.947 

■1.577 

l.(,H3 

•0.948 

-1.4', 4 

1 .668 

■0.950 

I.Ul.O 

■2.045 

1,785 

•0.807 

l.')64 

1.763 

-0,807 

■I.H'20 

1.727 

■0.808 

■1.696 

1.712 

■O.HIO 

15  5.0 

■3,380 

1 .780 

().(i(l(, 

2.209 

1.757 

■0.666 

•2.06-I 

1.721 

■0.667 

■1.941 

1.706 

■0.668 

140.0 

3.530 

1,72:. 

■0.52'7 

■2.450 

1.702 

■0.527 

■2.30,5 

1.666 

■0.528 

■2.181 

1.651 

■0,530 

14  5.0 

3. 7.)  8 

1.(.22 

•0.,l')5 

■2.678 

l.liOO 

•0.395 

■2.533 

1.564 

■0..196 

■2.410 

1,548 

■0,;i9(! 

150.0 

•2,908 

1,470 

-0.274 

■2.888 

1,4.53 

•0.274 

■2.74,1 

1.417 

•0.2  75 

■2.619 

1 .4(1:.' 

•0,277 

Frequency- dependent  effects  on  travel  times 
and  waveforms  of  long-period  S  waves: 
implications  for  the  scale  of  mantle  heterogeneity 

Hoiiicr  Igcl.  Institute  cf  Theoretical  Geoiiliysics.  Canibi iclgc.  CK 
Oli  Gudniuiidssoii,  Research  School  of  Earth  Sciences,  Canberra.  Australia 


Submitted  to  ['EPI.  ,luly  lOhli 


Abstract 

The  sjiertral  proper  iies  of  upner-mantle  velocity  perturbations  arc 
controver.sial  While  global  .Earth  models  obtained  by  tomographic 
inversion  are  characterized  by  the  dominance  of  long-wavelength  fea- 
tMi'i's,  results  fro:;i  regional  tomography,  and  stochastic  analy.sis  of 
travel  times  suggest  the  presence  of  mere  power  at  intermediate  scales 
than  p''eflirted  by  global  models.  We  study  frequency-depcndonl  ef¬ 
fects  of  long-period  5-wa,ve  propagation  through  random  upper-mantle 
models  with  specified  spectral  properties.  Wave  propagation  i.s  simu¬ 
lated  by  a  finite-difference  approximation  to  the  a.xi-symmei rie  wave 
equaticMi  in  spherical  coordinates.  For  global  models  will)  cylindri¬ 
cal  symmetry  and  constant  angular  incretnen;  the  u,so  of  spherical 
enordiriates  leads  to  an  etlectivc  lateral  grid  spacing  (arc  length)  de¬ 
creasing  with  depth.  I'his  is  contrary  to  the  o  moiits  of  gloltal 
module  with  low  velocities  at  the  top  of  the  maiule.  wliich  mwessi- 
tate  a  dense  grid  spacing  at  small  depths  and  a  wider  grid  spacing  at 
the  ba.se  of  the  mantle.  We  introduce  a  grid  with  deplh-depeiulent 
lateral  grid  spacing  to  overcome  this  inconsistency.  Our  .simulations 
suggest  that  (1)  the  properties  of  power  spectra  of  travel-lime  i-i'sidt;als 
are  frequency  dependent;  (2)  power  spectra  of  models  oblaiued  from 
long-period  tomography  may  considerably  underestimate  the  power  at 
intermediate  scales;  (.3)  fretjueacy-dependent  effects  on  the  waveform 
are  sensitive  to  the  sca.les  and  amplitudes  of  perturbations  present  in 
the  upper  mantle. 


Keywords:  nmntlr  structure,  finite  difference.s.  .b'-waves 


Introduction 


Determining  the  spectrum  of  heterogeneities  in  tlie  EtUlh’s  iiiiuillc  is  impor¬ 
tant  to  constrain  the  types  of  convection  present  in  the  mantle  and  to  under¬ 
stand  tile  origin  of  modeled  seisniic-velocitj'  anomalies.  .Spectral  pi'operties 
of  mantle  structure  can  be  estimated  from  three-dimensional  (.'3D)  models  ob¬ 
tained  by  global  tomography  or  they  can  be  estimated  directly  from  observed 
travel-time  residuals.  A  detailed  comparison  of  power  spectra  for  different 
global  and  regional  models  can  be  found  in  Pa.ssier  and  Snieder  (1E)95).  A 
di.scus.sioii  of  global  mantle  models  is  given  in  Ritzwoller  and  Lavely  (1995). 

The  spectral  properties  of  mantle- velocity  perturbations  determined  by 
various  methods  often  differ  substantially  (Passicr  and  Snieder.  1995).  While 
tomographic:  images  obtained  from  global-data  inversion  suggest  a  sharp  on¬ 
set  of  spectral  decay  at  harmonic  degree  I  8  (e.g.  Su  and  Dziewonski, 
199'2;  Zhang  and  Tanimoto,  1991,  1993:  Pollitz.  1994)  results  from  regional 
loniograjiliy  (e.g.  Snieder,  19S8;  Zielhuis  and  N'olet.  1994;  Spakmaii  c't  ak, 
1993)  and  stochastic  analysi.s  of  travel  times  (Gudinundsson  et  ak.  1990; 
Davies  el  ak,  1992)  indicate  tliat  there  is  more  power  in  intermediate-scale 
heterogeneity  tha.ii  predicted  by  global  inversions.  'I'liese  difiereuces  arc  at- 
tributcci  to  (1)  varic'd  effects  of  the  inversion  sclicmes  (e.g.  I..aske.  Mastens 
and  Zuni.  1991;  (2)  complexity  of  the  regions  wliere  regional  tomography  is 
carried  out  (Passierand  Snieder,  1995):  (3)  filtering  effects  of  finilc-freciiu'iicy 
piupagatioii  '..ii  l.avcl  limes  (Gudinundsson  and  Cummins,  1991). 

'I'o  understand  frequency-dependent  cfTi'cls  on  tlie  del  ('rminal  ion  of  spc.’c- 
tral  propel  ties  we  perform  full  wavefield  calculations  for  mantle'  models  witli 
known  spectral  iiroperties.  The  numerical  method  we  u.se  is  an  e.xlc'nsion 
of  the  finite-difference  (ED)  approach  by  Igcl  and  Welx'r  (1995.  ItJfXi)  to 
grids  witli  depth-de[)endent,  lateral  grid  si)acing.  Grid  refineiiK.'iit  near  tlic' 
Earlli's  surface  considerably  improves  the  performance  of  tlie  ED  algorithm 
because  (1)  low  velocities  (small  wavelengtlis)  at  tlie  top  of  tin'  mantle  can 
be  sam[)led  more  efTiciently,  and  (2)  liighcr  veloi  ilie.-  in  ,  .le  i..A\ei'  mantle  are 
sanijiled  less  densily,  improving  .str- u.lity. 

lire  purpose  of  this  [lapcr  is  to  present  an  ED  algorithm  will:  depth- 
dependent,  lateral  grid  spacing  for  global,  toroidal  wave  [)ropagatiun,  to  dis- 
cu.ss  its  accuracy,  and  to  simulate  frequency-dependent  waveform  effects  for 
uppei mantle  models  with  specified  spectral  properties. 


Numerical  algorithm 

The  algorithm  presented  in  this  paper  is  an  extension  of  ilu’  high-order  FU 
scheme  of  Igel  and  Weber  (1995)  to  grids  with  vertically  varying  grid  spac¬ 
ing.  The  equations  of  motion  in  sj)herifal  eoordinate.s  r.  arc'  soivc'd  for 
wa\’efielcls  and  media  invariant  in  y-,  thus  rotalionall\'  symmetric  with  re 
spec't,  to  tlur  axis  0  —  0.  For  spherical  or  polar  grids  with  ecpial  spacing  in 
the  angular  domain,  the  grid  is  coarse  for  large  radii  and  dense  icear  the  cen¬ 
tre.  T  his  is  coiitrary  to  what  is  required  for  global,  cartesian  Farth  models. 
The  velocities  near  the  Earth's  surface  arc  considerably  lower  than  deeper 
within  the  Earth  and  therefore  require  dense  rather  than  coarse  sampling. 
In  addition,  heterogeneities  are  clearly  strongest  in  the  lithosphere  and  up¬ 
permost  mantle  (]3ossibly  with  the  exception  of  the  lowermost  mantle.  D") 
which  conseciiK’nlly  should  be  sampled  with  a  dense  grid. 


Theory 

.'Vssuming  iin'ariancc  in  y-  the  cquaiions  for  loroichd  motion  ic,.  arc  Laj)- 
wood  and  Usami,  19.SJ) 
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f)  I -I  tlic;*  ma.s.s  density,  id  a  VuEiuc  .w.  *.  dO  the  c  i>i,ijKnii.‘nts 

of  the  stre.ss  tensor. 
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II  being  tlu'  shear  modulus. 


Discretization 

All  space  dependent  fields  -  stress,  sources,  and  displacc’nient  arc'  dc'lliu'd  on 
a  spherical,  staggered  grid.  .Staggered  grids  hn\  c'  been  used  wich'ly  in  isotro[)ic 
(e.g.  Virieux,  1984;  198(1)  and  anisotropic  (Igel  ct  al,,  199.'))  I'D  wavefiedd 
calculations.  In  staggered  gricls,  the  elements  of  disi>lacemenl .  slre.ss.  and 
strain  are  not  defined  at  the  same  locations,  allowing  the  first  derivati ve.s  of 
those  fields  to  be  centered  in  between  grid  points.  i)no  to  the  anlisymrnetrj' 
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of  t.he  clifforence  operator,  this  iin|)roves  the  arcitrary  of  the  fiiiite-dificreuce 
calculation  for  a  given  space  incn'inent. 

The-  spherical  grid  is  shown  schematically  in  Figure  1.  '1  he  grid  is  defined 
bclwccri  the  Earth's  surface  and  the  core-mantle  houndary  (CMB).  These 
two  boundaries  are  modeled  as  free  surfaces,  do  im]>leinen1  this  houndar\' 
condition,  the  grid  is  e.\tended  Lj'i  +  1  points  above  and  below  the  physi¬ 
cal  boundary,  A  being  the  length  of  the  FD  difTorential  (and  inter]'>olation) 
op(’rafor. 

In  the  toroidal  case  with  eonstant  grid  sj)acing  in  both  th<'  angular  and  the 
rarlial  direction  this  leads  to  a  consistent  scheme  without  rec|uii  ing  addit  ional 
interpolations.  Howe\'er,  to  link  the  two  domains  with  different  lateral  grid 
spacing.  intei|)olat.ions  arc  recinired  at  or  around  the  depth  of  the  grid  change. 

The  free-snrface  boundarv  condition 
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is  iin[)leiiiented  by  imposing  symmetry  of  u^.  and  antisyminelry  of  the  Pr,- 
coniponeut  with  respect  to  tlie  free  surfaco(s). 

.^s  shown  ill  Figure  1.  the  elements  defined  at  —  0  are  t  he  (lisplticcmeni 
(/„.  and  the  stntss  element  a^o-  At  0  =  U  the  ecpiatioii  of  mol  ion  is  singular 
and  the  dis|)la(:eineut  field  can  not  be  evaluated.  Ilowevcu'.  by  s\'nunetry  we 
impose  ic,.  =:  (Tro  =  0.  .'\nalogous  to  the  situation  at  the  free  surface  the  grid 
is  ('.N;tenoed  A/2  i  1  .•  Ijc-tindary  f  'J  .  :  ,1  tl,  b  aunlaiy 

condition  is  im[)osed  by  anti-symmetry  of  both  stress  and  disi)laceinent  field. 

The  source  region  near  d  =  0  is  schematically  shown  in  Figuri'  2.  Onr 
goal  is  to  model  point-source  like  behavior.  To  achieve  this,  sonicc's  ai'c  input 
at  a  grid  point  close  to  the  axis  of  symmetry.  A  directional  force  would  thus 
result  ill  a  ring  source  (pure  toroidal  motion),  '['lu'  radiation  patteini  of  .such 
a  .source  is  unrealistic.  Nevertheless,  it  allows  us  to  stiidv  moiiel-dependent 
wave  phenoiueua  vvherc  the  source  radiation  pattern  is  irrelevant. 

The  motivation  for  changing  tlie  angular  grid  spacing  was  given  aho\’e. 
.Ja.stram  and  Tessmer  (]99'1)  and  Rodrigues  (1993)  both  ini roiliict'd  FI)  grids 
with  vertically  varying  grid  .spacing,  condensing  the  (cailesian)  grid  near  the 
surface  by  a  factor  of  3.  1  he  stability  of  explicit  FD  algoriihuis  is  generally 
of  the  form  5-  const.,  where  c,„aj.  is  the  maximum  veloiily,  and  dt 

and  dx  are  lime  increment  and  grid  s]niciug.  resju'el  i  vedy.  Dejiending  on 
the  actual  velocity  mode',  the  time  slop  has  to  be  decreased  in  accordance 


willi  tlie  chango  iti  grid  spacing  to  keep  the  same  level  of  s( ahilitw  We  lind 
a  change  of  lateral  grid  spacing  by  a  factor  of  2  optimal  for  gh>l)al  models 
taking  account  of  computation  time. 

The  partial  differentials  and  the  iiilerpolations  are  calculated  by  high- 
order  operators  of  length  L.  The  weight.s  of  the  operators  were  (d)taiiu'd 
following  the  approach  of  Holberg  (1987).  The  time  evolution  is  cai  i  ied  out 
by  a  Taylor  expansion.  Details  of  this  type  of  FD  algorithm  ap])lied  to  wa\'c 
propagation  on  a  3D  cartesian  grid  can  be  found  in  Igcl  et  al.  (199b).  In  all 
simulations  carried  out  in  this  paper  we  use  8-point  space  oiierators  (/.  =  8) 
for  both  derivative  and  interpoiation.  The  t ime  extrapolation  is  accurate  to 
fourth  order. 

Comparison  with  other  methods 

Del'ore  a[)i'lying  our  algorithm  to  arbitrarily  heterogene<uis  models  wi-  loiii 
pare  FI)  synthetic  seismograms  with  tho.se  obtained  In'  an  r.iacf  mi’thod 
(Direct  Solution  Method  (DSM),  (jeller  et  ah,  1994;  Cuinmiiis  el  ah,  1091) 
for  sphci'icrdly  syumicTric  media.  The  FD  grid  is  divided  into  two  domains 
with  depth  ranges  0-320  km  and  320-2891  km,  willi  grid  sizes  20  18x01  and 
ll)28xbl2,  rcspectivels'.  The  aiigulai'  domain  is  O-tt  ami  the  tinu'  incix'iueuL 
is  0,3  seconds.  The  model  is  the  i.sotropic  part  of  PUF.M  ( Dzu'Womski  and 
Anderson.  1981).  jMd 

In  Figure  3  seismograms  are  compared  for  two  epicenlral  dislances  at  a. 
periorl  of  20  seconds.  T’he  .sourer'  (toroidal  ring  .source)  is  al  200  km  depth. 
Tl;e  overall  agrerum'iit  between  the  entirely  dilhM'C'nl  mimerical  lecImiepK'S  is 
excellent.  Small  pha.s(’  differences  are  attril)uled  lo  the  discreieiiess  of  the 
I'T)  grid,  which  can  be  suppressrrd  by  further  grid  relim'mciil s.  l  lie  relalir'e 
ami)litude  of  the  phases  as  well  as  the  waveforms  are  well  modeled  Ia'  the 
IT)  approximation. 

How  docs  the  grid  refincineiit  near  the'siirface  im|)rove  the  perlormaurr'  of 
the  FD  scheme?  To  demonstratt.*  tl:e  performancr'  of  the  suggested  algorilluii 
we  imdei'takr;  a  convergence  test.  \\'e  compare  seismograms  from  llirer'  dif 
ferent  simulatioii.s:  I.  Kefcreiicc  seismogram  obtained  witli  a  very  dense  grid 
(bl20xl024);  IT  Seismogram  for  a  1024>'!bC  grid;  Ill;  same  a.s  II,  l.nil  with 
a  2048x32  refined  grid  at  the  top  of  the  modr'I. 

The  algorithm  allows  us  to  generate  a  syiilheiie  Clreen’s  function  whicli 


wc  can  later  convoh'c  with  a  source  wavelet  of  the  desired  freciiK'iiey  baud. 
In  Figure  4  seismograms  are  shown  for  three  different  dominant  periods  at  a 
distance  of  GO"  .  At  the  c.oiisidered  frc'tiiienc.ics  the  coarse  grid  eoiiv'erges  ()nl> 
at  a  period  of  60  seconds  while  the  grid  refinement  enables  the  calculation 
of  accurate  seismograms  down  to  a  |)eriod  of  =^26  seconds.  I’hc'  incia'asc'  of 
(,'PU  time  for  the  simulation  of  type  III  with  respect  to  type  11  was  dO‘X. 
However,  this  increase  leads  to  a  substantial  improvement  in  accuracy. 

A  further  dc'inonst.ration  of  the  aecuraey  is  given  by  a  eomparisou  with 
i'a\ -theoretical  arrival  limes  (Witte  et  ah,  determiiieil  for  a  2d)  model 

with  lateral  perturbations  added  to  the  isolro|>i(:  ijart  of  PRHM.  I’lie  lateral 
■S'  velocity  perturbations  are  .shown  in  Figure  .'m.  The  i>erl  nrhations  are  co¬ 
herent  down  to  a  de])th  of  2o0km  where  they  are  tapered  to  zero.  IVavel  tiinys 
for  the  .S'  (.Saar)  pha.se  from  the  FP  seismogranis  were  obi ained  by  maximiz¬ 
ing  cross-correlat  ion  belw(H;n  perturbed  and  unperturbed  seismograms  for  a 
period  of  30  seconds.  ( krmpaiison  of  the  travel-time  ixn'Lurbation.s  obtaiiu'rl 
!'■■  ihe.se  diircrent  approaches  ar'-  shown  in  Figure  Td).  The  sampling  rate  of 
the  FI)  .seismograms  is  2/3  secoiuls  and  the  receiver  sparing  is  1/3  dt'grees. 
'1  li(‘  maxiimmi  dilferc'iice  Ix'twr'cn  the  ray-th<H)retical  and  I  I)  travel  time  is 
%1  second. 

Since'  (he  direct  S  (Sdiir)  phase  is  arriving  at  a  very  slee'j)  angle',  we'  e-xpe'e  t 
anti-ceu'!  ''' ■’* '  •"  I'-e .,.11  vrdeieSity  anil  travel-time  perlnrhal  ion.  whie  li  e'an  be' 
appre'e  iate'd  Ijy  e  ermparing  Figure's  .aa  and  .ab. 

TIu'  lu.nj  e/.  e/.ir  e.lgejri; h:i'.  a  fai-  as  wavef'.'vm''  anel  ! r.-'.'.i'!  'imes  a'.'e 
(onee'i'iie'd  give's  ns  e:ou[ide:ne'e'  I  hat  we'  e'an  apply  onr  FI)  me'lliexl  to  study 
iiiiili'.  fre'e|uency  amplitude'  anel  wav('fe)nn  e'lfe'e'ls. 

Spectra  of  mantle  models 

riie  aim  e)f  this  stuely  is  to  investigate'  the  fre'<|ne’ne'y  ele'pe'nelcnec  e)!'  liavel 
time's  anel  wave.'fonus  eif  ,S'  anel  SS  wave!s'Tor  ui>pe'r .mantle'  mexh'ls  willi  vai'- 
ieel  powe-r  s|x'e;tia  anel  the'  implications  fe)r  the'  se:ale“  ed  mantle'  hei.e'rejgene'il y. 
A  se'he’inatic  leinese-nlation  of  the-  powe-r  spe-ctrum  eif  mantle-  \ehxity  ix-r 
tiirbations  is  shown  in  Figure  6.  1  he-  twee  ty])e;s  e>f  spe-etra  slienen  in  this 

Figme'  (tvpt-  1  and  type  II)  r(!])re-'sent  extr<.-me-  e:ase's  of  Passie-r  anel  .Snie-elei 
( IDl).')).  The-y  showed  that  at  iiiternie-diate:  se;aless  (1=30)  the  spe-etra  ejljtaine-el 
by  global  and  regional  inversions  dilfe-r  by  as  much  as  a  factex'  ed  10-30.  One-  of 


their  argumenls  for  th('  prcseiu'o  of  signifiranl  inlornK'clirilc-.srale  lii'icrogeiie 
ily  i.s  the  eomirlexity  of  long- period  surtare-wavc  ,s]K'ctra  clii<’  to  mull  i|iiit  hiiig 
and  the  associated  intcrfercnctM'irocis. 

We  generate  models  witli  type  1  and  II  statistics  shown  in  I'igiiii'  (1.  l)e 
tails  of  the  model  gcmeiation  are  given  in  tin'  Ap]K'ndi\,  higiii'e  7  shows 
7)  raiKlom  model  realizations  of  each  spectral  tj'pe.  (dearly  spectral  l>pe  1 
contains  mainly  continent.-scalc  features  while  type  II  contains  a  signilicaiit 
amonni  of  energy  at  scales  down  two  100  km  wavelength. 

'lu  ti'st  the  spectral  prop<Mii('s  of  the  synthetic  models  we  aNcrage  the 
powei  spectra  of  several  realizations  for  each  type.  Nmnerical  tests  show 
that  at  least  .7  power  spectra  ol  diih'rent  realizations  have  to  he  summed  to 
achieve  an  ticeuracy  of  5%  for  the  decay  rate.  Hi's.dts  are  shown  in  I'igui'e 
,S.  '1  he  power  sjiectra.  disjday  the  characteristic  comer  frcipiency  and  the 
s|)('ciral  decay  whic  h  was  Impo.sed  in  the  gc.neration  of  the  models. 

2-1)  models  of  ,S'-velocity  perlnrhalions  are  ohta.inerl  by  supei  imposiiu'. 
se\eial  latei.d  11)  functions  woighii.'d  by  dej)tl'.-depeiideut  orilioguii.d  lime 
lions.  Details  about  the  21)  model  generation  are  given  in  the  .Ap|)eiidi\. 
Sections  of  modcl.^  thus  ohlaiinxl  are  shown  in  hij'pire  f).  1’he  .'i'  velocity  per 
I  iirbatioiis  for  all  the  simulat  ions  dc'serihed  Indow  are  tapered  down  to  zero 
In'tweeii  .htltlkiii  and  (iOOkm  depth.  ’!  here  are  no  lateral  velocity  pertiirha- 
tions  below  (iOUkm  depth.  It  is  lm|iortani  to  note  that,  the  2- 1)  perl  iirbat  ions 
vanish  iiea.r  the  a:\is  of  rcjtational  symmetry  0  —  0  t(»  avoid  focusing  edects. 
'I  lie  perturbations  are  iiii|)Ose<l  for  II  >11)"  . 

Ill  the  following  we  iii\'esligate  sviithelie  seismograms  oblaiiied  for  tlie 
model  types  just  described.  All  models  have  2  !)  .S' velocil y  ix'rl urbatioiis 
according  to  spect  ral  types  1  or  11  from  f  igurc'  (i.  '1  he  iiia.xiiinim  amplil  iide  of 
the  ,S- velocity  perturbations  varies  between  0%  and  12%  (root  inean  scpiare 
s'ariat.ioii  between  0  and  11.2%). 

Syuthtitic  sttisinograma 

W  e  (-,.1 .  nlati  at  hel  ic  seismograms  for  pertiirluilioiis  of  types  1  and  II  adch’d 
to  the  isotropic  jiart  of  I’Rh'.M  with  the  following  set.iip.  flic-  mantle  is 
samph’d  in  dc-pth  by  578  grid  points,  leading  to  a.  dept  li  spacing  of  5,t)km.  The 
np|)er  part,  of  t.he  mantle  (0-32l)km  dept  li)  is  dehiied  cni  a  grid  of  size  25111)  xti  l, 
w  Idle  the  bottom  ])art  (320-2891  km  depth)  is  deliiied  on  a  1280x512  gricl.  W'e 
focus  on  ,S  and  SS  phases,  which  are  fre(|uently  used  in  global  tomographic 


si  lulic's.  Scistnograiiis  of  l('iip;l  li  ]  lunn-  ( 12001)  I  iiiic’  sl('|).s  with  I  imc  im  i'(‘im;nt 
(It  =  0,3.s)  aio  ralcnlaU’d  foi-  cijiroiitral  clistaiu'os  50  150"  with  a  rcccivoi 
sparing  of  l/.'l"  .  Windows  containing  .5'  and  5', S'  arrivals  arc  ('Nlractcd  using 
rav  lh<'orcl.ical  travel  t.iiiK's.  A  sinnilalion  wil.li  the  iiarainctcrs  di'scrihcd 
al)ov('  takes  ~15  minutes  CIM  time  on  a  r2S-proc<‘.ssur  (tM5. 

Seisinugranis  for  S'aiir  «'  I'eceiver  sampling  of  1  degree  ai'c  shown  in 
h'igiire  10  lor  a  doniinunt  period  of  25  seconds  for  epiceniral  distances  100 
150"  .  In  all  simulations  the  wavehu  is  the  first  deri\'ati\'e  of  a  (iaiissian. 
The  seismograms  shown  in  I'  igure  10  were  obtained  for  a  2-1)  realization  of 
spectral  ty|)es  I  (top)  and  II  (bottom)  for  dilfereiit  niaxiimim  perturbation 
amplitudes.  In  t  he  contc'xt  of  seismic-velocity  aiiomtdies  ami  their  relation  to 
tectonic  features  the  small  and  large  perturbation  amplitudes  mr.y  repre.seni 
icgioiis  of  litth'  (('.g.  cra.tc'iis)  and  high  (e.g.  island-arc  regions)  tectonic 
acii\'ity. 

Coucerning  the  A'aiii' pha.se  we  ciui  make  the  following  (d)servat  it)ns:  (1) 

’  ill  model  type  1  till'  A'  wavi'form  i-  essentially  undisturbed  even  lor  large 
perturbations;  (2)  for  model  ty|)e  II  th<‘  wavi'form  begins  to  be  alfecled  bj’ 
t  he  .si'at  tering  loi'  jieii  ui  bat  ion  aniplit  mb's  g;realer  than  ■I’/. 

,  I'loin  the  .same  set.  ol  sei.siiiogratns  we  extrael  time  windows  coni, tuning 
t  he  phtise.  'I  hese  tin'  shown  in  h'ignre  1  I .  for  modi'l  1_\  pe  1  similar  ohsi'r- 
s'iitiuns  Ciui  be  made  tus  lor  the  .S'  plituse.  'I'lie  wtivefonn  in  the  splierica.lly 
s.vimnetric  c.-ise  tlie  Hilbert  transform  of  tlie  S  phtise  (('boy  tmd  Hiclitirds, 
Il)’i5)  is  litirdls'  allectod  by  the  velocity  pertni'litil  ioin;  eien  at  Luge  pertiir 
bation  iunplitiides  (l2'/().  lloweN'cr.  the  .S.S-wtiveforins  of  model  i>'pe  11  tire 
se\'erely  dislol'ted  for  pei'lurbtd ions  huger  than  •!'/.  It  is  worth  iiolin,!!;  tlnd 
most  of  the  ,S,S  coda  in  h'ijpire  11  (huitoin)  is  precnr.sory. 

How  due.s  the  wtivefonn  distortion  depeml  on  the  doiniiitiiit  period  of  the 
was'elel  used'.'’  h'or  long-ixn  iod  wavelets  (7.5  seconds)  the  sit  util  ion  is  (piit  e 
diilerent  (see  figure  12).  Although  there  still  is  ti  considri tilde  dilfercnce 
between  the  two  inodid  l\'|)es.  the  wavel'orius  for  model  type  II  tire  inucli  le.ss 
(listoi  ted  t  lian  iit  <i  domiiiiiiit  pi'riod  ol  25  seconds.  We  1 1 w  to  (|ii;iiilily  l  lii.s 
|;i'(|iien(  y  ilependent  SS  wtiveform  dislortion  in  the  next  section, 

I'YcMiueiicy-dopoiideiit  oil'i'chs 

lo  describe  the  lre(|ilency  dependenl  .'.ctit  lering  edect  we  process  the  SS 
windows  in  the  folluwiiii’,  wtiy:  (1)  hind  the  Lime  delay  r  of  the  perturhed 


phase  willi  i(’sj)er1.  to  a  refeixjiicc  (PREM)  wavelet.  (2)  Ceircct  the  time 
delay  and  ealeidate  a  seattci'ing  index  .s(7’)  (iiormali'/ed  root -niean-squme 
dispacemeiit  misfit  in  a  time  window  of  length  27',  T  being  the  dominant 
pei'iod)  aeconling  to 
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I'lii.s  piocessiiLg  is  carried  out  for  the  seismograms  shown  in  Eiguri's  1  1 
and  12.  for  periods  from  20-110  seconds  and  00  receivers  between  .jO"  and 
lotV  epiceiUral  distance.  I’Ik'  results  aix'  shown  in  Eigiire  l-'l.  for  spectral 
ly|)e  I  (Eiguri'  lii,  left)  even  for  large  inode!  perturbations  (12'X)  the  aver- 
ag('  waveform  distortion  does  not  exceed  50%.  In  the  case  of  s]K'clral  tyix' 
11  (i'igure  lii,  right)  there  is  a  rapid  iiicrea.se  in  waveform  distortion  at  a 
period  of  wiiU  seconds.  I’lie  eouvergence  of  the  jicrturlx-d  and  impei  1  in  bed 
wavi'fonns  is  more  rajiid  Idr  models  of  type  !. 

\','e  now  low!;  at  the  I'reiinency  dependence  of  the  speclrmii  of  .s', s'  iia\-el 
lime  ix'il nrbal ions  measured  at  tlu‘  surface.  The  model  considered  is  ol 
speciial  Ivpi'  II  (see  I'igure  (i).  'Idle  tnaximum  ])erturl)al ion  ampliliidc  is 
(i'X  .  wliicli  is  a  likely  uiiderestiinate  of  I  lie  streiiglli  of  lielcrogcneil y  in  I  lie 
n|iper  maul  le  judged  liy  results  of  surface- wav<'  tomography  (e.g.  '/iclliiiis 
and  Nolet,  llliJ  1). 

Till'  dense  and  regular  receivi'r  saiiipling  (l/ii"  Ix'tweeii  50"  and  150"  ) 
allow.s  us  I o  achieve  high  accuiacy  in  tlic  tras'cl- 1 ime  siiecl  ra  at  intermediate 
scale.s  (I':  iiOO).  d  lie  travel  limes  are  determined  l>y  iiiaximi/ing  the  c  ros.s 
correlation  lielween  perturbed  and  unperturbed  .S'.s' wavefoi'im-.  I'raN'el  limes 
are  iiicki'd  from  .seismograms  of  three  dilferiMil  domiiiaiil  periods  (25.  50,  and 
1011  seconds),  I'lxaiuples  are  shown  in  Eigiiri'  M.  Wliih'  lliere  i,s  oi’erall  a 
good  correlation  bi'tweeii  the  travel  time-  residuals,  they  lend  to  be  larger 
and  show  more  details  at  shorter  pi'riods. 

To  obi  .-un  si  able  result  s  for  the  Irai'el  time  sped  ra  we  sum  power  sped  i  a 
I'roin  .5  (liH'erent  realizations  for  each  of  the  three'  I'reepieiicy  hands  consid¬ 
ered.  d  he  decay  rate  of  the  stacked  power  S|)ccl,a  is  deleniiilied  b\’  lillnig, 
a  liiK'  through  tli<’  spectrum  between  liariuouic  degi'ce.s  /  -  211  and  /  -  dlKJ. 
Tli<'  results  shown  in  Figure  15  demonstrate  a  clear  ireipiem  y  dependence. 
At  high  frecpieiicies  (dominant  period  25  seconds),  tlie  traiel  time  residuals 
rellect  well  the  iiiojiei  ties  of  tlii'  siiedrimi  of  the  uiideriyiug,  \elucily  |iei  lm- 
bations.  However,  already  at  50  seconds  dominant  [leriod  the  slope  of  the 


(Ic'cay  is  nearly  2  powers  difrereiil  from  the  one  in  the  model  (see  figure  S). 
Furl hermore,  the  spectral  corner  is  shifted  to  lower  harmonic  degrees,  which 
i.s  even  uicu'e  sigiulicant  at  a  dominant  pr-riod  of  100  seconds. 

At  I  =  100  (ssdOOkm  wavehuigth  at  the  Earth’s  surface),  the  sjx'ctral 
[lower  determined  at  periods  25  and  50  seconds  varies  by  a  factor  of  RilO, 
This  clearly  indicates  that  for  the  model  considered  in  this  e.sainple.  the 
S|)eetra  ot  .S' .S’  travel-time  residuals  relh'ct  the  actual  sjiec.tral  [iropei-ties  of 
t.he  ui)j)er  m.intle  only  at  [uu  iods  close  to  or  smaller  than  '25  seconds. 

Discussion  and  Conclusions 

Tilt'  goal  ol  this  study  v  as  to  aualy/e  fiill-wavefield,  synthetic  ,sei.sinograms 
lor  random  up|)('r-mantle  models  with  specified  spectral  propertie.s  and  to 
investigate'  Irc'cpK'iicy-depc'ndeiit  ('Ih'cts  on  wa\'eform.s  and  travel  time's  eif  S 
and  .s', S'  wave's. 

Tlu'  s^’iit  helie'  study  we;  have;  imde'rtalven  is  statistical  in  nature  and  re- 
epiires  a  numhe-r  ei|  wave'  siinidations  le)  he  e’affie'd  out  iu  eerih-i  to  aeliiee'e'  a 
stable'  e'.st  iniaie  ol  spe'e't.ra.  'rhe't'eleire',  an  e-die  ie'nt  leerward  le’ehnie|in'  edlow- 
Ing  us  te)  meiele'l  e  eniiideete  seismograms  for  arliilrary  homejge'iie'eeus  si  rue't  lire's 
was  m'e-e'ssai'v.  .Ae'e'iirate'  5d)  ghihal  simulal ienw  witlu)ul  assiim])i ions  on  the 
amplitude' eil  [lei  l.urhatieni.s  are-  still  teie)  e'xpe'iisive  ceempuliitionally.  We  eoni 
promise'  ill  using  a  FI)  eepproximat  ii/ii  lei  the'  wave*  e'epial  ie)ii  In  spherical  e'er 
eu'dinale.s  where'  all  fie'lds  are*  invarnmt  in  y-.  'lids  h'ael.-,  Ie>  a  '!  1)  |iieeMe'in 
w  hieh  e'aii  lie'  solve'd  e'lllciently  on  pre'.se'iil- elay  paralh'l  e  einijiute'is. 

We  e'xte'iiele’d  Idle  h'l)  algorilliin  Migge'sle'el  by  Ige-I  anel  W'l'be'i'  (1111)5)  to 
grids  with  vi'rtieallv-varying,  lateral  griel  spae'ing.  'Idle'  grid  re'linemeni  imar 
the'  I-;,',  's  surlae'e'  is  ui'ceissary  freim  an  algeirilhmie'  peiiiil  eif  \'i('w  (the  ai'c 
le'iigtli  eii  the  angular  grid  inereineiit  incre'a-se's  with  distanec  from  I  lie'  e  e'ii 
l.er  eif  the  npliere')  anel  ll'om  ve'loe'ily  meiele'l  coiisiele'ratieilis  (low  ve'loeil  ie's  ill 
idle'  e'l'ust  anel  litluispliei'e;  ro(|ui;'e'  snialle'i' grid  spac  ing).  1  he  e'hang;e  in  grid 
s|)aeiiig  ha;;  been  implemented  using  a  higb-ofele;r  approxiiiial  ion,  I  he  rela 
ti\e'l,\'  thin  layer  (<  2()0km)  wil.li  relined  grid  spacing  lead.s  lo  an  aehlilional 
lit."/  in  computation  time  (as  compared  to  a  simulation  wit  hunt  reliiienieni  ). 
Ilowee'cr,  l.o  achieve;  the;  same  accuracy  without  grie.1  re'liue;me'ul .  simiilal  ioie; 
about  lour  tiiiK;.s  longer  would  luue’  heen  neres.sary.  This  i;.  erueial  lor  this 
st  udy  whe-re  the,'  investigation  of  random  structure-'s  and  llie'ir  spect  ral  prop 


oriic.s  |■<’C[^liI■os  itiaiiy  Kitiiiilnl ions  Lo  Ix'  carried  out. 

The  scale  ol' s('isniic'V(.‘lo<:ily  perturhalions  in  llie  Kiirlli's  maul  le  is  im|)oi  - 
taiit.  in  many  asijerls  of  goodynamic.s  (e.g.  inanlle  coin'ectioii.  snhduct ion. 
hot  .si)ols,  etc).  In  several  .studic's,  eouclusions  about  mantle  coiivect ion  \ver(' 
drawn  from  power  .spectra  of  global  tomographic  models  (e.g.  Monlagner 
and  Taniiuolo,  li)!)!;  Monlagm'r  or  even  dinrclly  from  the  ])ow('r  spec¬ 

tra  of  .'s',S'-li'av('lt,im<'  residuals  (('.g.  Su  am!  Dziewouski.  Iflhl.  1!)I)L’).  In 
a  la'cent  .study  I'as.sier  and  Snit'der  highlij','iled  the  disciepancy  be¬ 

tween  the  scab's  oi  velocil.y  perturbations  obiained  from  glolial  and  regit)ual 
tomography. 

Our  synt  liefic  mod('ling  suggf'sLs  I  hat,  as  far  as  S\  SS  and  .‘'  .'s'  N  si, tidies 
are  concerned  l.he  discrepauc:i('s  re])orled  by  Passier  and  Snieder  (190,"))  are 
at  least,  in  part  due  t  o  linil('-lr(!(|ii('ucy  I'llecls  ol  vvavt*  propag,al  ion.  II  t  he 
IraM'l  l.ime  perLiirbal  ions  observa'd  at  t  he  Karl.h’s  surface  do  not.  I'ellei't.  Lhe 
scale  ol  the  underlying  maiiUe,  it  i.^  mii>ossible  for  delay  lime  lomograpliy 
to  model  t  lie  powi'i  spei  tnmi  of  mantle  helerogeiieily  coirecllv.  Similarly,  if 
phase  velocities  nieasnred  along  the  |•'.,■u■lll's  surface  do  not  relieci  tin'  details 
ol  the  sti'iiclnre  ol  lhe  mantle  bt'iiealli.  tli.il  detail  ca'in  ,:  bi'  !ri()\-ered  b\' 
invei'sioii,  Waveforms  ami  amplitmies  are  necessary  to  model  IS'irlh  siriiel  ure 
in  detail. 

lhe  results  sliown  hi  figure  Id  : . ’  '  tin'  f>'.',e!eiics’ dependet|i 

elh'cls  on  the  wavelorm  a.re  st'iisilive  to  the  .scale.'^  ,ind  amplilmles  ol  per 
lurbal.ioiis  present  in  lh<'  up|)er  manile.  I  lii^  poniis  to  possible  proees.sing 
techniques  which  shoidd  be  used  in  coiubinalion  with  lrav>‘l  I imi' processing 
when  esi  imal  ing  the  power  spectrum  of  mantle  vehicity  perl  ml  ml  ions.  (1ml 
mundssoii  (l!)9(i)  studied  the  ellects  ol  lre(|uency  on  waveform  ill:, tori  ion  and 
travel  limes,  I  le  concluded  that  depending  on  I  he  baud  widi  h  waveform 
distortion  could  be  minor  while  dillraction  ellecl'i  on  limiii".  are  siguilic.ml. 
Ills  results  and  our  synthetic  modeling  suggest  that  even  though  was'eform 
dislort.iou  may  be  signilicanl,,  il  may  not  be  clear  from  a  single  seismojpam 
(e.g.  figure  1  1 ). 

The  powei  spectra  in  figure  l.a  demonstrate  that  the  characierisi  ies  of 
the  ])ower  .spectra  ol  SS  travel-time  residinds  are  strongly  Ireiiuem  y  depen 
dent.  I  he  slope  oi  decay  extracted  from  the  powt'r  siiecirum  \ari<'s  with  the 
irequency  band  at  wliich  the  |»ower  s|)ectrum  was  deteriiiiiied,  l  or  exanqile, 
al  around  /  H]()  (ssdOOkni  wavelength),  the  power  npecirum  determined  at, 

'd.u  seconds  dominant  period  contains  around  10  limes  more  energy  than  the 
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one  d<'U‘rniino(l  at  a  dominant  i)erio<l  of  r)0  seconds.  'I  liis  implies  dial  con 
elusions  on  the  scale  of  mantle  heterof;eneitie&  from  long-period  s('i, sinograms 
may  he  erroneous  if  these  linit-frcejucncy  effects  ar<'  not  taken  inl.o  aeeounl. 

Onr  results  suggyst,  that  the  travel  time  lai.siilual  speelra  lelermined  hy 
,Sn  and  n/iewonski  (lOh'J)  do  not  ridlecl  the  scale  of  mantle  hel  erogeneit  ies, 
in  part  .it.  Iea.st  due  to  liltering  effects  of  finitc-frer|nency  inopagat  ion.  and 
lh.it  tlu'  s))('clra  they  estimated  an-  ronii>a.til)le  with  mantle  models  with 
ronsiderahly  more  miergy  at  inlennediati- scales. 

It  is  important  to  note  that  since  onr  iiiodi'ls  are  iiu'ariaiil  in  y  and  thus 
two  diiiu'iisional  the  ellects  we  descrihed  in  this  |)aper  ;ire  i  oiiserv.it  ii’ely 
desc.rihed.  for  If  l.t  models  with  the  saiiie  sirectral  in'opi'rties  we  expei  l  the 
liltei'ing  ell'ects  to  In'  even  iiion'  pronounced. 
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Appendix:  Generation  of  velocdy  perturba¬ 
tions  * 

W'c  1-1)  iiiodols  with  ])o\vi‘r  sjx'Ctral  |)ri)|><’rti('.s  s.iowii  in  l'‘if!;iir('  (i. 

Uir  rc^lulioii  bcl  ween  waviMminlx'i' /r  and  liariiiDiiit'  order  /  ^  -  /  1  1  /-. 

U  beiiiK  Uie  radiu;;,  \V('  describe  the  aiiipliludt-  speclnuii  as  a  liiiiction  ol 

h.irnionie  ui'dei'.  Ih'rtnrbatioiis  are  llien  calenlaUal  accordiiiK  lo  tiu’  deiined  I 

S|)e('lruiu  as  follows 

(1)  (Generate  a  random  (wliilt')  phase  spcvirnm  |  —  7r,;r] 

(2)  iiiodnlal.e  eoiuph'x  ])haii<' spi'clruin  with  pre<lelined  am|ilitude 

speelrmii  t,o  deline  the  homier  speclniiii  ► 

(b)  inverse  I'  l'T  lo  obta.in  space  domain  representa.tion 
(■))  scale  to  desired  maximum  ix-rhirlial ion. 

(."))  iiiodidal.e  perturbation  willi  a  prescrilied  function  of  <lepth. 

» 

(  'lioise.s  of  amplil  mle  (power)  speci.ra  and  t  he  ilepl  li  lilodnlal  ion  are  dr'scribeil 
in  ihe  main  body  of  tin'  (taper.  hh\am|)ie:i  of  (lerlni bations  I'alcnlated  in  this 
way  ai'e  shown  in  l''iy,uie  7.  (The  analof;y  in  ill)  s|)herieal  cooi'dinales  wo\dil 
be  lo  (trech'llne  tin'  .sitln-rical  harmonic  am|>lilnd('  s|)<‘<'trimi  and  distribute 

the  power  within  each  def;ree  randomly  amon(>  all  ils  orders),  I 

ho  make  the  models  inoia'  rr’alislic.  incoherence  in  d<'|)lli  is  inliodnced. 

1  ,a  I  ei  al  perl  nrbal  i(ins  arr'  weiyjited  wit  h  I  he  sel  ol  de|>l  li  depem  leiil  I  inn' I  ions 

-  1.  "  T). 

/„(•.')  '  sin(''-p),  I)  :=  l,;t . 2m  1  1 

<-os(’'^^-' ).  11  ---  2.  1 . 2m, 

and  2  I)  velocily  |iei  tm  bations  are  oblained  as 

I 

in 

.\v,{zj))=  )_2 (b) 

lO  I) 

where  a.re  the  lateral  (lertui bations,  and  ii'„  aie  weij'b's,  I’kxaniples  r)l 

2  1)  perturbations  obtaimnl  for  the  spectra  in  h  ikin' (i  ar<' shown  in  h'i/pu'e  1).  ^ 


» 


» 


•  • 


•  • 


m 


\V('  UHC  111  =  5  for  all  2-0  models.  For  .sju-clrai  type  I  the  weights  in  (’<|uatioii 
(5)  arc  ivi  =  =  1.0,  W3  =  ic.,  =  0.8.  n;.i  =  w;,  =  0.4,  and  for  spectral  ty[)e 

II  all  weights  are  unity. 
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Figure  1:  Stivggerecl  grid  used  in  the  FD  scheme.  At  the  top  of  tlie  model  tlu' 
angular  grid  si)aciiig  Ad  is  increased  by  a  factor  of  2.  The  radial  grid  spacing 
Ac  is  constant  throughout  the  grid.  The  F/arth's  surfacr;  and  tlic  (1MB  are 
modeled  with  a  free-surface  boundary  cotidiiion.  At  the  level  of  changing 
angular  grid  spacing,  Lj'2  levels  of  interpolated  OrO  componejits  ha\'e  to  bi' 
added  to  connect  the  Iwo  domains  (L  being  the  length  of  the  s])ac('  operator). 
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Fi}j;iir<‘  2:  'I'lic:  ytaggcuocl  FD  grid  close  to  the  axis  0  —  U.  TIu’  eqiiidious  ol 
motion  are  singular  at  6^  =  0.  At  llie  axis  of  syininclry  •«,.  -  1)  and  a,._.  -•  0 
fire  imposed  by  symmetry.  Note  that  a  directional  force  iiiimt  at  grid  iioini 
(H-)  next  to  the  axis  of  .symmetry  Icad.s  to  a  toroidal  ring  source.  'I'lie  gi  id 
is  exlended  Lj'l  -f  1  levels  beyond  (?  —  0.  where  L  is  the  length  ol  the  spac<' 
ojjcrator,  to  enforce  the  boundary  condition. 


•  • 


Time  (seconds) 


Figure  3:  Comparison  of  DSM  (bold)  and  I’D  seisinogratiis  I'oi'  llie  isotropic 
part  of  I’REM.  The  dominant  period  is  21)  seconds.  Some  .arrival  limes  of  S 
]diascs  are  marked  by  arrows.  Top:  Epiceiitral  distance  8U  degrees.  Ilotlom; 
140  degrees. 
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Figtiroi  4:  The  eifect  of  tlio  change  in  grid  spacing  al  tlie  lop  of  tlic  niodid  for 
uiffereut  frequencies  (2o,  35,  anrl  GO  seo  uds  tiominani,  jieriod).  I  lie  source 
is  at  t,he  surface.  The  model  is  the  isotropic  part  of  I’Hl'hVl.  Arrival  limes  of 
some  pha.ses  are  marked  by  arrows.  I’or  each  frequ(Mic\'  band  three  seismo¬ 
grams  at  an  eplccntral  distance  of  80“  an;  shown;  [.  If.eferenc.e  seismogram 
obtained  with  a  very  dense  grid  (5120x1024);  I!:  Sei.smogram  for  a  1024x250 
grid;  III:  same  as  II,  but  with  a  2048x32  refined  grid  at  the  top  of  the  model. 
The  accuracy  of  the  surface  waves  is  con.siderahle  improved  by  the  grid  re¬ 
finement. 


> 


« 


4 


*14 


Epicentral  distance  'deg) 


f);  'lop:  1 /.'il.cral  velocity  perl, iirbat ions  wliiclj  vvci'*'  added  co  IMiMM 
at,  t,lie  lop  ol  l.lu'  iiiaiiUc.  'I'lie  [)('rl.\irl)al'.i()iis  ar<>  colicieid  down  In  a  depl.li 
of  2 50  km.  lioUoiu:  'IVavol  lime  |H'iiml»al ions  for  S  (.Siar)  oblained  l)\'  Uie 
IluyKoiK'i  iiiclliod  (.solid  lino,  V\ ilte  ol  al.,  l!)l)5).  and  llioso  oblaiiu'd  uilli  the 
I'd)  molliud  (dots),  d  lic  sam|)lii0',  of  the  I'd)  seisniojp'iuiis  was  '2/d  second:;. 


b  velocity  perturbation 


2000km 


l''i(’,uic  'I  I'wo  iliiiii'iiiiiiiiial  motli'ls  Itir  iln'  Iwn 
perl  in  III!  I  .ire  liiuil, cd  t,u  llic  lop  oOl)  |;iii.  '1 
Sci'  IcM  loi'  ilclails  alnail  (lie  dcpi  li  ilcpi'nilrn' 


Epicentra!  distance  (deg) 


Epicentral  distance  (deg) 


I'i.mifc  14:  SS  p<-‘rl urljatioiis  (lc!(.criniiic(l  IVoiii  sci^moufiims  of 

wii'icd  doitiiiiiuit  period  (2"),  50.  and  100  sec.ond.s)  Ixd.wecMi  50  .uid  151)  d('Jr^^'c'^^ 
cpic  rnlral  distaix'e,  T.'hc  rocc’ivcr  ;spaciiix  is  1 /.'t  dcgn'c.s. 


Log  power  power 


15;  I’owoT  Rpccti'a  of  travd-liiiK'  residuals  dctcrmiiiril  in  dillereiil. 
frcciiiency  baiuls.  In  earli  case  S[)cc1ra  for  5  raudoiii  r^■aliza^ions  nsiii”  Larlli 
models  with  identical  characteristics  arc;  stacked.  The  model  lipecira  art'  as 
in  l'’igur('  8  (lower).  'I'he  maximum  model  pertiirhation  Is  iWi .  The  sjfeclral 
decay  determiiied  from  SS  residuals  is  ciearlv  fret|uenev  tlenentlenl . 
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Automatic  seismic  event  recognition  and  later  phase 
identification  for  broad-l3ancl  seismograms 


(i.l  INTRODUCTION 

The  ccjutiiiuous  siiTain  of  data  recorded  by  a  modem  lu;di  fidelity  .seismic  st.at.iim 
i;;  only  u.sei'id  wlioii  the  portions  of  the  record  associateil  witli  dillcreiit  eveiiti;  ciiii 
hi'  isolaU'd  for  furtlier  analy.sis.  Subsc'c]uc'nt  interpietation  depends  nn  lieiii!':  abl" 
to  rliarar'tei  i.'ie  tie'  event  'generating  the  patteiii  ofoliseived  an  ivals.  With  a  iimhi 
station  network  an  approximatt^  riassilicalion  of  an  event  in  terms  of  epieeiitial 
di.stance,  azimuth  and  depth  can  i.irovirh'  a  good  starting  point  for  leliiiemeiit  ol 
event,  location  (see  e.g  Kennett.  Iftfl-a).  Particuhuly  foi  a  sparse  (dolial  iietwoik.  ,i 
is  planned  for  flie  mniiiforiiig  of  a  comptehensi\i'  lest,  han  treat)',  tlm  ijnality  ol  tie 
initial  location  estimate  is  ci  itical  to  the  (|ualily  ol  event  ehai  :\etei  isat  ion  In  sie  li 
a  context  it  is  very  important  t.o  he  able  to  determine  whi’thei  the  anivals  delerl'  d 
at  dili'erent  stations  are  likely  to  have  arisen  from  a  .singh,'  evi  nl  or  from  two  or  nioi' 
geographically  distinct  eveirts  in  a  limited  tiiim  period. 

The  procedure  described  here  is  designed  to  recognise  the  ])atLern  of  arrivals  as 

sociated  with  an  event  from  a  single  three-componeut  broad-baird  recorti,  and  to 

inovidc  a  preliminary  estimate  of  the  eiricentral  distance,  depth  and  a'/Imulh  to 

that  event.  The  method  builds  on  the  work  in  the  previous  chapters  d  and  5  where 
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we  have  introdviced  an  automatic  system  for  seismic  pliase  detection  and  analvsi-. 
which  produces,  in  real-time,  a  stream  of  phase  attributes  characterising  each  cl  - 
tected  phase.  In  this  chapter  we  show  how  these  attributes  can  be  exploited  wiili 
the  aid  of  a  new  automatic  reasoning  method  (the  “assumption  tree”  method)  t" 
combine  the  information  from  many  phase  detections  to  form  event  segment. s  whi'  ii 
can  then  be  characterised  in  terms  of  the  properties  of  the  source. 

G.2  IdECOGNTTION  OF  A  SEISMIC  EVENT 
If  we  are  to  be  able  to  characterise  portions  of  a  seismic  data  stieam  in  terms  ol 
events  we  must  be  able  to; 

( 1)  recognise  related  seismic  iiliases  as  components  of  an  event,  and 

(2)  provide  i)ussible  interj)ret.ations  of  tlic  nature  of  the  phases  based  on  their  observed 
att  ribules. 

.^s  a  workhiy;  delinition  of  a  seismic  event  we  have  used  the  expectation  tlmt  there 
will  be  ;i  number  of  /’  idiases  followed  by  a  number  of  S  phases  during  a  certain  tine 
interval  (sa.y,  21)  mimiles),  'llie  ijrcK.eduies  suggested  in  <-hapler  h  ])ruvidi.'  a  mean, 
of  recognising  tlie  cliaracLer  of  teleseismic  phasr-s  in  term;,  of  /'  and  .S'  arriv.als  nsiii;, 
for  example,  t  lie  V  and  6'  detection  sclumie  ba.sial  on  the  ielat.i\’e  distrilml  ion  ol  tlr 
energy  in  the  wavef(.)rm  Iretween  the  vertical  and  liorizontal  comjnment s,  lloweV'  ;. 
suf:h  a  measure  of  waveform  cluiraeler  ha.s  to  be  suirpileiuented  I'V  oilier  altribiil' 
In  suggest  possible  identilications  for  the  seismic  pha.se.  in  chapter  5,  we  lia',' 
shown  liow  a  waveform  segment  associated  with  a  phase  detection  can  be  specili'": 
ill  terms  of  a  .set  of  jiaramelcrs  based  on  a  model  of  the  seismic  wavelet  .  Tliese  l)a:.i. 
paranu'teis  can  tlien  lie  used  to  (.-xtract  live  tiliribntes  lor  each  delecied  jihas'’,  i  ; 
for  the  Jtli  jrhase  we  would  have: 

(1)  b  -  the  arrival  time 

(2)  a,  -  the  amplitude 

(3)  -  the  local  frequency 
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(4)  0,  -  the  azimuth  in  the  horizontal  plane 

(5)  01  -  the  angle  of  incidence  to  the  vertical 

We  can  represent  the  stream  of  detected  phases  as  a  collection  ol  specilicatioii.'.  i  I 
phase  attributes, 

,  0,,  0,),  F,+i(t,.n, fl.-i.  0.41).  --  (Ih  ' ' 

from  which  an  event  is  to  he  constituted.  Thtis  if  there  are  A',,  pha-ses  detected  loi  a 
presumed  event  we  have  a  set  of  5A’,,  features  (5  attributes  per  phase)  as  the  input 
to  the  classification  process. 

The  pattern  classes  we  have  to  recognise  are  seismic  events,  and  .since  dili'erciii 
epicentral  distance  or  ue])th  can  produce  totally  different  sets  of  observed  iiliase.-, 
we  will  endeavour  to  cla.s.sify  the  events  by  range  and  depth,  In  order  to  proviiK  ,i 
reasonably  com[)rcl'.ensivc  coverage  of  possible  wtivefields  we  work  with  DU  jio.ssibli 
distances  (dividing  the  range  from  0’  to  180'’  into  2°  intervals)  and  a  selection  o! 
4  depths  (0  km,  100  km,  300  km  and  GOO  kin)  as  in  the  tes/nOf  tiascl-timc 
(Kennett  1991),  In  this  way,  with  .360  pattern  classes,  wo  are  able  to  provide  a 
acterisation  into  shallow,  inlennediate  :uid  deep  events  ami  to  provide  a  reas  iiinhl' 
sampling  of  the  deiteudetice  on  opiceniriil  clisttince. 

.  I'he  pi  tblcin  of  event  recognition  can  then  be  viewed  .as  a  classilical ion  jiiebl'  in 
in  a  5A0  diiiicnsiorial  vector-space,  with  3G0  possible  destination  i  lasses.  In  uul -i 
to  uncler.sl  and  the  complexity  of  this  probhaii  we  have  to  ieco;',mse  that  th"  piia  • 
.summari  s  extracted  from  the  oliserved  dtita  m;iv  be  disttale.l  due  the  pies'  ii-  i  .■ 
noise  or  incomplete  (e.g.  the  amplitude  in.ay  lie  below  the  detectinii  I  hi  r'siml' 1 ;  1  -i 
eacli  of  tiic  360  pattern  classes  we  ran  I'.enerate  complete  informal  i'Ci  foi  ea  h  "I 
the  phas"  attributes  for  a  bioad  range  of  possible  phases  using  the  nr./c'd  innd'  i 
However,  the  A'j,  pha.ses  will  frequently  lepiesejit  only  asub.set  ol  the  jjo.s.sible  iilnr  ' 
and  so  we  will  be  faced  with  a  classification  Ira.sed  on  iiK'.omjjlel.e  data. 

(Dne  obviou.s  approach  to  the  c\ent  classification  problem  would  be  to  set  up  3611 
model  patterns  and  compare  any  set  of  observed  phase  attributes  with  all  tins'- 
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model  iuit,l.crns  to  lind  tlie  best  iiirilrli.  'I'lu  phases  associated  with  a  [)arti(nl:ii 
Koiirec’  dc[)tli  and  distance  will  have  dideieiil  degrees  of  import mu'e;  thus,  somr 
])hases  must  be  observed  Ibi  this  combinalioii  to  In:  iclentihed  wliilsl  others  mi”, In 
or  nii{’;ht  md.  be  oliserved,  'I'he  attributes  lor  each  plitisc:  must  tilso  be  re;’, aided  a- 
Imvin;’,  dill'erenees  in  signilicmice,  since  e,;',.  the  relative  amplitudes  ol'  phase:-,  u  ill 
vary  de|)endiug  on  the  source  mechanism.  W’lieii  the  coinpai  ismi  belwi'i'ii  the  pha.'-e 
attributes  and  the  values  lor  a  model  source  is  made,  the  \vei;',ht.s  mid  iec|uii'i’d  level 
of  lit  for  dill'erent  attrilniles  need  to  be  flexible,  d'hi.s  means  llial  tlieic  will  be  o| 
the  Older  of  liA'j,  X  flGl)  features  to  be  individually  considered,  i.e,  tyiiically  several 
llumsand  cases,  In  addition  any  niodilicatioii  to  the  sysleiii  (e,;;.  to  inciease  the 
samplin'';  ol  tlw  model  s|)ace)  will  impose  a  recoieiiderat  ion  ol  an  aheaily  i:oniple,\ 
ab’,nril  Inn. 

Instead,  we  employ  a  new  strati’;', y  based  on  “assiimpt  ion  liee;;”  wlii-ii  allow;,  tin- 
pi  oicressi ve  iiu'lusion  of  seismolo;',ieal  expertise  lo  i.solat.e  the  mo:. I  piobable  conilii 
nation  ol  ejiieeiit.ial  distain  e  and  source  di  pih  for  the  available  pha.se  attribute:, 
'rhis  procedure  provid'":  a  more  )’,raee|ul  and  elileient  ploce.lille  lliali  exlnm..li\'' 
malehlli;',, 


ti.ii  'I'lii'i  usi'i  Ob'  y\ssuMr'j'i(>i\  'rmiMs 

Ib.'t.l  lindimioii  of  seisiaologie.al  ex|)erlise 

In  ordei  to  ledlice  t.he  computational  elloil  i.  -nmeil  U.  Imd  I  he  .ipplopliate  some, 
distance  mid  dept.h,  we  t.iy  to  build  In  as  mueii  ■  i  iiiologieal  evperli:,i'  a.;  po'ciible 
111  parliculai  we  lu'cd  to  take  account  ol  lie  way  in  wliidi  lie  p.illeins  of  :,ei  . 
niolo;',i-'al  phases  vary  with  epiceiitral  dii-lmi'i  boi  exani|)le,  n|i  In  SI'  tin:  maj"i 

observed  jihases  :ue  /’  mid  .S',  but  beyond  81  .S/\'.S  ]irecede;.  .S'  and  .appeals  a.s  lie 
liist  ])liase  'with  an  “.S-wave”  chai.ictei  In  oiliei  c  ,ises  M;',mlicain  inioimatioii  i- 

hi;’;hly  distance  siiecilic;  thus  lia.s  a  I'oi  ns  near  12(1’,  wliih.l  near  lli'd"  ,S7\’/' 

and  I’KS  are  m-mally  the  hugest  phases  in  the  eaily  jiart  of  the  .';eismo;';ram. 

Ill  order  to  identify  i)ha.ses,  we  need  iuforniatiou  about  the  epiceiitral  distance. 
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I'iiblc  G.l.  F-S  models  and  related  ranges 


P-S  iilciitity 

Time  iiilei  val  laii"'-  (niiiintes) 

DLstaiuitf  ( '  ) 

P  ■  S 

2.18  ■  10.4.3 

12  -  85 

P  -  SKS 

9  .37  -  10.0.3 

82  -  99 

P.iM  -  SKS 

9.39  -  10.09 

100-  12'.l 

PP  ■  SKS 

3.00  -  7.22 

82  -  129 

PP  -  PS 

9.20  in.M 

lO-i  -  12 .5 

PKP  -  SKS 

5. SI-  7.21 

li  t  -  143 

PKP  -  SKKS 

7.97  -  12-1.0 

120  -  LSI) 

PKP  -  SS 

20.00  -  27.GG 

130  -  ISO 

I'KIKP  ■  PKS 

3.35  ■■  3,00 

120  -  141 

Uii  tlir  oilier  hand,  lu  oliLain  the  estimation  ol'  epiccntral  disl aiice  i'lom  a  single 
us'oi  d  we  have  to  idenliry  at  least  Uvo  phases  -  we  can  I  hen  use  diflcreiil  ial 
lime  hetweeu  the  ideutilied  jihases  to  determine  Iho  distance.  In  order  to  make  llu 
process  tiaetalile  we  make  a  set  of  hy])otheses  aliout  the  idenlities  of  two  obseived 
|,ihnses.  'I'lms  we  .select  the  lir.st  plnise"  .and  the  first  'aV  phase'  in  a  ,'ieqtieiice 
a;.,  the  k''y  jdiasi'.s  d'ho  |)hase  chara-'ter  :'.ssi''iunenl  will  he  hier-rl  on  nmasures  such 
a:,  /’  and  iS'  as  disen.ssed  in  chaiiler  .a.  We  iIh'ii  try  t"  a.ssoeiati’  the  pail  ol  pluev. 
with  (he  most,  commonly  ohseived  eomhinations  of  /’  ntul  .'i'  phases  lot  dilTercn: 
epiceiit  I  al  distaiu  e  ranges  as  listed  in  'i'ahle  (i.l .  'i  lie  de|a'nd'  la'i'  ol  the  dilfei  cut  ial 
t  iiiieii  foi  these  /’-.S'  coiiihinations  on  ejiieent ral  distance  ami  depth  ai i'  illitst rat ed  in 
Id;', till'  (i.l.  We  I'an  .see  that  most  of  these  jiaiis  of  pliase.s  have  sulliiieiit  vaiiatimi 
ill  dillereiit i.al  t.ime.s  to  proviile  .some  eunsU.iint  on  ran;',''  :md  d'  plh  llowevei  ,  tie  i' 
will  oftciii  he  amhi,f',uit,_v  in  t.lie  approjiriaii'  comhina.t  ion  of  epieenlial  di.stauce  aii'l 
depth  a.ssociated  with  a  partiimlar  dillc'cnl ial  time  which  ran  only  he  re.iohad  h;. 
till'  mat  I  liiiie;  of  olhei  jihase.s  in  the  lecord. 

We  will  employ  this  .set  of  nine  /'’-.S' pairs  :is  the  lir.st  sta<>e  of  hypotlie.si.s  (  est.inn,  ac 
indicated  ill  I'd/pire  0.'2,  which  j eirre.sents  the  classification  hieitucliy,  .Since  there  i;. 
some  decree  of  overlap  between  the  operative  distance  intervals  for  a  iiiiuilrer  of  tlu' 
I’-S  ])airs  the  tree  .structure  in  Fieure  0.2  is  not  a  .strict  classification.  The  iimnher 
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of  leaves  on  the  tree  is  actually  more  than  360,  whilst  the  tolal  iiimihei  of  diflnci,' 
pattern  classes  that  they  represent  remains  3G0.  In  the  followin;';  sia  i  ion  we  devcl-.j' 
the  new  strategy  of  an  assumption  tree  based  on  the  stiucUnc  illusliati'd  in  I'igu!’ 
G.2. 

6,6.2  All  a.ssiuniitiou  tree 

IJefore  we  inl.roeUice  our  '  assumiilion  tree”  we  will  brielly  review  llie  ''decision  tii  i 
•sLiatej'.y  which  is  frequently  used  in  i^atlern  classilication.  Decision  trees  ust:  a  ■ 
quential  decision  making  strategy  to  cla.ssify  a  vector  of  features,  At  e;irh  liranchin;’, 
point  in  a  deiusion  tree  the.  branch  with  which  the  feature  vei  ioi  will  be  associated 
is  (leLenuincd  by  the  evaluation  of  a  lest.  A  coimnonly  used  test  is  the  hyperplatii 
le.sl  (  b,  Hreimnii  el  al.  IDU'I)  in  which  the  fcaturo  vector  >;  is  lesled  with  lesp.  ,i 
to  a  hyjierphuie  in  model  aimce.  The  l.irancli  i.s  then  delcnaineil  by  which  ,‘;ide  "I 
I, he  hy[)erplaue  x  lie.s,  leading,  to  a  binaiy  tree  structure.  Buildiii!',  a  decision  Ire.- 
retiuires  the  consti uctioii  of  appropriati;  hy))erplane;;  at  tin'  luam  hin;;  |H)iuts. 

In  our  ss'slein  a  multi  Inanch  tree  is  constructed;  this  tree  eunv.'-.  its  hrc.iiclic,'. 
a  set  of  tissuiuptiuns  ;uc  te.sled  aii  1  :ui  \m:  call  it  an  ''as  'm-pii-C!  I:',’'',  .Such  m 
assumption  tree  i.s  desij'.msl  to  (h'al  with  a  chc.silical  imi  p'  in  ivhich  thci'c 

a  immliei  of  factors  to  he  ideiililii'd  fsee  k'igurc  6.3). 

At  each  level  of  the  assumption  tree,  a  set  of  a.ssiiinpti  ii-  about  a  sin;’;lc'  laep  ; 
are  tested.  'I'lie  set  ol' stems  sjningin;;  liom  the  same  1  lilic,  p'lini  v.il'  1 1  ■ 

tlie  full  raiig.i:  of  possible  hypotheses  connected  v.ilh  thi:-  I.e  I'U  l.\'  iy  biaiidm. 
point  (nude)  in  the  tice  has  an  associ.'ii.ed  jjackaj'c  ul  r  !■  e  m  ■' j.  ,;i  ;\i  tie'  i,,  . 

h'vcl  the  inroMiiatioi)  consi.sl.s  of  just  the  oliscried  data  im'-'.  node  in  I:.' 

assumption  l.ree  will  inherit  .'dl  the  iuronimtioii  of  its  parei,-  n  ef  ',  and  iiielude  e.\ii. 
iiifoi  Illation,  d'lie  new  iuforuiation  includes  the  premise  whieh  lia.e  jie;t  been  nm  1 
to  produce  this  riotle,  and  the  inferences  based  on  the  jiarliculai  asi.niuiiLioii  iliawn 
from  the  siumuaiy  of  expert  luiowl^lge  held  in  asepaiate  knov.led;',e  haul;.  '1  li' 
information  jiackagc  at  a  node  is  examined  to  look  for  any  conl.i adiclory  rc.sulls.  II 
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any  cont.radiction  is  found,  the  bianrhing  associated  wiUi  the  cunenl  node  will  h- 
I ciiuiuat.ed.  Otherwise  the  braiudiing  inoeess  will  continue  from  the  noch'  and  alli'’.‘. 
the  evaluation  of  further  hyirotheses.  At  any  node,  the  new  information  introdiued 
from  the  knowledge  hank  may  limit  the  choices  for  other  currently  undetcrminid 
factors.  As  a  result  the  pattern  of  branching  at  a  node  will  depi’iid  on  the  curn'ii; 
information  package.  When  no  further  assuniiUions  are  to  he  tested  we  reach  a  leaf 
of  the  tree  rejneseiil  ing  a  solution  based  on  a  set  of  hypot Ihcsi's  alunil  the  fai'tou. 
dc.scribiii'’;  the  data. 

The  difference  between  an  assumption  tree  and  a  decision  tree  is  that  at  each 
ste|)  the  decision  t.i'ee  is  const.ruetcal  using  a  specific  ch-cision  ciiterioii  (knowledjte). 
wlu'i  ea.'i  ill  ,111  a.ssiuiii)l  ion  tree  the  cunstnictiou  depeiid.s  on  the  testing  of  hypotlii’se.-. 
In  the  decision  tiee  strategy  ex|)ert  knowledg''  is  irsed  to  determine  wliieh  branch  li 
choo.se,  wdiile  in  an  as.sunipl ioii  t.ree  such  exiiert  knowledge  is  aiiplled  to  l.eniiiiiale 
a  liraiuii  associated  with  a  ]);uiicular  inoiiiise.  A  diudsion  tree  is  pre  designed, 
that  adding;  a  new  decision  criterion  (additional  kiiowledgi')  reipiires  niodilying,  lie 
ili'cision  t.ree  and  hence  niodilying  the  pattern  classilical i' m  program  with  ollni 
pot  eiil  ial  side- ('I feci  s,  l'o|-  the  assumption  tree  iiroc,  dm  ■■  empl'gvi  ,1  ;■  i  lii ,  i  i  : :  ; . 
.elding,  an  addil  ioind  item  of  expert ise  rei'imii  .^  i  he  addit  ion  i ,1  an  item  m  i  In'  sepmsc  - 
kiinwh'dio'  baiilr;  this  pmeess  is  nnnii  .siiiipl''i  ami  sai'i  than  i  ec  oiisi  i  net  ing,  t!i' 
iniigram 

II. d  IMPLEMENTATION  OF  d'llL  ASSUME  I’lON  TiU’.k; 

I’KOCd'lDUE.E 

(i.d.l  I’liase  deti'ctioii  and  fi’alure  extrm  tioii 

In  (liaptci  ,h,  we  have  descrihed  a  niimtu'i  of  techiin|iies  lor  si'isiim  ]ihasr  lealni- 
extint  'inii  parlicidai  for  tlireocoiiipoiient  lecords.  Also  desi  lihed  is  ,i  new  teclini(|m 
based  on  ;i  coiuirlexlty  iiieasiiii-  which  can  detect  an  iiicicase  in  l're(|uemy  content 
for  plumes  whose  signaht.o-iioise  ratio  less  than  one. 

In  the  apidication  of  these  tcchiiiciues  to  the  problem  of  event  iilciitihcation  we  li\ 
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to  ensure  thrt  the  maximum  opportunity  is  made  to  extract  the  relevant  ])hases.  I'lii 
first  stage  uses  a  detection  system  trased  on  the  compicxity  measure  to  detect  tli- 
a])proximate  position  for  a  pliase  package.  Then  the  oiigina'  ■■  sinogram  is  divided 
into  a  number  of  segments  wliicli  contain  different  phase  packages  for  each  of  ivliich 
the  local  freciuency  can  be  estimated  as  described  in  chapter  b.  Earh  seismogram 
segment  is  then  filtered  n'o  a  high-frequency  trace  and  a  low-frec|Liency  trace  uith 
filter  parameter.s  ba.sei.  on  the  local  phase  frec|uency.  Finally,  a  set  of  detections  air 
made  on  boih  the  low-jrass  and  high-pass  lilteied  traces  using  adapti\e  STA/LTA 
(leteetoi's  which  coiiibine  energy  measures  on  one  or  more  components  of  ground 
moT-  'u  with  a  local  frequency  analysis.  To  simplify  the  input  for  event  ideni  nication. 
v.a-  choose  to  detect  only  the  first  phase  v.diere  there  arc  complicated  phase  packages. 

In  order  to  be  able  to  gain  icliable  information  -rbout  a  phase  wc  want  not  only 
to  de' ect  the  arrival  hut  also  to  characterise  tire  phase  in  terms  of  a  phase  vcctcu 
coin]'irising  a  summary  of  the  iiiformatioir  m  the  arrival.  By  filtering  segment  hy 
sf'gu.ent  we  avoid  problems  associated  with  a  vide  variation  between  the  frequeiirie.' 

(  ii  the  three-components  and  by  using  a  common  filtering  for  all  thrt.’e-comiJoiu-mt s 
■  !'i[)r(i\’e  the  rehai'/iiiiy  of  the  phase  \  ector  iiifoniiation 

<').4.2  Pi onrnresshrg  -  separating  redevant  phases 

!■'  chairter  5,  '..'e  have  shown  that  P  or  S  cliaraclei  can  be  as.-iigned  to  a  di'lecUM  pliasr 
b'.,!'.  .;()  on  the  relative  strength  of  the  cnetgy  on  the  vertical  conipoiienl  comirarcd 
V  ,th  the  total  ene’gy  on  all  three  comiroiients.  For  a  lelesei.sniic  ex’ciit,  a  dctecli  d 
phase  which  has  more  tli.ui  cue  third  of  it.s  total  energy  on  the  ccitical  comptmeiii 
c:an  be  recognised  a.s  a  '  F  phase" ,  while  otherwise  it  will  be  assigned  a.s  an  “S  phase'  . 

'mce  an  event  on  n  seismic  record  can  be  rec'  prised  as  a  nLiiiibei  of  P  phast',< 
followed  by  at  least  one  S  phase  ar.-'  possibly  a  number  of  otlier  P  or  S  phases 
within  a  30  minute  winuuw,  a  seismic  event  can  be  formally  represented  as: 

(an  eve'-’t)  — ■,  (first  P  phase)(P  phase)*(firsl  S  phase)((S  phase)|(P  phase))*, 
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where  tlie  notation  (A'”)'  means  that  X  may  appear  many  times  (inciuditig  zetd 
while  A''|V'  means  ‘'either  X  or  Y 

This  event  definition  is  used  to  assemble  a  set  of  detected  phases  (and  their  ass'  - 
ciated  parameters)  into  a  list  which  will  be  used  as  the  basis  for  cu'ont  classificat  r  n 
and  further  phase  identification. 


G.4.3  Seismic  event  interpretation 

The  system  of  event  processing  is  based  on  an  assumption  tree  with  a  similar  struc¬ 
ture  to  that  illustrated  in  Figure  6.2.  The  assumption  tree  grows  at  run  time  a.s  we 
test  a  set  of  hypotheses  about  the  observed  seismic  phases,  branches  that  lead  tn 
cocitrc.dictory  results  svil!  be  terminated  and  not  giow  to  leaves  reprc.senting  a  vial  l 
event  classification. 

6.4. 3.1  Stagr  1  -  c.'to'cc  of  P-S  pn'r 

The  ini* in'  information  package  comprises  the  s'>t  of  informal  ion  on  the  detect''.! 
lihascs  at  the  root  node.  Nine  stems  are  then  grown  to  cliihl-nodes  repre.-jentii,., 
different  assumptions  about  the  identity  of  the  combination  of  lirsi  P  and  first 
pliast’s  (as  indicated  in  Figure  6.2).  The  couseciueiices  of  these  assumjitions  are  tin  t; 
te.sted  against  the  expert  knowledge  held  in  the  separate  ktiowlcdge  base. 

If  th.;  assumption  about  the  P-S  phase  pair  is  true,  the  event  must  lie  within  tie 
disttiticp  range  for  this  combination  (see  Table  6.1),  and  the  (liU'ereiit  i;il  time  bciwe'':: 
tlie  first  P  and  tire  hist  S  must  he  inside  the  pos.sible  P-P  titne  intcr\-al.  The  possil.'l' 
distance  range  is  thereby  reduced  from  0''"  •-  180'  to  the  dtstance  range  ajjpropriat' 
to  the  particular  choice  of  phases.  The  observed  P  to  S  diffei cntitil  time  can  tlu  ii 
be  compared  with  the  expected  range  for  this  phase  combination.  If  the  observed 
differential  time  doe.s  not  lie  inside  the  expected  range,  we  liave  a  contradiction 
of  the  assumption  about,  the  pha.se  pair,  mid  the  branch  of  tlic  a.ssumption  tree  is 
terminated.  The  phase  pairs  in  Table  6.1  have  some  overlap  in  differential  times  so 
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Table  6.2.  P-S  pairs  and  iclated  feature  properties  for  earthquakes 

• 

% 

P‘S  idcntily  Expected  featiircii 

(») 

1’  ■  S  The  fiequency  of  S  should  be  lo'vci'  lli.-in  tlii'  fri'qiiLiu  y  of  P: 

1 

P  should  not  be  too  steep  (a.s  PKP]: 

there  should  not  be  any  low  frequency  S  pha.'i'  befme  the  S. 

JtC! 

P  -  SKS  SKS  should  not  be  a  high  frequency  phase 

Tiiff  -  SKS  P  lia-s  a  smaller  amplitude  compared  with  SKS 

» 

PP  -  SKS  PP'.  long  period  and  shallow  incidence 

PP  -  PS  PP\  long  period  and  .shallow  incidence 

yVv'/’  -  SKS  PKP:  short  period  and  steep  incidence 

PKP  -  SKKS  PKP:  short  jieriod  and  steep  incidence 

PKI  -  55  PKP:  short  period  and  steep  incidence;  SS:  long  period 

» 

‘ .  :• 

I'KIKP  -  PKS  PKIKP:  shoVt  period  and  stec.  incidence 

T  lit?  "srei'o"  e?-  ”cn.,iPy.“  p.-  fpp  c.r  lli(. 

smaller  t  '  ■■  i;.'  id  -r.'-.e  an;.  '".  tle-'.-P  •: .  •  'PP',.:.  PKKP:.  -  ■  ch"’;  f.c.-  \  i'ry  distimt 

events  for  whicli  the  incideiiee  angle  is  very  stei'p.  noniiallv  le.ss  than  10'".  Thus,  an 
olj.served  pha.  e  s’-.e.r.'..!  l.axe  an  incidence  angle  i-i'ialler  than  10'  to  he  matched  with 

I'KP  01-  PKKP.  Tlie  value  of  IG"  has  been  obtained  by  trial  and  error,  and  is  subject 

to  furtlier  adjustment. 

» 

• 

tluit  at  ino.st  5  possible  assnrniitions  about  the  /''-5' pliri.se  chnrarU'f  ran  survive  thi- 

test . 

At  this  stage  we  can  apply  additional  scismological  iufornialion  to  test  the  din'cicnl 

premises  about  the  P-S  pair.  The  jiroperlies  which  have  employed  are  tabulated  in 

Ttible  0,2,  ■ 

I  '.'V 

G.4.3.2  Stage  t’  -  the  depth  of  the  event 

t 

.A.t  this  stage,  every  surviving  first  generation  node  grows  four  new  .stei  is  to  child 

.y. 

— 

nodes  (see  Figure  6.2),  where  each  child-node  is  associated  with  a  differeril  assuni];- 

tion  about  source  depth.  Based  on  the  assumptions  about  the  identity  of  the  P  and 

S  phases  and  depth  and  the  observed  P-S  differential  time  we  are  able  to  isolal'- 

a  single  possible  epicentral  distance.  Nine  differential-time  tables  (one  for  each  of 

the  P-S  pairs)  connect  the  P  to  S  time  interval  with  epicentral  distance.  A  search 

» 

h 

— 

• 

•  •••••• 

W 

m 

0 

'  i 

1 
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in  the  appropriate  table  will  return  the  distance  whicli  corresponds  to  the  small'  -■ 
difference  between  the  theoretical  and  the  observed  difl'crential  tiine. 

Since  we  find  the  distance  by  table  search  we  do  not  implement  the  third  gen-:- 
ation  assumptions  fur  distance  (indicated  in  Figure  6.2),  and  as  there  arc  no  mo:  ■ 
assumptions  to  test,  we  have  reached  a  leaf  of  the  assum]rtion  tree. 

With  the  estimates  of  the  depth  and  distance,  we  can  compare  the  patterns  of 
observed  phases  with  those  expected  for  the  nine  different  assumptions  on  the  P-S 
]rairs.  As  described  ab'ove,  at  the  most  b  models  for  the  P-5  ])air  can  survive  from 
the  first  generation,  Once  four  possible  depths  are  included  there  would  be  at  mo-i 
20  live  nodes  in  the  second  generation,  jince  20  is  not  a  large  number,  the  arialysi.-> 
procedure  accepts  all  the  second  generation  nodes  as  possible  interpretations  aiid 
th.cn  we  tiy  to  a.-ses.'-  the  i'kplihood  of  each  in^erpretatm:i  !'v'.h';'r  than  look  1 
contradictions, 

6. .{.3. 3  J (kntification  of  other  frequently  observed  phases 

So  far  we  have  used  seismological  information  on  the  likely  character  of  the  tm. 
main  phases  (the  first  P  ]iha,sc  and  the  first  S  phase)  which  are  required  to  form  a:; 
event  pattern.  For  a  choice  of  source  distance  and  depth,  we  can  make  use  of  a  s- : 
of  ‘'oilier  frequently  observed  phases"  as  summarised  in  Tabic  G.3.  the  arrival  tini-  - 
and  other  properties  can  then  be  calculated. 

Then  by  comparison  between  the  sequence  of  observed  jihases  (expression  tl.l'. 
and  the  phase  predictions  for  the  a.ssuined  distance  and  depth  we  ma>-  be  abh 
identify  phases,  and  further  the  presence  of  expected  pha.ses  can  jirovide  support  I  : 
the  assumed  depth  and  distance. 

The  specific  choice  of  .P-5'  pair  and  the  other  expeclod  phase.s  constitute  the  i:,- 
forination  set  to  be  compared  with  the  observed  set  of  ]ihases.  We  therefore  need 
to  match  the  observations  against  the  expectations  for  the  proposed  distance  and 
depth.  For  each  expected  phase  we  have  an  arrival  time  determined  by  the  pha.se 
identity;  and  each  detected  phase  has  a  set  of  observed  phase  attributes  which  in¬ 
cludes  the  arrival  time.  The  first  P  and  first  S  phase  in  the  observed  data  ha\e 


«  •  »  •  • 


O.Jf  Implementation  of  the  assumption  tree  procedure 


HI 


Table  6.3.  Expected  phases  for  choice  of  PS  pairs,  for  a  certain 

distance  and  depth  ranges 

P-S  identity  Other  oxpectod  yili.i.sp.s 


P  -  S 

P  ■  SKS 


Paiir  ■  ShS 
PP  -  SKS 


PP  -  PS 
PKP  -  SKS 


PKP - SKKS 


Di  ef)  event:  dejitli  plifvses  s/’, 

25"  -  05":  PeP.  SvP.  P(:5:12"  -  39":  ScSi  31'  ■  79b  PP. 
For  nil:  PP: 

O.j"  .  99b  PKKP. 

For  oil:  PP: 

115°  -  125°:  PKKP.  PS. 

115°  -  125°:  PKKP.  PS:  95'  ■  90°:  PKKP: 

Fur  nil:  SS. 

For  oil:  PKKP.  SKS,  SKKS. 

Dee])  event:  dejjlh  piloses  pPKP,  sPKP\ 

For  all:  PP. 

135'  -  143b  SKKS:  115°  -125b  PKKP,  PS: 

125°  -  113°:  PKS:  114°  -  131°:  SS. 

D  .".Viit:  <le;  ‘b  ;  '  •  --  pFKP.  sPKP\ 

Fur  all:  PP: 


125°  -  139b  PKS:  137°  •  180b  SS:  159"  -  181)°;  PKP,c... 
PKP  ■  SS  Dec])  event:  dejith  pliosos  pPKP,  sPKP: 

1  11-  -  1.SO-:/'/’.  , .  isilb  PKP„i,. 

PKIKP  -  PKS  Deep  event;  deptli  ])li.r.ses  pPKIKP.  .iPKIKI  : 
lor  all:  PP.  PKKP.  SKKP: 

1.37=  ■  143°:  SS. 


alread}'  been  associated  with  two  jihases;  \vc  can  Uierefoii'  establish  a  .siiiijjle 
mapping  belween  tlte  observed  dalr.  and  the  ex)5ccted  arrivals.  W’e  search  aiii'-:; 
the  observed  arrivals  lor  candidates  for  an  exjjected  phase  whose  arrival  time  la 
with.!]'  a  j.iredefiried  tolerance  of  the  expected  time.  The  ob.served  a.ttribntes  of  tli' 
phase  are  then  compared  with  a  summary  of  phase  properti  in  tlie  knowledge  ba.^" 
(see  Table  6.4).  If  there  is  no  match  then  that  hypothesis  for  the  phase  identification 
is  rejected,  but  if  more  than  one  phase  match  is  possible,  the  one  with  the  .smallesl 
time  difference  between  observed  and  expected  times  is  selected. 
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Tabic  6,4.  Expected  properties  of  other  frequently  observed  phases 


Plia.'io  identity  bxiicttcd  features 

pP.  sP  pP.  sP  arc  similar  lu  P  both  in  freaiueucy  and  incidence. 

PcP.  SeP  PcP.  SeP  are  hi”h  frequency  P  plueses  wliieli  is  .mee|)ei'  than  P. 

PP  PP  i.s  shallower  than  P  and  lower  frequency. 

PKKP  PKh'P  is  a  steep  and  hi,!;h  frequency  P  phase. 

PKP  PKP  is  a  steep  P  phase. 

Sr.S  SiS.  high  frequency  steep  S  i)hase,  strong  on  langeiitial  cnniponeni. 

PrS  PcS:  high  frequency  steep  S  pha.se. 

PS  PS  is  shallow  and  low  frequency  S  ph;\sc. 

SS  SS  is  a  low  fre((uency  S  phase. 

SKS.  SKKS  SKS.  SKKS  can  not  he  high  frequency. 


6.4-S./1  Estimation  of  likelihood  of  hypotheses 

In  order  to  estimate  the  likelihood  of  any  particular  interpretation,  we  considi’i 
both  the  ob.scrycd  data  and  the  expectation  for  a  partierdar  coinbiiiatiun  of 
pair,  distance  r.nd  deptlr.  When  we  compare  the  two  sets  of  ]/.ia.,e.s  tve  ha\c 
thicc  different  cases;  matched  pha,ses,  unmatched  fthases  in  the  observed  data,  and 
unmatched  pha.se.s  in  the  expeci.ed  data. 

For  every  matched  phase,  a  time  error  c,-  is  used  to  it’iuc.senl  the  inismalch  helweei. 
the  corresponding  observed  and  expected  phastrs.  We  set 

e,  =  ((i.P) 

‘  .<  I'V 

where  i„  is  the  pliase  time  in  the  observe' I  data;  is  the  idiase  time  for  the  e.q)ec  li'.  1 
phase;  t.,  is  the  observed  time  of  the  key  S  phase;  and  f,.  is  the  observed  time  th 
the  key  P  phase.  For  every  unmatched  phase  (both  in  llic  observed  data  and  lie 
exirecled  phase  list),  we  set  e,  =  1. 

In  order  to  provide  an  overall  measure  of  the  quality  of  the  match  between  the 
observed  phrrscs  and  those  exi)ected  for  the  particular  distance  and  depth,  we  weigh 
each  observed  phase  with  its  amplitude  (i.e  we  place  more  weight  on  large  ar¬ 
rivals).  For  those  expected  phases  which  have  not  been  siicce.ssfully  matched  witli 
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any  observed  phase,  the  smallest  amplitude  of  the  matched  cnmiioiu'nt.s  is  aiipre  d 
as  the  weight.  We  denote  the  weight  for  each  phase  a.s  ii.',,  and  tlu-n  define  the  mislc 
niea,surc  <5  as 


6  = 


(ti.:;. 


where  the  sum  is  taken  over  n,  v.’hich  is  tire  miinber  of  phases  in  all  the  tlnee  ra.si's. 


6. .{,3. 5  Details  of  the  implementatioti 

We  take  a  pragmatic  approach  to  the  seleclion  of  phases  for  the  key  y-’-.S’ pair  .siiic  i 
we  have  used  the  most  commonly  observed  jrhascs;  rather  than  insist  on  tire  “hr.st" 
a-rival  of  particular  type  we  look  for  a  jnominent  arrival,  e.specially  for  S.  Aflei 
the  detection  of  a  “first"  /’phase,  we  choose  the  first  S  phase  encountered  in  llu 
processing  sclw.ue  n-s  the  candidate  for  the  key  S  phase.  Ilowmci,  if  a  .second 
phase  orcur.s  wiifiin  -h  minutes  witlinnt  an.v  mlcrvening  phase  and  is  larger,  it  will 
replace  the  fir.sL  candidate 

In  general,  as  frointed  out  in  chapter  4.  a  low-passeil  soisinograni  gives  a  bd 
ter  definition  of  phase  attributes  than  the  corresjronding  liigh-pn.ssed  seismngraiii 
'I'iieiefoi  c,  when  wc  select  the  tv.'  key  phases,  we  choose  low- frcipiency  detecl  ion - 
We  use  the  azimuth  estimati'  for  the  key  /’|ihase  as  tlie  aximiil  h  of  the  whole  e\eni 
I'oi'  very  di.stant  events  from  110'  u]iward.'  the  liisi  /’airi\'ai  i‘.  usuallv  scmiewhm 
weak  and  the  resulting  azdiniith  estirnale  ■  somewhat  douhtliil. 


(i.5  EXAMPLES  OF  EVENT  RECOGNITION 
'J’lie  event-recognition  .system  des  abed  in  the  laevi.  us  section  ha.s  licmi  apiilied  to  a 
wide  range  of  observed  seismic  records  with  considerable  suece.'.s,  Wc  have  select ''d 
a  set  of  four  examples  which  repre.seni  events  from  (liffereiit  distami's  and  depths. 
We  use  seismic  records  from  portable  broad-hand  inslnimeiits  dc|doyed  hi  tic 
Northern  Territory  of  Australia.  The  data  weic  recorded  on  H.eflek  72A-07  di.'v 
recorders  (24-bit  resolution)  with  Giiralp  C.MG-3ESP  seismometers  (flat  to  ground 
velocity  from  0.03  -  30.0  Hz)  with  Omega  liming.  The  sampling  rate  was  25  samider. 


6.5  Examples  of  event  recognilwn 


Table  6,5.  Events  used  for  illustration 


Event 

YortI' 

Day 

Titnc 

Station 

UislaiieL'l''| 

Do|ith  [kin] 

A/iiniitli[°) 

,M0 

A 

199t 

281 

21:4f:55.K 

SC()3 

15.4 

17 

175 

(7  1 

U 

1991 

22(1 

21;15:3().0 

scot 

52.0 

122 

133 

(1  11 

t: 

1994 

244 

irn2().'18.8 

SOUS 

113.8 

1(1 

201) 

(1.1, 

IP 

1994 

too 

00:40:55.9 

SC()9 

137.1 

G31 

213 

7.1) 

pel  .seeond.  Such  data  is  of  good  ciunlily  hut  not  of  ohservaiory  standard  and  se 
re])rosentK  a  useful  test  for  both  event  recognition  and  phase  identification. 

Idle  four  events  we  u.so  for  illustration  arc  listed  in  Table  (>.5.  ^Vc■  note  that  llie 
'^'•ent  location;:  for  these  events,  as  given  in  the  table  had  been  detonnined  usinr. 
inforuiation  from  many  stations  with  a  broad  azimuthal  coverage  wherea.s  oui  es 
liniales  are  bai-ied  on  just  a  single  three  componout  record.  The  time  listed  in  tlm 
table  IS  the  beginning  time,  of  the  record  from  which  the  automatic  sy.stcm  is  applied. 

In  every  ca.se,  the  proces.s  of  phase-detection  and  feature-extraction  is  repie.sentecl 
in  a  secpience  of  display  panels; 

(al  'I'iie  oi'igiur.l  hroad-’^’”' '  ■  'ord  is  ilisplayed  together  with  tlm  mii\al  tinv,,-  I  e 
die  major  jihasas  imedictod  from  the  iarpOl  travel  time  t.ahles  as  a  n'ferenc"  !■': 
the  idiase  interpietation  process. 

(h)  The  STA/LTA  detectors  generated  for  the  three  low  jiass  liltens’,  component .s 
,N,  b),  I’hc  arrows  indicate  where  phase  detections  !i.iv<'  been  made,  I'm  lieii 
l)hase  detection  the  phase  features  are  extracted  a.s  de.sciihed  in  chajitei  5.  1  Im 
relevant  waveform  segments  and  the  extracted  pha.se  attiihule.s  ate  disphiycd  in 
(-1) 

(c.)  The  bTA/bTA.  detectors  on  the  three  high-jta.ss-liltered  traces  are  displayed,  aloii;;. 
Nvitli  till:  arrows  which  indicad.c  phase  delerlions-  The  numbering  of  the  plnisi- 
detections  in  panel  (b)  and  panel  (c)  is  linked  to  panel  (d),  so  that  the  sequence 
is  in  order  of  arrival  time. 


(d)  The  last  panel  shows  the  waveform  segments  corres)  onding  to  each  detection, 
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together  with  the  ainijlitude  and  and  other  attributes  of  the  phase  (see  chapt'  i 
5).  In  the  waveform  segment,  the  time  interval  between  tiie  two  labeled  tine  - 
represents  the  assigned  period  of  the  detected  phase. 

The  figures  for  the  phase  detection  processes  are  accomiranied  by  a  summary  tab! 
of  the  stages  in  the  event  recognition  and  phase  identification  [uoceduu’.  This  tabl  ■ 
i.s  organised  in  terms  of  the  four  major  steps  in  the  analysis 

1:  identification  of  the  key  /'and  5  ])hases 

Step  '2:  comparison  of  the  differential  time  between  the  key  ]rhases  with  the  expected 
times  for  P-S  jjairs. 

Sir}  .‘I:  checking  of  cx]rected  /'-S pair  jrhase  attrilrutcs 

.^trp  inatcliing  of  observed  phases  against  e.K]3ecled  ]ratteriis  for  ijitiei'ont  models  (d 
event  character. 

The  liiia'  .step  is  selection  of  the  combination  of  distance  and  depth  which  leads  t" 
the  least  misfit  between  the  observed  and  expected  phases. 

(i.ri.l  Event  A  Ilulinahc"''  (  j 

'I'lie  event  rdiown  in  h'igure  (),4  is  a  region:’.!  {’wnt  at  lli”  distance  of  bV  .  In  ])am'i 
(a),  the  original  sci.sn.iic:  u-iouvl  i.s  dominated  b',  /‘ri  and  .S';i  ivickages,  haler  ]dias'  ■■ 
aie  mixed  with  surface  waves.  Since  the  /’  ariival  show.s  little  energy  on  llie  /. 
componc-mt,  the  seisniograni  is  almost  natnraily  polarised,  i.e,.  the  /•.' com])onent  i-- 
the  tangential  component  which  imeals  57/ w.ive  arrivals. 

Ill  panel  (b)  and  (c),  the  i‘r.  and  Sn  arrivals  are  detected  on  both  high-fretiueii' 
and  low- frequency  lets.  Since  ScS  is  strong  on  the  tangential  eonqroiumt,  the  hlgli 
freqiumcy  set  // componerit.  clearly  reveals  tlie  .Vc.S' arrival  (see  (])lia.se-,S]  in  panel  cl. 
'L’iie  detections  [phase-lj  and  [ph,'.se-2)  occur  at  almost  the  same  time  on  the  taco 
ddl'eiently  filtered  sets,  so  also  do  [phase-3]  and  (phase-4]  The  automatic  system  uses 
tne  aTributes  for  the  phase  on  the  low-frcqucncy  set  from  each  pair,  we  therefom 
cihoosG  [phase-1]  as  the  key  P  phase  and  [pliase-3]  as  the  key  5'  phase  for  the  difh : 


iDa.12n'.21 -IB.SS.S-IGr.COS  Depth  17  AppiO!t_distenc: 
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Fig.  6.4.  PlKise-dctection  and  featurc-cxlraction  for  event  A  -  Haliiiahera  (Indonoaia), 
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Talile  G.6.  Processing  sequence  for  Pvent  A 

Step  1:  Choice  of  Key  PIuiks 
Ki-y  P.  phase  1,  Key  S:  ijhase  3 

Stcji  li  'J'esL  (in.  Piffcrenliid  Tint'. 

ICT.S  s  -  sui  viviiiy  inoUcls:  I’-S 


Sleji  S:  P 

Suix  ivill" 

Step 

S  Pci 

moilcls: 

Pniliiic 

PS 

Mu  tel. 

K ft: III  Patlcni  Match 

P-S  p;ui' 

Depth 

Dist. 

A'/iin 

Identified  Phase.s 

Misfit 

r  s 

II 

1-1 

,1S0 

phase  1:  P,  phase  .3:  S',  phase  0:  SiS 

O.Oil.'l.'i 

11)0  . 

13 

180 

phase  1,  /',  jiliase  3:  A' 

0.333.1 

300 

10 

ISO 

pha.si.'  1:  P,  jihase  3:  S 

0. 3,1.1, 1 

GOO 

IG 

180 

l>hase  1:  P,  phase  3:  S 

0  3333 

niUal  lime  aiialy.sis.  'lahlc  G.G  suniinariscs  the  progress  of  the  event  inifiprct.  -: 
scheine.  After  the  lest  (jf  Ihe  dill'ereLitial  limes,  only  the  A’-5pair  sun  ivrs  'l  i;;  / 
lutulel  also  suivived  i.he.  lest  on  tlin  expected  /’-5 1’eaUii'es.  '1  hf'ielon-,  th'  i'  ,!;■ 

I'luu  nodes  in  the  second  generation  corresiioiuliii'';  to  difl’emni  inleiptetatioii  t 
event  depth, 

h'or  (O'ciy  elioice  of  seismic  source  inleipretat  ion  h)  .ste|)  '\  of  lahle  (i.''.,  tlm  ilnt  ;ir  • 
laiige  and  7’-.S'pair  load  ns  to  expect  thcarrival  of  tlie  jiiia.'v  Ac.''.  Im  tlm  liisl  ■ 
with  a  surface  source  (dejith  0  kin),  the  exiiected  time  I'oi  .S'r.V  conospunds  to  In 
oliservcd  (];liase'5]  and  the,  e.xiiccled  pha.sc  features  for  .S'cAaie  matched  with  tin 
attribute.s  of  [pliase-f)].  Therefore,  in  addition  to  the  two  key  pliase.s  a  luither  jilm' ' 
is  identified.  I'oi  each  of  the  other  tlirec  cases,  the  exireclcd  time  range  for  ScS  doe,, 
not  i:ovcr  any  ofiscrvcd  iihases,  as  a  resuli  the  expected  ph.asc  ScS  r-  not  found  and 
the  misht  measures  are  larger  than  that  for  a  .surface  source. 

When  we  compare  the  best  interpretation  (Distance  14^  Depth  0  km.  Azimuth 
180°)  with  the  event  specification  given  in  Table  6.5,  we  can  sec  that  the  estimations 


C.5  Examples  of  event  recognition 

of  dist;inco,  dnpUi  aiul  ;v/inmtli  are  very  close  to  tlic  Inic  values.  Al.  the  same  tiiu'  . 
the  aulomalic  system  Iku.  cotrectly  identified  three  phases:  [p!ia.se-l]  a.s  f,  [phase  - 
a.s  S.  and  [phasc-h]  a.s  ScS 


.IT.'  11  M'.  a.iiMi  a  Hiu  in. 


a  it<’d  liv  ch-ai  /’,  hi'  rnul  core  leflcc 
•  '  >;i  i  .‘  .nn'wd  .  are  detected  on 

■li'-:  i-i''  h-'iis  arc  [phase-d] 
:•  •  '  1  ,i  I'.-'  .ilso  indicate 

•  -  .  ■  -  ■  ■  a:''  la'i  fla^yr’'  I 

.:  ::.:'--hi>l  1  h.i 

i  htth  t"  ’ 


:  t : ,  •  • .  , .  •  1 M  ■'  n  i  ' 

'  tin  ti'iUni  i  I  i-r.t  <i:i  1  s  I  1  li  'i '  le ■■  h-  ■  in  t  le-  i .  i:i  1  eni,  i  .n ii >ii  in  1 1  ;• 

:r  auiiptiiiii  tree  eon  espondiii"  to  ditieieiii  i  ■iiiiliinatioii-.  of  di-t.ine,'  and  (li'])t!i  .i 
well  a.s  phase  choice. 

Idle  snudlesl  nicasuic  of  misfit  occurs  for  the  second  choice  in  step  A  of '.[’able  G.’^' 
rvliich  con  esponds  to  an  interpretatiou  of  the  key  phases  as  P  and  S.  The.  csthnalcs 
of  the  event  paiameters  arc  epicentral  distance  52",  deptli  100  km,  and  azinmili  510 


O 
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I 

I 

Tabic  0.7.  Processing  scguencc  for  Event  B 

Ste]>  I:  Choice  cf  Key  Plidsi-s 
Key  /':  pluvsr  2.  Key  S:  G 


■$' 


■K} 


Stc]>  2:  Test  on  Diffcrcnliiil  Tiro 

431.3  s  -  sui  viviii^  nuicJcls.  P-S  PP-SES.  PKl'-SKS 


5  p  ,1.  P-S  Pair  Prntuir 
Surviving  models:  P-S. 

Sop  / 

Mulch 

P-SKS 

EiLiit  Pullcrn  Match 

I'-S  |i:iil  Di>pth 

D>si 

Axini 

lo-.-ntifiocl  Phaso.s 

Misti: 

/’  S'  11 

5(1 

31C 

phase  2:  P,  phase  3:  PeP,  pha.sp  6:  S 

0.00100 

ino 

0'2 

3tG 

pha.se  2:  /',  phase  3:  PeP,  phase  4:  PeC, 

phase  C:  S 

0.000  ;■ 

aun 

51 

31C 

phase  2:  .P,  phase  3:  pP,  phase  4:  PcS. 

pha.se  0:  S 

o.oon.i 

GOO 

5S 

310 

ph.ve  2:  P,  phase  3:  PcP.  phase  4:  PcS. 

pha.se  G:  5 

0.00!  O' 

IT  5/\.S  0 

81 

31G 

plu'-SQ  2;  PP,  pha.se  G:  SKS 

0.33.TV 

mi.) 

82 

31G 

pha-se  2:  I'P.  [ihaso  C:  SKS 

0.333:; 

3ui  1 

82 

310 

p’l;-.--  2;  PI',  pha.se  G:  SKS 

0.33,1.3 

GOO 

S'2 

310 

p'.e.~e2:  P/t  |)liase  G:  SKS 

0.333,')' 

(which  al-s'..'  iui[)hc3  13(T),  which  nuiichs  the  iiifoi mation  i)roviclccl  in  'labile  G.-'i  \ci} 
well. 


(3. lAcnt  C  off  coasi  of  Northern  Clalifornia 

Event  C  is  at  a  much  greatci  distance-  and  so  is  in  tlic  core  sliadov,'  for  }'■  'I  lr 

expected  fir.sl  arrival,  difl'racted  F.  is  so  weak  that  it  can  not  be  seen  on  the  seisnr' 

gram  (Figure  6.Ca).  Core  phases  (e.g.,  FA'F  and  SKS)  and  surface  multiples  (c.g..  ® 

PP  and  PS)  are  the  main  features  on  seismic  records  at  this  distance  range. 

The  set  of  detections  displayed  in  panel:  6(b)  and  6(c)  pull  out  all  the  major 
arrivals  on  the  seismogram.  We  notice  that  the  detection  [phase-?]  and  another  high  > 


O 


ir))4  ir.20.’ '  fi-IC  sc:0  Ocplh  10  App.fOx  rJ-stancc  ns  so 
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G.i  Exarnjtlcs^  of  i:vent  rei'Oi)nihon 


I'j:. 

rreciuciir^'  S  ai  '  ivat  liavc  '.o  obvious  cur.oi'ii  ic.i  v.  jtli  llu'  phases  oxpeci'  ' 

for  this  even  and  ii’.ay  v.ell  come  from  some  differoiU  souiti.. 

The  deteclion,  on  tlic  low  frcciuciicy  set  [phasc-l]  is  selected  as  (lie  key  7’ plia' 
Note  Liii’J  s', lice  [pliase  ll  is  bihacr  '-han  ipliasc-5]  and  iheie  T  no  detection  ’in  I,. 
UvoeU;  [ph;‘se--'ij  is  selected  as  the  hey  S  pliase  in  place  of  [pliase-h;  Tlte  progrt  .-- 
of  die  aul.jniated  event  recognition  s’,  rleiii  is  presented  in  Table  (IS,  In  the  lir.-:-’ 
geiicratiiii,  in  the  assiiiiijition  tree,  five  set  of  F-S  pairs  sinvived  the  test  on  tie 
d'Tereiidal  time  between  the  key  p'lases.  For  the  choice  P-S,  we  do  not  expect  te 
sec  ancither  low  frequency  S  phase  before  the  key  S  phase,  therefore  this  choice 
uiisuilabh  ,  In  the  c,a.sc  of  P^i-SKS,  tire  amplitude  of  the  key  P  phase  is  cxjrect;  ,1 
to  br,  mu'll  stiioller  than  the  key  .9  jrlia.se.  However,  the  chsorved  btl;ar'ic.ur  is  tl;;  ' 
the  key  /-’phase  is  almost  comparable  with  the  key  S  jrhasc,  therefore  ihi.s  choice 
also  rejected.  For  the  choice  PKP-SI\b\  tlie  angle  of  incidence  of  the  key  P  jrli:’.: 
is  not  suitable  for  <i  PL'P,  as  also  ocen.-red  for  event  B.  Finall}',  tin  rurvir  ing  pli,';  ■ 
jKiir.s  ;it  the  fir.st  getic'i  alion  nodes  are  P-Sh'S  iUK\  PP-FS 

111  step  ■)  of  Talilc  OF  we  note  that  tire  most  likely  intcrpietnt  mu',  aie  for  soin-  ■ 
ateitiici  I  he  surface  01  100  km  depth,  at  a  distance  of  1 12  with  /’/'-/’.S' idciilific.l 
tlie  jiair  ol  le  y  jihn.se.s.  Tin,  l■cvl■l  of  iiiisfil  i.s  \eiy  close  I'e;  ih'  tv,-  ,  e..,.  i  Wl'en 
roniptrre  the  cvcirt  parameter  estimates  with  the  event  sjierilicatiou  in  'J'.rble  G.,F  v. 
r;ui  see  tluit  the  distance  is  very  well  detefmiired.  the  e.stiiuatn  of  deiitli  i.s  re.asomd  1 
gniocl,  while  the  estirmited  a'zirmitli  is  far  from  the  true  value  bccjuise  u!  the  di.’,i;u,  ■ 
llownvcn,  it  i,,  interesting  to  see  tliat  tlie  aziiiuitli  2G1F  for  [phase^.T  in  I-'igurc  ('■  '■ 
(identified  .a.s  PI\[\P'j,']o  very  close  to  tlie  true  a'/iurntli  200'. 

d  he  larger  misfit  for  event  C  compared  with  the  previous  cases  is  ca.tsed  by  tin 
uncnrrelated  big  arrivals  (i.e.,  [jrliase-T)  and  |pliase-9]),  which  most  I'kely  come  fion 
an  overlapjriiig  event  (note  that  there  is  no  obvious  association  with  the  phas-.- 
jrrcdicted  from  the  iasp9J  travel  times). 
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^Fahle  6.8.  Proccsstnrj  sequence,  fo"  Event  C 

Stigi  I:  Choice  oj  Eci;  p/iu.'ir.' 

Key  E:  jihaJir  1,  Key  .S':  phase  -i 

Slcj'  2'  'jest  oil  Dijfcrciilmi  Tim- 


ri75.0  ' 

SUl  vivill;. 

inodols 

:  P-S.  P-SKS.  Pi 

e-SKS.  Pl'-PE.  PEP 

SEES 

S{i  j’  .7. 

r  s  Piiu 

Feat  tin: 

Mate’ 

Survi  vin 

j,  models 

P-ShS 

PP-I’S 

Step  4: 

El  l  lit  Pultcrn  MalcS 

P-S  pail- 

Depth 

DLst 

A'iiin  l'ici'.lil’..::cl  Pliases 

Mi.sFit 

/'  SKS 

0 

45  l>!ias»*  J: 

P,  phase  4.  SKS 

0.8339 

ion 

82 

45  phaMi  1: 

P.  pha.se  4:  SKS 

0.8399 

300 

po 

V  1: 

r,  jihas" 

0.S399 

GOO 

92 

45  phae  1: 

}'.  jihase  4:  SKS 

0.8399 

I'P  PS 

0 

112 

45  |.h.as'.-  1; 

PP.  phase  3:  SKKS. 

phn.'"  4: 

PS.  phase  0;  PKKP 

0.4602 

100 

112 

1.5  plia.-;  1: 

I'P.  pb.ase  3:  SEES. 

pi  KOI'  P 

PS.  ph;-,.-'  i  Eh!' 

0  4G.5.S 

300 

110 

PP.  y'y  <‘i  />' 

0,7101 

GOO 

lOS 

4-5  pli:.-  1 

PP.  i.li.i-.-  4:  /D 

0, ril'd 

Li. 5. 4  D  -  Northern  Bolivia 

I'lic  cpiceiitial  di.stnucc  i.s  even  larger  for  lliis  event  ami  the  scisiiiograni  i.^-  ir''. 
ilominatccl  by  core  phases,  such  as  FEE,  FKEl\  SES,  and  SEE'S  o.lc.  Ihgtiie 
shows  two  complicated  packages  of  care  .irrit'als  in  the  e.aily  part  of  the  broad-bon  : 
record. 

I'lic  phase  detection  systeii'.  is  de.sigMed  to  pick  out  liie  lir.st  onset  of  a  plia.se  pari, 
age.  I'lu'rcfore  the  detections  [phase-2j  and  [pliase-d]  arc  indicaU'd  at  the  bcginnii: 
of  the  two  packages.  hVom.the  low-frequency  detection:-.,  [j)ha.se-2]  i.s  selected  as  ill- 
key  P  pliase  and  [phase-6]  as  the  key  5  phase.  As  before  we  can  follow  the  progress  ol 
the  automatic  system  in  Table  6.9.  Three  choices  foi  the  pair  of  key  phases  surviv.- 
the  test  of  differential  time.  Since  the  incidence  angle  for  [phasc-2]  is  9,78“,  which  i.s 
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'rab,!c  6.9.  I’roccssing  sequence  for  Event  D 

Stiqi  1:  Ckoirr  of  Koj  Phase s 

Key  P\  plia.'^e  2,  Key  S:  plia.sc  C 

stcgi  'A:  Test  on  Differential  'I  tv 

371.5  s  '  .siu  vivin"  models:  P-S.  PP-SKS.  PKP-SES 

Step  I'-S  Pair  Feature  iMiilrl 
.'■■11  vieiiij'  nuitlel.-,;  I’KP  SKS 

Step  4:  Freni  Pullrrn  Match 


f'-S  pair  Depth 

Uisi 

Azi.'ii 

I'lentiliecI  l’has('.s 

Mi.slii. 

PEI'  SKS  0 

114 

U'.() 

liliaso  2: 

J’KP,  jiluise  C: 

SKS 

0.6.1 '2  3 

100 

114 

1 00 

pluwe  2; 

PKP,  pluLSO  0 

SKS 

0.5703 

yv) 

114 

) 

.  0; 

PPr.  pl.a  -  0- 

SKS 

('.,5763 

cm 

128 

10(1 

plia.'ip  2: 

I'KP.  plia.se 

pPKIKI'. 

phase  5; 

sl’KIh'P,  ])lia.si 

•  (>:  SKS 

0  2(.).5.> 

'  1'/  l/-.  i.  2  Jli  liil 

;  ,...i 

tlie  hia;  '. 

■  ('idine  til 

is  ti'i'iuin.'iUii.  Siinili'.vh ,  the  iiv'i-:!- ive  :in"le  aii'i  th-;-  fi enn-ii.-i  d  [].lir'-' i-  ii^e 
.siiil.ahle  i'ui'  a  }‘P  in  (.lie  /’/'-.S7\’,S'  paiv.  'J  l!■'rctllle  (inly  (ine  ehniei'  nl  j’-S  phase  paii'i 
model  l^KI’  SKS  i.s  ’el't,  after  tin;  I’-S  feal'.ne  test. 

'Die  most  likely  interiirclation  idenfil!.'.!  in  .step  4  of 'I'uhle  6.9  is  Ih''  la.st  cluiiis 
of  a  \'(;iy  dce]i  event  at  a  distance  of  IdK'.  Once  afjain  the  laiee  epic  cut  ral  distair  ■ 
niean.s  that  tlie  aziinnth  of  the  liift  /'phase  i.s  not  a  icliahle  csliniat's.  W'lieii  v.s 
compare  with  the  event  spec:ili(alion  in  'fahle  6.5,  (he  eslnnated  depth  is  (|nil' 
accurate.  IIo^^■ever,  the,  estimated  distance  is  9  away  fioin  the  true  \alne,  whii  li 
is  not  as  ();cmd  a.s  in  the  previous  examples,  'i  lii^-  dilference  is  actnall}’  causnl  hy 
the  juocednre  de.scrihed  in  ,section  G.4.3.2,  beicause  the  niajipin";  helwecn  diliereiii  ia! 
time  and  distance  is  not  one-to-one  for  the  case  of  SKS  J’KE  (see  l'i;;nie  G.l).  ^VheIl 
only  one  choice  of  distance  i.s  made  from  the  mapping  (differential  time  to  distance), 
the  true  distance  can  be  missed. 


•  •  • 
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6.6  Discussion 


1 

G.5.0  Pcrformanco  of  tlio  cNcnt  l  ecoguit.ion  systnoi 
As  can  Ijc  seen  from  these  examples,  the  automatic  system  can  .succi’.ssfully  ie.sl:i>  : 
attention  to  a  limited  number  of  possible  interpretations,  aii'.l  then  ;i  c|uantit;;ti\' 
assessment  can  be  made  by  using  the  misfit  measure  between  tlv  ol.i.seivtcl  ai,  ! 
ex])ectod  ])haso.s.  h'he  mi.slil  incasaire  i.ioriorms  very  well  in  solectiii;’,  (In'  nio.--t  ijla.  1;, 
inter[)]etation  and  results  of  tlir  phase  identification  jirocess  give  c.stiniaic.s  for  tli' 
di.staiice  and  depth  of  the  event  wliicii  v.ould  be  very  valuable  h'r  pri  htninarv  c\ei’.' 
loeation 


(•..(i  DISCUSSION 

(i.d.l  Extc'u.sions  of  tho  jn'oco-ssiiig  sc.heino 
'.['Im  present  e\’enl  recognition  sy.'tein  ha.s  been  dc.-igned  to  be  simph,-  iiiui  robirs:  s 
that  \vi'  can  concnit  i  ale  on  the  iroce.^sinc  stratep.y.  W'c  recogni  that,  ;.i  ic’.''  nw 
W"  have  not  inclnd'^d  the  full  ran"e  c!  availalde  f ''ismological  infoir.uit inn  nn  tie 
e\]iecled  cluuacier  of  the  si'isniic  wa'ccli’‘ld-  hor  exauiph',  hu  some  so\iire  laienta 
ti  ..  I  a/ii,i>n'.,.  1.  je,. '  .  /.  I'.i,  d  pth  phases  e.j’,.  pi’,  .s/'an  1  ;ss’,  ,'..S'tn  b’ 

siilestaiilially  lari',er  than  tin-  ('f|invaf  in  fhr.et  phti.ses  and  sn  a  d''plli  plri;-''  coul  1  b- 
pi'da'd  11''  a  key  ]ihase  with  conseiinent  di  lortion  o!  the  pier.eni  pine"  i  lent ilicat in; 
piocess.  I'm tmiat ely,  ;ui  iinpoitant  f"a;me  of  the  assninjttion  lice  ai'pinadi  islli;.’ 
it  i.s  ridatively  easy  to  make  revi;-,ion,s  and  impiuvement.s  in  the  liiLiiii 

I'iistly,  the  seismolnjpial  cxpeili,--  fa  th"'  natuie  iunl  piopeiti.-:  ol  innimoid. 
obneiaed  |)ll.'.se;,  (Ttilile  ti.d)  is  kept  sepaiiUely  from  tlu'  id'enillim  'I  hits  ill- 

addit  icin  m  moditicatimi  ot  tin;  summaiy  of  this  seismologierd  expertise  will  not  alfei  ■ 
the  algmithm. 

More  .si;-mh(  aiitl\,  the  information  leqiiired  tor  coiistrnet ine,  the  assumption  tiee  i 
.also  held  seji.aratefl  fi'oiii  the  algorjlhtn  itself  I-’nr  example,  the  pntjwtli  of  the  .sleni'. 
to  the  first  generation  nndc.s,  representing  the  sot  of  /-’-.S'  jiairs,  is  imiilementei  1 
hy  remline  a  sejtarato  list  of  / ’-5  choices  (d’ahle  0.1).  Tin:  automatic  system  will 
construct  a  first  generation  node  for  everv  item  in  this  F-S  list.  'I'hiis  tlic  inclusion 


6.6  Discussion 


of  further  I’-S  pairs  can  be  readily  a-.xoinplislied  by  adding  iuot  ili■’l[|^  to  tin; 
list  A  similar  apjrrouch  can  be  applied  to  du'  second  generation,  to  irt  ise  o;  exi'  r.  i 
the  [rossihle  source  depths. 

d  he  sampling  in  di.stancc  could,  be-  invrca.'--d  Viv  a  employing  a  dens'  i  travchtii; 
talile  for  each  depth.  'I’lie  current  r.ii!'..h-;:  of  clepth/dist anci;  cinnhinatism 
ha.s  been  chosen  to  allow  rapid  testing  for  a  possible  real  tluie  eiivii oiunent  wliih' 
.still  prm  idmg  a  useful  dehnition  rd  tlu'  i.atu;-'  of  an  e%'ent  foi  jneliinin.'uy  location 
1  lo'.vi'vei .  the  nnnibei  of  classes  coul;!  'j-  increased  without  ;i  very  huge  jirocessiir.' 
overhead  to  provide  an  increased  density  of  states. 


(l.ti.2  I’urtlu'r  (l('veloj)nient> 

till.:.'.!  lu'J'" '■uyi  of  soiirc'(’\=Uini’  nro!  (‘'-nth 

.'\s  (.lesciibcd  in  section  G.d.3.2,  the  sotirco  distance  is  obtained  by  inference  from  th. 
!'-S  ni'i'.l  1  th.e  difh"  -ntial  tin'.'  inti:.:.'.  For  a  tyrS--:.!  r'V.i'.l.n-p, ,  '-.S' plia.-.' 

1  III.'  ,  the  difi.'tetit  lal  tilin'  hni'  an  alni"  '  1;:. '.n  dependein  e  cn  th"  s  . un  ■  disl.aii.' 
’!  In  piiipeiiy  is  exploited  in  tin-  s  iuii.  t  ■  tiinl  th"  closest  nniich  I"  lli  '  iili:."ivi 
dili'eiential  time, 

'I'liei  e  ai  e  tw Cl  factoi's  which  need  to  I,-,  tr.ken  into  consi<!ei  a(  ion ,  lii  s'  1\'  v. e  Iium'  t 
'  .sjus  i  s.iiu"  clili'erenct-.s  between  the  fjb.'-'i vt-.l  difiei riilial  time  ami  th.at  fo:  an  Mail  : 
iiiocle]  such  as  laspflj  and  sti  wo  need  i"  aliow  for  soni"  toh  innco  v.li'T.  niat'liii; 
tiiin's,  .Secondly,  the  mapping  between  differential  time  aii'l  distiuici  is  not  alv.;i'. 
om'-lO'une,  ns  can  be  seen  for  tlie  {I’KI'.  ShS}  atui  {f'Kll\I\  I’h'S]  pairs  in  I'lgiii 
li  1 ;  till.'.  C)C  c.i.sioin.l  [tiobleni  con!.!  V  ;•  :  ..■.'!v••d  by  n.sing  ;i  tlnrd  gi'iniaiiun  in  ih 
:ir.'a.mi})tion  tree  (i'igiue  6  2)  to  determine  distance, 

'1  li"  distance  and  dcijLh  e.stiinate.s  can  f.r  rcline'l  after  t  he  main  event  recognili  ' 
inoce.ss  by  introducing  a  local  assumption  tree  to  search  for  the  best  .solution  in  th 
distance/depth  neighbourhood  of  the  postulated  event  [rarameters  . 
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6. 6.2. 2  Conlnnioxis  operation  and  overlapping  events 
As  we  have  described  in  suclion  0.4.2,  the  jdiri.se.'  roi rosijoiulmj;  to  nii  event  scciuen  ■ 
aie  recognised  and  sejraralcd  from  the  ininil  s'rernn  l.iy  rompaii;;on  wiih  a  cert:;', 
(•.xirecled  jnittr’tn  for  .an  event.  'I'ho  input  is  a  r-jiitinnous  stream  of  jdnise  detecti.e 
so  ihal  we  need  to  also  be  able  to  run  the  event-intnrjneLation  ])r(.Hi’  .>  eimtiiaunr  ■ 
'I  hc  normal  procedure  in  jrattern  recognition  is  use  information  in  time  order  t-. 
.t’.'par.'d.e  eiou]).'.  of  phase:,  a.s  events.  I'vent  by  (  vent,  witliout  going  track  or  reusi:: 
tlie  jnevious  data.  However,  in  the  ca.se  rif  seismic  event  recognition,  taking  ini" 
eonsideration  the  likeliliood  of  arrivals  I'voni  diflerent  events  overlajiping  in  time,  v.- 
would  suggest  tliat  the  rmnnitclmd  arrivals  in  jiiovioiisly  observed  ilata  he  combiii'  1 
v.'illi  tlm  le.st  of  tie'  input  slieam  .ami  reused  1  sr  recognising  an  eta'iit  s-rjiicnce,  i 
veld  tin.’  lajiut  to  the  next  event  process  In  tiiis  way,  we  will  be  iihh  to  .separao 
o\rrla.ni);u!'  events  to  allow  lull  ])hase  u;cogniti.''n  and  inbupietalioi: 


Automatic  seismic  evenl  recognition  using  multiple 
broad-band  stations 


j  liv  auiuiiiaLcd  analysis  oi  scismograuis  provides  useful  infunualion  from  a  singiv 
sO,,-,inogram  bul  is  of  greatest  utility  when  applied  to  a  network  of  stations.  Fo; 
global  monitoring,  the  network  of  stations  is  sparse  and  the  ambient  condition.^ 
eaeli  site  will  vary  so  tliat  tuning  of  signal  detection  thresholds  will  be  required  s  • 
tiiat  compaiablc  operation  can  be  achieved  at  cadi  site. 

.\  1  elated  issue  is  the  robustness  of  the  event  recognition  procedure:  u.siiig  simil:',: 
ti  i.ta  Tills  issue  has  been  addressed  by  applying  the  automated  analysis  procediii' 
described  in  the  previous  chapter  to  multiple  recordings  of  Ihc  same  event  at  clifl'crci:' 
portable  broad-band  stations  in  Australia.  The  stations  cover  a  .significant  sinui  n. 
epicentral  distance  (although  generally  the  azimuths  will  be  similar). 

The  use  of  a  numbe'-  of  differert  recordings  for  the  same,  events  proMcles  ;;  n-  ' 

of  the  consistency  achievable  with  the  automated  analysis  procedure.  'Ihe  portal  1 

sites  cover  a  range  of  geological  conditions  and  ambient  noise  environments  aud 

so  provide  a  useful  test  of  the  reliability  of  the  algorithms.  With  accurate  GP.^ 

timing  and  positioning  the  extraction  of  information  on  range  and  depth  can  b- 

very  effective.  However,  a  limitation  of  the  use  of  portable  sites  is  that  the  accuracy 

of  seismometer  alignment  is  limited  and  so  there  can  be  systematic  (site-dependenl ) 

13:; 
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Table  7,1.  Events  used  for  iliuslration 


Event. 

7'cai 

Day 

Time 

Latitude 

Longitude 

Dejttl; 

.Mil 

A 

1994 

281 

21:44:07.2 

-1.20 

127.98 

17 

0.4 

n 

199-! 

220 

21:08:31.0 

24  72 

95.20 

122 

G.U 

C: 

1991 

244 

15:10:53  0 

40.40 

-125,08 

10 

GO 

D 

1994 

IGO 

00:33:10.2 

-13.84 

-07..50 

031 

7.0 

r: 

1904 

231 

17:20:37.5 

-11.51 

100.45 

142 

0.2 

F 

1994 

232 

04:38:50.5 

44  Go 

149.18 

24 

0.2 

c; 

•  O'',  ' 

18:37:20.0 

44.78 

150.00 

19 

0,1 

11 

1993 

133 

11:59:49.2 

55.18 

-100.46 

32 

0.4 

errors  in  azimuth.  By  varying  the  sets  of  stations  it  is  possible  to  get  a  good  coverag'' 
in  epicentral  distance. 

Figure  7.1  shows  the  stations  which  have  been  used  for  the  events  discussed  in  tin- 
section.  Most  of  the  events  use  the  stations  SCOl-SClO  in  the  Northern  Territory 
of  Australia,  in  conjunction  with,  station  YBOl  in  the  soutl;s:i‘!!.  One  event  uses  th.i- 
S.h  stations  in  Oticensland. 

Figure  7.'.’  shows  tiic  locations  of  events  whicii  have  been  used.  I'he  event.s  l.n\. 
been  ch.o'en  willi  a  minimum  of  foui  recorcl.s  and  normally  at  lea.st  ;i  5  degree  sjeu- 
in  epicentral  distance.  Tlie  events  are  displayed  in  Table  7.1. 

The  first  group  (Events  A-D)  are  those  from  which  indisidual  statioins  have  been 
illustrated  in  the  previous  chapter,  'i'he  second  group  (Invents  E-11)  covers  a  rant;- 
of  epicentral  distances  not  represented  in  the  earlier  set  (30-40',  00-70’,  80-01'  ' 
The  same  set  of  trigger  levels  have  been  used  for  all  stations. 

The  results  of  the  automated  procedures  are  displayed  in  a  set  of  summary  tnhl-  - 
(Tables  7. 2-7.0)  in  which  the  estimated  ranges,  depths  and  azimuths  are  cornpaic-i 
with  tlie  values  determined  from  published  hypocentral  information  derived  from 
global  observations.  Note  that,  the  azimuth  estimate  for  the  key  P  phase  is  used 
as  the  azimuth  of  the  whole  event.  Since  the  estimated  azimuth  actually  implie.- 
two  possible  directions  (with  a  180'  difference  in  between),  both  the  two  possible 
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Table  7.2.  Event  A  -  Depth  17  hm 


StaLion 

Plncsc  pair 

Esi  Range 

£■■1  D._'|iili 

EsL  A'/imutli  Ren-.'- 

.O’/.iinnUi 

SC0,3 

P-S 

PI 

0 

0.  ISO  13 

173 

sens 

P-S 

20 

GDO 

170.  3-00  21 

171 

sc  It) 

P-S 

20 

0 

103.313  2; 

IGc 

PP-SKS 

12S 

li 

YUOl 

P-S 

40 

0 

I.'IS.  3.3S  4-1 

153 

Table  7 

.3.  Event  B 

-  Depth  122  l:m 

Sta  ti-m 

I'liase  pail 

Esl.  Range 

Est  Depth 

Est  Aziiiuith 

Rang'' 

A'/iiiinLli 

SCOl 

I'-S 

62 

100 

141,  321 

62 

133 

SCOl 

P-S 

6-1 

3.-H1 

130,  310 

C-'* 

13  1 

SCO  2 

P-S 

61 

lOO 

133,  313 

132 

SCO.i 

P-S 

6-1 

300 

160,  330 

.53 

13S 

values  for  nziinutii  are  gisen  in  the  Tables  7. 2-7. 9.  In  the  case  when  the  value,'  o; 
til'.'  ‘'luisfit''  nieasMre  (.^ee  chapter  hi  cl-'  not  show  a  clear  jnef''! li'.h.ween  tw  ■ 
inlcipretaticms,  both  of  them  are  iiiclu'.lerl  in  the  suinin.iry  taM-  [s"e  e.g,,  in 
i.e.  i.ii'.i'.  '.l!  *.  I V.  o  mlei  pi  >.  tiiLi^-iis  ivi  staii-'ii  Sv.  It'i. 

■  The  twu  closest  events  arc  lilveiu  A  (Table  7.2 j  and  hvmi  I'.  ('Itihli'  i  : 
v.'liicli  the  phase  pair  used  for  event  thaiactcrisalion  is  /'  -.S'.  I  cu  n.  .-V  tin:  d'-j 
c.stiniale  war-  good  when  an  additional  phase  Pel’  (SCtlh)  o:  .Ve.S  iS’Ilill)  couM  I 
associated  with  the  main  arrival.s.  For  station  SCX)8  the  c.\]rert(d  arrival  time  of  /\  !' 
is  \'cr_v  close  to  that  of  the  S  waves  interacting  with  I  lie  upper  maiille  Iriplicai  idu  - 
and  cannot  be  resolved.  For  SClO  the  discrimination  prorided  by  just  the  F  ■  .' 
differential  time  was  limited  and  the  possil/ilily  arose  of  a  much  larger  ]hcciiir,  ! 
distance  which  can  be  discounted  by  comparison  with  the  other  Elatioiis,  For  c\en' 
F  the  range  estimates  are  quite  good  and  the  event  is  recognised  to  lie  at  dc)itli 
even  though  the  tendency  is  for  too  great  a  depth  to  be  estimated.  The  rather  large 
depth  at  station  SCOG  arises  from  an  attempt  to  associate  PcS  very  close  in  lime 
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Table  7,4.  Event  C  -  Depth  10  km 


.Station 

Pliase  pair 

Est  Range 

Est  Depth 

Est  A'/iniiitli 

Kang. 

A'/ininth 

sco.s 

PP-PS 

112 

100 

-io,  22o 

11-1 

2cn 

SCIO 

PP-PS 

IH 

300 

32.  212 

11-1 

2.37 

SCOl 

PP-SKS 

lOS 

0 

55.  235 

l(i> 

2GS 

P-S 

'12 

100 

SC03 

PP-SKS 

112 

0 

33,  218 

112 

2GG 

I'liis  CVL'IU  .suffers  from  iiilerfcrcnce  fiom  another  event  of  larger  amplitude. 


Tabic  7 

.5.  Event  D 

-  63l  km 

.Station 

Phase  pail 

Est  Range 

Est  Depth 

Est  Asdmuth 

Range 

A'/.iinnth 

SCOO 

I’Kl'-SKS 

1.23 

GOO 

100,  3  10 

137 

213 

SCOS 

PKP-SKKS 

132 

COO 

1G2,  342 

MO 

20S. 

scni 

r!-C''-SKKS 

M3 

101,  311 

1  17 

210 

SCOG 

PKP-SKKS 

131 

100 

1-12.  322 

144 

200 

SCO-! 

I'KP-SKKS 

1.3-1 

iOO 

MG,  32G 

1 

213 

SCOG 

PKP-SKKS 

134 

100 

MG,  32G 

142 

217 

SKS  V 

i.'i'./l*’  ill  { :  1  0 

l-lUi  1  ■  ». 

I.-,;-  1 .  ..ni  lai;:-' 

ennugh  U.  meet  llic  | 

tiig.ger. 

'lablc  7 

,i.j.  Event  E 

-  Depth  14 

2  km 

StiUinn 

Ph.a.se  pair 

E.sl  Range 

Est  Depth 

IP;  A'/.inmlh 

Aviiiinth 

SCOG 

P-S 

32 

300 

SC.  2GG 

30 

2G0 

SCOl 

P-S 

31 

300' 

91,  271 

33 

2G7 

SCOl 

P-S 

31 

101) 

87,  207 

3  1 

203 

SCOS 

PS 

30 

300 

87,  2G7 

3'. 

SCOG 

P-S 

3C 

100 

70,  25G 

30 

258 

SCOO 

P-S 

42 

300 

.S  I,  204 

38 

2G-1 

SCOO  r-s 
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Table  7.7.  Event  F  -  Depth  2.-,'  kmi 


Strition 

PhiisL'  jjaii 

bst. 

Est  Dqith 

E.'ii.  .^.zliiiulli 

Aziinutl, 

SCOl 

P-S 

GO 

100 

10.  19G 

Cl) 

200 

SCOj 

p-a 

C2 

0 

G,  18C 

C.3 

194 

SCO! 

P-.S 

C-1 

100 

0,  ISO 

G.’. 

197 

SC0.3 

P-a 

f>4 

100 

IS,  19S 

01 

201 

SCOG 

P-S 

GO 

luO 

9,  1S9 

GO 

202 

SC03 

P-S 

GS 

0 

7,  1S7 

CO 

19S 

P]\'P-SKK,S 

130 

100 

SCIO 

P-S 

72 

ion 

173, 353 

71 

195 

YBOl 

P-S 

S4 

100 

1G3,  343 

81 

ISO 

Tabic  7 

.8.  Event  G 

-  Depth  19  km 

Slat.ioii 

PIkvsp  pair 

P,st  Paiigi! 

Pst  Depth 

P.st  Azimuth 

liail;.'!' 

Aziiniitli 

SC04 

P-S 

(■■! 

109 

4  1.  224 

C3 

197 

HCOl 

P-S 

GO 

100 

PO.  195 

03 

201 

cr',  •. 

P-S 

I"  ' 

11.  PC 

('■' 

sc:o3 

P-S 

0-4 

100 

19.  199 

G-1 

202 

a  coo 

P-S 

GO 

100 

20.  209 

07 

203 

acoa 

I'-S 

Vu 

i-.-l 

G.  160 

(lO 

199 

PS 

7S 

GOO 

sc  10 

r-s 

72 

100 

I? 2 

71 

1  99 

I'ablc  7.9.  Event  H 

-  Depth  3'. 

?  km 

Stalum 

Plia-SG  |)iur 

Esl  Ilaiige 

E'  1  l)ei)lli 

Pst  Azimuth 

K;ui;',i' 

Azilimlh 

.SA03 

PKP-SKKS 

15-1 

100 

IGG,  31G 

89 

232 

P-SKS 

98 

0 

SAOG 

PKP-SKKS 

154 

100 

21.  201 

89 

220 

P-SKS 

9S 

0 

SAU2 

P-SKS 

98 

100 

27,  207 

SG 

232 

SAOl 

PKP-SKKS 

150 

100 

31.  211 

81 

231 

P-S 

84 

100 
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to  S.  It  Ilia}'  be  desirable  to  modify  the  ynocedu  e  to  try  to  force  a  minimum  tin.- 
separation  between  associated  phases  to  minimise  such  j)roblems. 

Event  B  with  ranges  close  to  50"  gives  a  good  inatcli  of  epiccntral  distance,  th. 
dc’j)lh  estimates  are  mostly  constrained  by  matching  the  expected  lime  of  PcS  (whir: 
might  in  fact  be  SeP  with  somewhat  differoiu  implications  foi  depth]. 

The  two  sliallow  events  E,  G  (Tables  7.7,  7.8),  winch  mostly  lie  between  CO"  an  1 
7fr  fiom  the  sourct:,  arc  generally  well  inaiched  and  there  is  goo.d  consistency  acro=' 
the  sets  of  .stations.  We  note  that  once  again  there  is  a  tendency  for  the  dej^th 
to  be  overestimated.  For  shallow  sources  this  can  occur  when  a  feature  on  th? 
rccoid  i.s  a.ssociated  with  a  near  source  reflection  (e.g.  pP,  c,F)  since  these  caniK'i 
be  correlated  .,  'th  a  surface  source.  Because  of  the  sparse  saniiding  in  depth,  th. 
l)rorcdurc  is  forced  to  make  a  decision  between  0  km  and  100  km.  Gomplcxity  i:. 
the  source  pulse  (in  the  case  cT  these  two  examples,  the  interference  from  the  dosciy 
following  PcF)  can  be  mistaken  for  a  depth  yjhase  and  hence  .suggest  a  deeper  origi;; 
than  aj)i')ro}niaLe.  Also,  in  order  to  match  the  observed  dili'crential  lime  with  tin: 

i'tiiiij.hn^  c>r  th'..  n  1..,  il  r.u.y  ,  ■  t*,.,  iiii|.)i..'Vt  lir  iiulAu  l.'V  a.ljustii.. 

the  dcptli.  I'urther  the  uispOJ  Irav.-l  tiin'i'  ta!  F.s  have  been  used  in  tlie  analysi.'  f  • 
which  .5  is  now  recognised  to  la-  a  iitil-,'  slow  (Kennott.  Fngdahl  iV  Buland,  Ih;'" 
and  this  c.mi  be  compensated  by  increasing  tlie  depth  of  the  souice. 

A  detailed  inspection  of  the  atilornaled  jmocess  for  (■\i  ]il  1'  fnmi  f.m;  slatio;,- 
(SCOl,  SGOG,  SGltl,  and  YBUl)  i.s  presented  in  ITgures  7.M-7.ti  ;mcl  'Itihl'.s  7,10-7.1' 
he  procedure  of  [ihasc-dctectiou  and  chavacteristil ioii  is  lu'pieseiilc.l  iu  the  stu: 
form  ol  lig'tre  with  four  panel.-,  a.s  d■•scribetl  in  rhaplci  G  (seel  ion  G.o).  Etich  of  th’  - 
liguros  is  accompanied  by  a  .summary  table  of  event  inlorpretal  ion  in  llu'  same  fon.. 
a.s  the  tiible.s  in  chapter  G.  The  figures  arid  tables  are  associated  with  the  indivithf: 
.stations  (SCOl  -  Figure  7.3,  Table  7.10,  SCOG  -  Figure  7.4.  'iable  7.1 1 ;  SGIO  -  I'iguu 
7.b,  Table  7.12;  SCOl  -  Figure  7.G.  Table  7.13).  'Fhe  original  seismograms  in  ji.an- 1 
(a)  in  Figures  7. 3-7. -5  show  that  the  waveforttrs  are  very  similar  at  stations  SCO). 
SCOG,  and  SClO.  Every  of  the  three  .statioirs  recorded  a  clear  jrhase  P  followed  b;. 
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7'abie  7.10.  Processing  sequence  jot  Event  F  at  station  SCiOl 

Slc]>  I:  Choice  of  Key  Phufcs 
Key  P-.  phase  1,  Key  S:  phase  G 

Slc]>  2:  Test  on.  Differential  7  ime 

■ISS.O  s  -  surviviiin  models:  P-S.  Ph'P-SKKS 


Stcji  ,7.'  P-S  Pair  Feature  Match 
Sui  viviii"  UKidels:  P-S.  PKP-SFKS 

Step  J,:  Ei'ent  Pattern  Match 


P-S  pair 

Depth 

Dist 

Aziin 

Identified  Pliasc.s 

Misfit 

P  S 

U 

58 

IG  or  19G 

phase  1:  phase  C:  S 

0.515S 

100 

GO 

1C  or  19G 

phase  1:  P,  phase  3:  pP. 

phase  0:  5 

o.isr: 

300 

G4 

IG  or  19G 

pha.se  1;  jihase  3:  }\:P. 

jjliaso  G:  S 

0.2M.S 

000 

GS 

IG  or  19C 

phase  1:  ph.ase  4  pP. 

pl\a.sc  G.  8' 

0.2903 

PKP SKKS 

0 

12G 

10  01  19G 

phase  1:  PKP.  pha.se  O-,  SKKS 

0,5C1G 

Hill 

12G 

10  or  19G 

pha.M'  1:  PKP.  pha,-o  3:  ]iI'Kl\:. 

ph, ■».=''  0:  SKKS 

0.2-tjO 

:’iOii 

120 

10  or  19G 

lihii-  1;  PAT.  idia.-'  0;  SKKS 

0.4S1S 

000 

128 

10  oi  190 

phase  1:  PKI'.  ph.-e-  -I:  pl'KPc,- 

pha.'-e  G:  SKKS 

0.2.:iKi 

.*f.) 
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a  stronger  dcptli  arrival  in  about,  bo  seconds.  The  deplh  aiiisal  conu's  a  lilUc  laK  : 

Uian  expected  from  the  inspOl  travel  time  table.  At  t.he  station  bdlfll.  the  eM.-n’ 

► 

distance  is  somewhat  greater  than  for  the  oilier  three,  riu  l  .signal  to  ii  ji.-j  ratio  ; 

also  lower  (around  1/1  for  the  phase  /')  so  that  the  phase  /'  can  only  be  pic!-''i 

out  by  the  high-frequency  detector.  The  irhase-detection  works  very  well  in  the  foe: 

cases,  and  has  aetually  detected  PKiKP  at  1  101  and  SKiKP  at  YBOl;  these  coi'  ^ 

phases  liave  not  been  identified  by  the  intei])retation  system,  simply  because  they 

have  not  been  included  in  the  preliminary  pattern  (i.e..  Table  6.3). 

The  performance  of  the  automated  procedures  is  very  good  out  to  about  80  °  when  I 
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Fig.  7.6.  Phase-detection  and  feature-e.xtraction  for  event  F  at  the  station  YBOl. 
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Table  T.ll.  Processing  sequence  for  Event  F  at  station  SCOO 

Step  1:  Choice  of  Key  Phases 
Key  P'.  phase  2,  Key  S:  phase  5 

Step  2:  Test  OH  Differciitiat  Time 

C25.5  s  ■  sui  vivini;  luodoLs:  P-S.  PKP-SKKS 


Step  J.  P-S  Pair  Featuie  Match 
Siuviviiig  models:  P-S 

Step  f:  Event  Pattern  Mulch 


P-S  pair 

Depth 

Disl 

Aziin 

Identified  Phase.s 

Misfit 

F  S 

0 

GG 

9  or  189 

phase  2:  P,  plia.se  5:  S 

O.oOVli 

100 

6G 

-9  ( iT89 

fihasc  2;  P,  phase  4:  pP.  phase  .7; 

S  O.2o4:i 

300 

70 

9  or  189 

phase  2:  P,  phase  5.  S 

0.4918 

000 

74 

9  or  189 

l)ha.se  2:  P,  phase  5:  5 

0.491  s 

Table  7.12.  Processing  sequence  for  Event  F  at  station  SClO 

Sie)i  1:  Choice  of  Key  Phases 
Key  P:  phn.'io  2.  K-  y  j.hasc  o 

Step  2.  Test  on  Difjerentiut  'Jinu 

s  !i viviii;.;  moili'ls:  P-S.  PKP-SKKS 

Step  !I:  P-S  Pair  Fculurt  Match 


Surviving 

Step  E 

models; 

P-S 

Event  I’altcrn  Match 

F'-S  pair 

Uejitli 

Disi 

1 

Idei.'.ified  Pliascs 

P  S 

0 

70 

.73  or  353 

phase  2;  P,  pir'..."  5.  S 

(1.59 ; 

100 

72 

J  ’  3  n;  353 

pliase  2;  P.  idias'-  4.  ;'P,  plia.'e  5; 

S  0.23' 

300 

74 

173  cr  353 

pliase  2:  P,  phase  5:  .S 

n.4G' 

GOO 

80 

173  or  353 

phase  2:  P,  phase  5:  5 

0,505 
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Table  7.13,  Processing  sequence  for  Event  F  at  station  )'BOI 

Slc]i  1:  Ctioir.i:  of  Kcij  Pluisc.i 
Key  /':  phase  1.  Key  5;  ]>hase  C 

Step  2:  Tcfl  on  Differential  Tinn 

G13.7  ii  -  suiN'iviiiM  iiiodel.s:  P-S.  f'-fKS.  P,\-.i--SKS.  PKP-SKKS 
Step  3:  P-S  Pair  Feature  Matcl' 


StUA'ivilVJ, 

nuid'ls: 

P-S.  P- 

SKS 

Step  4: 

F.ienl  Pattern.  Match 

P-S  pair 

Deptli 

Dist 

Aziin 

Identified  Phases 

Misfit 

]'  S 

Q 

S2 

1G3  or  343 

phase  1;  P,  phase  C:  S 

0,CC30 

100 

St 

i63  (Ji  343 

p'li.ase  1:  P,  jjhasc  3:  pP,  phase’  G 

;  S  0,2975 

300 

81 

1G3  o;  343 

phase  1:  P,  phase  C:  S 

0.5102 

000 

S-i 

1G3  o:  343 

piiase  1:  P.  pliase  G:  S 

0.5 ',02 

P  SES 

0 

12G 

1C3  o:  343 

p!mse  1:  P.  phruse  C:  SKS 

0.5G93 

100 

12C 

1G3  ov  34? 

plia.se  1:  P,  phase  G:  SKS 

0.5003 

300 

120 

103  o;  3T 

pha.se  1:  P.  phase  G:  SKS 

0,.'102 

GOO 

12S 

1C3  0!  343 

phase  1;  P.  pha.'^o  0:  SKS 

0.5102 

» 

» 
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in  general  the  phase  pair  employed  F-S.  Once  5AT' ovortakcb  6' the  conibinatien 
P-SKS  has  le.ss  clisLauce  and  depth  disrniiiinatiuii,  as  can  be  seen  in  the  rather  jjo  e 
results  for  event  H  (Table  7.9).  Tlii.s  can  he  expected  from  the  rat licr  flat  SEs 
/•’  difl'erential  time  curve  as  displayed  in  Figure  G.l.  As  di.scnsscd  in  tlie  jjicvi':''!-- 
chapter,  the  resolution  could  be  improved  by  chosing  a  number  of  choicf.s  of  distr;:,' ' 
■y  ithin  an  error  tolerance  rather  lhar.  chosing  only  one  in  our  present  program 
At  distances  beyond  90’  such  cornbi nation?  as  PP-}’S  or  }‘P-SI\S  can  often  I- 
useful,  as  can  be  seen  from  the  results  for  event  C  (Table  7.4).  As  indicated  i;. 
thc  previous  chapter,  event  C  sufl'er?  from  uiterference  from  another  e\'ent  of  larg' 
amplitude.  At  the  stations  SCOl  and  SC03,  the  successful  identification  of  the  phar,- 
SS  helps  to  produce  quite  accurate  solutions  of  depth  and  distance.  At  the  station- 
SC08  and  SCIO,  the  low  frecpiency  pha.se  SS  is  contaminated  by  a  high  frequency 
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pliasc  from  tlic  otlier  event.  The  failure  to  identify  £'S causes  a  relatively  poor  depth 
resolution. 

For  the  largest  ranges  useful  coarse  estimates  of  range  can  1«-  ohlaiiied  (event  D  - 
'J'able  7.5)  although  as  noted  in  the  previous  chapter  the  non-inonotonic  dependeiK' 
of  the  din'erenlitil  ti;  te  on  range  tends  to  force  too  small  a  range.  Anothet  contr.- 
buliou  to  the  jjoor  results  for  event  D  comes  from  the  lathei  complicated  phas'- 
jjackages  on  the  .seismograms.  The  situation  of  too  many  phases  interfering  with 
each  other  letids  to  a  relatively  poor  phase  detection,  which  consequently  affects  tin- 
resolution  of  event  recognition.  Improvements  for  phase  detection  in  the  prescnci' 
of  interfering  pliases  will  be  a  target  for  further  development. 

The  various  tests  indicate  chat  in  general  the  performance  of  the  automated  anuh 
ysi,.  lifoceduic  i.s  good.  However,  the  origiiial  set  of  trigger  Icvchi  pr;'ve  to  be  somc- 
v.iiat  c'ouscrsalivo  and  some  clear  arn-ais  can  be  iiiissed  which,  if  iiic.udcd,  won!  i 
improve  the  resolution  of  tlie  seismic  parameters.  A  refinement  of  tlie  eurrcnl  pro 
redure  to  includ"  a  pieliminary  analy.sis  with  standard  trigger  parameters,  followed 
by  a  second  pass  with  tuned  parameters  and  a  liner  sampling  in  eiiicentral  distaiu-. 

'  \  ]  *1'.  v'll  j  t'"''  1’''  P'lv;:*.;* 

Th"  examjvles  al.so  suggest  that,  a  siinultaneous  upi:>licalioii  ol  tin  •  analysis  jir-' 
ii'dme  to  a  seismic  array  wil.h  seveial  stations  could  im|)iovi-  ih"  ac.  inacy  of  t!: 
solution,  'Lhe  correlation  of  the  results  olrtained  horn  diderent  slate.'i..-.  can  a'co;  1 
]jrobloins  arising  at  only  one  or  two  stations.  The  full  power  of  tlu'  array  can  he  ti^-  1 
to  provide  accurate  slowness  and  a'/.imntlL  esliniales  for  eacii  cl  the  ideiit ified  pha,:' 
For  a  cli.strihuted  network,  real-time  event  recognition  can  miu'wed  hy  collectin'., 
rireliniiiiary  rc.suUs  from  each  station,  and  then  all  the  dill'crenl  iiussihle  sohitiun 
(an  be  fed  back  to  every  individual  station,  and  a  final  round  ol  nialching  of  lat''; 
pha.sc.s  cross  all  stations  can  be  msed  to  determine  the  best  .sohilinn.  In  such  a  we.y. 
cflicient  use  can  be  made  of  later  phases,  so  that  the-  numbei  ol  .stations  needed  i'. 
locate  a  seismic  event  can  be  greatly  reduced. 
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AI  techniques  applied  in  the  automatic,  interpietation 
system 


As  I'i  re^’icw  of  tlie  autoinfitic  iatcipretatioii  system  developed  from  cliajilci  4  t.hroi’.f'J; 
to  c.liapter  6,  this  chapter  provides  a  general  description  of  tlio  whole  work  fro:;; 
the  view  point  of  Artificial  Intelligence  (Al).  It  gives  a  syntactic  (IpscrigAion  C'f  a 
seismogram,  and  explains  how  the  AI  techniques  :ue  apjdieil  in  the  iiiiplementati  e. 
of  the  sci.-'moljgical  inlciprclalion  system. 

S.J  OVERVIEW  OF  THE  AUTOMATIC  IN'l’ElU'RF/.l'ATION 
SYSTEM 

liie  whole  system  developed  from  chapter  4  tlnough  to  ciiuiit'';  ii  i,-.  ilosigm''.:  < 

simulate  the  interpretation  by  a  hiinian  scisinologi.st . 

When  an  exirerienced  seismologist.  ii:ads  a  seisinogi ain,  s!ie  would  jirocerxl  in  iv. 

stages:  ]3hase  picking  and  event  interpretation.  Al  the  fir.st  stag",  she  picks  o'.c 

pliasGS  based  on  her  experience.  She  could  feel  the  rlilfereuce  in  lie:  apimai aiic;.' 

jihases,  and  tell  which  phases  are  similar.  This  kind  of  experience  or  skill  can  nc’ 

he  precisely  explained.  It  is  similar  to  the  case  of  speech  recognition,  which  i:-^  ;■ 

task  that  we  all  perform  extremely  well,  but  none  of  us  have  much  idea  how  we  d'  ■ 

it.  Then,  after  the  phases  have  been  picked  out,  she  would  associate  the  different 

1 


S.S  Pattern  recognition  techniques  in  the  work  oj  jihasc  chnractcrisntion 

phases  to  a  particular  seismic  event,  and  give  out  interpretations  on  tlie  locatinn  ol 
the  event  as  well  as  the  identities  of  the  phases.  In  the  second  .stage,  ;i  coiisidci ahl 
amount  of  expert  knowledge  is  applied.  The  expert  knowledge  includes;  frequently 
ob.sorved  ]:)hases  generated  by  an  event  in  a  certain  range  of  distance  and  depth 
dhferent  appearances  of  phases  associate  with  different  identities,  1  he  sei,sinolfwj<: 
might  also  need  to  check  seismic  travel  time  tables  during  reasoning. 

In  the  a  il(.)matic  .sy.stem,  for  the  first  .stage,  we  apply  I’attern  Recognition  tech¬ 
niques.  'I’he  waveform  of  a  phtise  is  recognised  as  a  hierarchical  structure  and  the 
feature.s  characterising  a  jthase  are  extracted  via  a  structural  analy.sis.  For  the  .sec¬ 
ond  stage  of  event  recognition,  wc  apply  Expert  Sj'stem  techniques  and  pa\'  close 
attention  to  the  iin])lrnicntation  dotaifs  of  how  expert  knowledge  is  repre.sLMiteil, 
.icres.sed  (or  organised)  and  applied. 

'J'he  interpretation  proceilures  i)crfonned  by  a  human  .sei.sinologisl  is  compaied 
witli  tliat  used  by  the  designed  automatic  .system  in  Figure  8,1. 

8,'^  I’ATTEIfN  RECOGNITION  'J’ECHNTQIIES  IN  THE  WORK 
CJi  I'liAciij  t.dIARAC  1  ER.tSAi ION 
111  the  work  presinited  in  chapter  4.  a  seismogiam  i.s  an.aly.s'd  a.s  a  hi'  i andiical  sliu- 
line,  'file  techiiiriue  of  sli  net  ural  pattern  recognition  is  also  knov.ii  a.s  ''.syni  a'.'t  i'- 
pat  tern  i  ccognilion" , 

8,2.1  I'lio  Imckgruund  in  syntac.lic  pat  t.uni  rcc.oguil  ion 
1  shall  giye  a  biiel  ilcscrijition  of  the  barkgroimd  to  syiilaclic  patleiii  leiaquiitio; 
based  on  the  textbook  by  lii  (IDfTj. 

In  order  to  repie.senl  the  hieiardiical  structnial  in'orniati.in  of  a  ]iriltcrn,  the  st  ii 
tactic  approach  draws  an  analogy  between  the  structure  of  the  iiattcrns  and  Ih' 
.syntax  of  languages.  Patterns  are  .specified  as  building  up  out  ol  subpattcrns  in 
\arious  ways  of  c,omi)o.sition  just  as  ])hrases  and  sentences  are  built  up  by  concale 
Hating  (or  linking)  chaiacter.s.  Evidently,  for  this  approach  to  be  advantageous,  tin 
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Fig.  8.1 .  Conipari.sou  uf  suisniogi  ain  mtcrprulalioii  i)Voc<;dui  i.‘s  ])crforiiic’ci  by  liuiii.iii  si'lsniologi- : 
and  by  the  de-sigiled  automatic  system. 


S.a  Paitcrn  recognition  techniques  in  the  work  of  phase  charactcnsotion  I",'; 

simi)lcsl  siihpatterns  selcctcid,  called  ‘"patterii  primitive's’',  should  h'’'  much  cus'u  i 
to  recognise  llinn  the  jratLeins  Iheiiiselves. 

The  ‘'langtiage’’  which  jjrovides  the  slnictiual  de.scrij)tioii  of  patterns  in  term- 
of  a  set  of  jmttcrn  priniilivos  and  their  composition  oimrations,  is  called  ’'itattci  ji 
(l(!scrii)tdon  language”.  The  rules  governing  the  composition  of  i)iimiti\'os  int  ■ 
pattern, s  arc  usually  specified  by  the  so-  called  ‘  graininar''  of  the  p.atierri  d('Scri])tioii 
langu;i;te. 

After  c:u'h  ])iiinitivc  within  the  ptiltorii  i.s  i.lentilied,  the  recogiiiliLin  pioccs.s  i- 
aecoinplislicd  b}'  performing  a  syntax  analysis  fpiirsiiipj)  of  the  ])atterii  to  detenniiie 
whoLliei  or  not  it  is  syntactically  (or  grtunmatically)  correct  with  re.s[)ecl  to  tie- 
s)'ecilied  jpaunnitu  .  In  the  mean  .lime,  the  syntax  antdy.sis  also  i)ioiluce  a.  .structura! 
de.srripl  ion  of  the  input  ])attcin. 

d.t.'.  i,  /  }’iiro.sr-strni:iurc  firarnmnr 

Tile  concept:;  of  phra.se- sti  net  me  g.nuunuir--  oiigiiuUed  in  the  p:n.siii;,  o,'  a  siinpl 
Idljdi.-ili  selitem  c,  (.'otusidel  the  lUitidtlg  of  th"  S"iilencc,  ‘’I'lie  lo.se.s  lil-mm  hitaitl 
t'litly’',  !,el  l!:-' r;,'mliol  men.ii  ‘a;-.!'!-  ;  e.  .  1 

l-'ititi'd  by  ] lerformitig,  the  followin'',  Irtuis.  iipti  (or  '■r''wi it ■ " ml--:- 

(sciit etice)  — >  (noun  |)hiiisi.-)  A'eib  plna.e-/ 

(noitn  phrase)  — >  (iulji'Ctive)  (noun) 

(vi'i'li  phrtise)  — >  (va-rh)  (.atlveih'; 

(adject  is'e)  — >  tlu' 

(noun)  :-  i..:.-;, 

(veil))  — >  hlooiii 
(iidverb)  — >  hixnritnilly 

'I  he  slnictiire  aViove  can  be  tihstracled  to  dehtie  the  formalised  jihrase-structuii- 
j'pamniar. 

Definition  2.1  A  phra.s(>.strnc:ture  grammar  6' i.s  a  four  tuple  CV  -  (I'T,  Vp,  J\  S] 


.*1 

m 

■-*) 


» 
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ill  whirl)  and  Vr  are  llic  noiileniiinal  and  terminal  vocalnilaries  of  G,  respei 
tivcly.  }’  is  a  finite  set  of  lr;uisc!i]ilion  rules,  called  ‘'procluetioiis’  ,  denolecl  I  ;, 
o  • — >  /?,  where  a  and  fi  are  strings  over  I’’  (T  =  IV  n  Vr),  and  with  n  involvin;.', 
at  lea,st  one  symbol  of  \V.  S  t  IV 's  Hie  starting  symbol  of  a  senlence. 

So,  in  the  exainjile  sentence  above, 

IV  -  {(sentence),  (noun  jrhiase),  (veil)  phrase),  (acljoctive) ,  (nauii),  (veib  ,. 
(acK'erb) } 

\r  —  {the,  loses,  bloom,  lu,\uriantly} 

S  “  (senlence)  (S.l: 

ami  the  ))i  oduc.tion  set  P  is  jusfthe  ,set  of  transcription  rules  given  ineviovisly, 

( ih'imshy  diN’ided  the  phrese-stmcUire  grannnais  into  four  l\’p''s  aerording  to  the 
toi  ins  of  the  pioductions.  In  type  (•  (unre.stiielc'l)  granunars,  there  is  no  lesliictioi: 
(111  the  prc'diictions,  which  makes  it  loo  general  to  be  useful. 

In  is'pe  1  (coutexL  fien.sil ivc)  j'.iannuars,  the  pioductions  are  resiricled  to  the  loini 

CidV  — >  Ci-bj, 

^\here  .'1  ‘  \',v;  V  1"  (1"  is  the  set  of  all  lmiledi'n,",lli  slum’s  oi  s>’mhi'l'.  ii 

th”  linit  set  of  symhols  including  A,  the  strin;;  nf  leni'.ih  II);  and  /' 

111  lyi'  (coiitexl- flee)  jpaiumais,  the  produ'iions  are  of  tie’  lorm  ‘  /I  •  . 
where  yl  (  \V  iv'xl  P  C  \'  *(V'  '  b*  -■  {A}), 

Itejgtdaf  grammar 

'I'lie  pii (iiluetioiis  of  type  .!  (regulai  oi  liihti  stale)  g.iaimmus  aie  ol  tlm  he  in 
.■1  '  /  uH  ol  A  -  ->  Ij,  when.’  A.Ji  C:  1'),'  and  <:.!■  t  IV.  N<.ile  that  ,d,  h\(iP'  all  an- 
-siii;';le  symhols,  'I'his  is  the  grammai'  we  use  in  ihi,,  papei  to  describe  the  stim  ini' 
of  .seismiogi ;uns. 

S.t’.  ! .!!  Recognition  devices 

A  (p'ammar,  as  introduced  above,  can  he  regarded  as  a  liuite  specilication  ol  a  laii 
guage  fiom  the  viewpoint  of  generating  sentences.  An  alternative  way  of  specifyin;; 
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I  'i;;.  8  i!  Gi  a|  illicivl' rciircni.'iiliitioii  of  tlii' li  Jinsiliiiti  .  o)  “  </; 

a  l;i.ii[';iui[;c  is  in  tinins  of  the  .sot  of  strings  that  are  ace<;i)ted  by  a  certain  jccogniti 
device  (a  lecognilion  vicwjioint),  Here,  we  introdueb  a  simple  reeaigniser,  called 
finite- state  automaton,  which  is  lelcvanl  to  oiir  work  on  pliasc  recogpiition, 
Jffd'M'l  ion  P.P  A  (deteriiibv';',!')  finiti'-stale  .■'utnnialor'.  .■'1  i:.  a  (|uintu])le 

■■  -  (^'  Q,  Co 

in  whit'l' 


1  i  j  M  is  a  finite  set  of  input  .syiiibirls  (alphaliel) 

[ (.J  is  a  finite  set  of  stati’s, 

fli)  h  is  a  inappin;';  of  Q  y  >1  info  0. 

f'l)  (/(I  C  Q  is  Hu:  initi.il  .slate,  and 

(fi)  /•'  t.  O  i.s  the  sol  of  final  sLaLes. 


A  convenient  rejn osentalioii  for  the  inappiii''.  'cAci.  "i 
sitiun  diujtrain  shown  in  ligme  82.  'Die  inlei  ptei  .ii  imh  i  : 
the  antoinaton  in  .stale  qi  and  scanniii;',  the  mpni  svin! 
the  inj)ut  device  inn\'es  one  sti.’])  to  tlie  iie.\t  positi  m  to,-.  ■  ’ 
.symbol). 

Tlieoreui  2.1  Let  G  be  a  regular  grariimat,  llcii  thej.  i  • 
tomaton  A,  that  the  languages  generated  by  (!  can  afw.c,  In 
automaton  .A. 


I  .  lie  sI.lI'  I: 
•ii  i  ■:  .1111  1  I : 


7.  ■ 


1 1 1 1 ,  ■  7  7. ; , 1 !  ■  , 

I  if]'!  'mI  b',  I 
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8.2.2  Syntactical  stiaictiire  for  a  soi.sniof;rani 
In  cliapter  4,  a  seismogram  is  recognised  as  a  seriuencc  o!  ‘  [ilia.se"  cji  ‘'non-jdi;;- 
seginenls,  vrliilc  a  jihase  is  a  strucLure  construrlcd  from  a  nmnlioi  c)f  'units''.  .X-  ■ 
tlial  a  ‘'unit”  of  seisniogi'am  is  basically  a  hill  shaiie  or  valley  sliaju-  .S";’,nif  nt  v.l,. 

)  L'prc.sents  half  a  local  cycle.  'I’hc  hierarchical  siriu  tiiie  ha;,  been  Iniill  lion;  ; , 
holtoin  (pattern  i)rimitives)  to  the  top  (a  coiniilete  jiattein  of  senanogi  am). 

X.i’..'.’.  y  rnllr.i  i}  prnntiifrs  in  n  si  isintH/rinn 
In  the  construction  ol  tlu'  descri|ition  structure  foi  a  .seisinograLii,  six  ])atlein  ]iri;;. 
ilives  have  been  used.  All  the  units  generated  by  an  infidligenl  seginent atinn  jiro.', 
lime  <  ;in  be  cl.assilied  int(i  the  six  groups  of  pat t eru  primitives; 

(1)  I’l,:  a  unit  ol' non  pha;.e  waveform,  i.e.,  backgunnul  noise; 

('.!)  l\,:  the  linsl  unit  wlieie  a  i»h:ise  wavelet  begins; 

(d)  I'l'.  a  null  on  the  left  side  of  a  phase  wavelet, 

(D  a  unit  on  the  li'thl  side  of  a  phase  w:ivelel, 

(.'i)  the  last,  unit  of  a  plnise  wavelet; 

(Co  I’  •  ;,n  1 1 i  r. iliei cut  mtii  elO'.li  ,.nd  il,<- 1  •  ■'  ^  ■  ...i-o  le;  i 

Cl  Miiplel  e  not  mail) . 

inlioduced  in  chanter  .p  phase-delectious  (b.ase.l  i.n  t;'l  ‘\  'bl.\  l!i""ei.'' 

■  out  iinioiisly  j'.eiiei Sited  in  ie;il  t  ime,  which  make.',  availabl"  I  he  infoi  mat  ion  of  w),,  ■ 
tie  I  mient  unit  be;’, ins,  a  ph:i..s('  wavelet.  W'e  us,’'  a  tw(e\';ilue  lla;;,  called  a  ''),li,. 


tia;',"  in  the  proceduie  of  piiniitive  idi’lltili’ at  ion.  Its  initial  wilue  is  'oil", 


indii  ates  no  phase  is  bein<,  :inaly:ed  at  the  curicnt  time.  Whim  :i  phas"  is  indr 
b,\  the  |iha.se-(let,ect,ol  ,  the  phase-llag  is  set  ‘dn",  which  means  a  ]4iase  is  Ir 
analv.sed.  ll  will  be  im  nc’l  ''oil"  later  when  ilu'  cm  rent  phas’'  wa\  el''l  i  .  ieco;’,ni 
a:,  ended. 

1  iussi'  pi iniitive.s  can  be  identilied  a.s  following, 

(1)  When  the  phase-fla;'  is  olf,  the  unit  will  be  recognised  a;;  If.  or  if,  aceordin;, 
the  t:un'enL  phase-detection  value,  if  the  phase-defector  indicate  a  detecLcd  pilia; 


.•.o'vib 


■' Vd.^.V'‘wT, 


"M 
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t.lie  current  unit  is  identifioci  us  U„.  and  itie  pliase-flag  is  turned  “on”;  otlicr\'.iS' 
tlu’  current  unit  is  identilicd  as  L'l. 

(2)  \Vlien  tlie  phase-  flag  is  on,  the  cuiienl  unit  will  be  recognis'^d  as: 

,  if  Uie  amplitude  of  curieni  unit  is  grc.ater  than  the  aiuplituth'  of  former  uni: . 
the  h)rmci  imil  is  ffj  or  atid  tlu-  duration  is  consistr-ni  with 
f'V  ,  if  the  cmrient  amplitude  is  smaller  compared  witli  the  foiinci  one  but  f.reatei 
than  backipouiul  noise,  and  the  duration  is  consistent, 
r,.  ,  if  the  (  urreni  amplitude  is  within  the  anifilitude  range  of  background  nois; . 
f-’„  ,  if  the  duration  of  the  cuireiit  unit  is  not  consislenl  with  i'a', 

Having,  the  pattern  juimilives  ileliiu'd  as  above,  \vc  can  jia\'e  .i  relatively  siiii]ilr 
jsitteni  g.iammar  to  desciibe  the  jiatlern  strueture. 

X.d.. rdttcrn  (imniTiinr  for  a  scifino^mn 
A  regulai  |',ramniai  (!  has  heen  constrnrled  to  'le-wribi'  aay  .s'-ismog,! am  vdlli  nii;. 
niimliei  of  phases  in  it . 

(t'  r-  (r.v.  \'-i  ,  }'.  (.'-'•-isrtiorii-a  rn)) . 
v.heie  (Si-ismugnun)  is  the  start  syinhnl. 


1,-,  {{S(  isinojnidi) ,  (/’/i'l.''/,  {I.'j!;,  ; /I'a, } 

i  ■  I  he  si'l  of  nonterminals; 


t 


iv  -  r,.  I  ,.  f,.} 

i,  lht‘  set  of  leuuinal  symhols; 

.'Mid  till'  ]irnd uetions  are 

l‘  { 


LSeisiiianrain) 

{Scisrnogmrn) 


- >  l'l  {Seisnmgrant), 

-  }  ('„  {I’hast), 


{Phase)  --o  Ui{Lefl), 


» 


► 


w 
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(Phase)  — 1  Ur  (Right), 

(Phase)  — >  Ua  (Seismogram), 

(Phase)  — >  U„  (Seismogram), 

(Left)  U(Left), 

(Left)  — >  Ur  (Right), 

(Left)  — 1  Uu  (Seismogram), 

(Left)  — >  U„  (Seismogram), 

(Right)  — >  Ur  (Right), 

(Right)  — 1  U„  (Seismogram), 

(Right)  — >  Ua  (Seismogram), 

}• 

Note  that  the  pro\-ious  identification  procedures  for  pattern  primitives  are  actually 
context-sensitive,  while  the  pattern  grammar  is  context-free^  Tire  relatively  simple 
context-free  grammar  can  be  easily  recognised. 

S  S.S.3  Pattern  recognition  device 

The  pattern  described  by  the  previous  regular  gr  lar  can  be  recognised  ly 
deterministic  finite-state  automaton  A. 

A  =  (  E,  Q,  (5,  7o>  P  ). 

where, 

E  =  {Us,  Ua,  Ul,  Ur,  Ua,  U,,} 

is  the  set  of  input  symbols  (alphabet); 

Q  ~  {90.  9i.  93} 

IS  the  set  of  states;  the  set  of  transitions  6  is  listed  in  Tab!-*  8.1;  go  is  the  initial 
state;  md  the  set  for  final  states  F  is 

^  —  {90.  9i>  92.  93}- 

The  state  transition  diagram  of  the  automaton  A  is  shown  in  Figure  8.3. 
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Table  8.1.  Transitions  in  the  automaton  A 


{  ''  [go.  th)  =  qii- 

(9ii-  P..)  =  9i- 

(5  (q,.  Ui)  =  q^. 

i5  (91.  Ur)  =  93. 

'1  {q\-  Pa)  =  (J'l- 

<1  (91-  P«)  =  9i). 

<5  ('/2.  P|)  =  <12- 

i  (92-  Ur)  =  93- 

(92,  Pa)  =  9l|. 

^  —  Qii- 

(92-  Ur)  =  93. 

<i  (93-  Pr.)  =  9'). 

(93.  Pu)  =  9o 

) 

Note  1:  a  seismogram  can  be  terminated  at  any  state  among  the  four  choices  in  Q. 

Note  2:  the  four  states  in  Q  actually  correspond  to  the  four  nonterminals  in  the  grammar 
G  which  is  given  in  the  previous  section. 

Note  8;  tlic  two  pattern  primitives,  U„  and  [/„  h;ivc  the  same  effect  of  ending  a  phase 
in  this  preliminary  implementation,  therefore  they  can  be  combined  into  one 
However,  tlie  identification  of  f/„  may  help  to  indicate  the  onset  of  an  interfering 
I)ha.se  in  future  developments  of  the  system. 

P'l  ..  .  '  -,«r.  yn^-n^niliop  system  as  presented  iii  chapter  -1  1ms  been  built  ba.ced  n;; 
the  automaton  described  above. 


8.3  EXPERT  SYSTEM  TECHNIQUES  IN  THE  WORK  OF 
EVENT  RECOGNITION 

A.s  described  in  chajAer  6,  event  recognition  is  a  comj^licatod  problem  which  in 
volve.s  the  application  of  an  incomplete  domain  of  expert  knowledge.  From  the  .A I 
point  of  view,  the  main  problems  in  designing  an  appropriate  expert  system  aif 
(1)  how  to  represent  the  expert  knowledge  (including  experience);  (2)  how  to  a])- 
ply  the  nnordcred  expert  knowledge  during  the  search  for  a  solution.  “Knowing 
what  one  knows,  and  knowing  when  and  how  to  use  it” ,  seems  to  be  an  important 
part  of  expoiti:-;;  this  is  usually  termed  “meta-knowledge”,  i.e.,  knowledge  about 
knowledge. 
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8.3.1  Representation  of  e  pert  knowledge  -  parameters  witli 
adjustable  values 

The  expert  knowledge  collected  in  the  automatic  system  developed  in  chapter  (. 
includes  travel  time  tables,  differential  time  tables,  and  expert  cxpeiience  on  tie 
character  of  seismograms  summarised  in  Table  G.l  -  Table  G.4. 

Some  of  the  knowledge  as  obtained  from  a  human  expert  is  not  immediatel}' 
suitable  for  representation  i-  the  comimter.  such  as  the  term  ‘‘steep  high  freciuenc)' 
[)hase”  (which  is  used  to  describe  the  appearance  of  phase  PcP).  We  introdm  e 
irararnetcrs  to  represent  this  kind  of  knowledge.  For  example,  as  indicated  in  Table 
C'i.2,  we  use  an  angle  of  incidence  parameter  with  a  preliminary  upper  boundary  value 
of  16“  to  define  the  term  "steep’''.  Similarly,  we  introduce  frequency  boundaries  for 
])articalar  phases  (c.g.,  PcP,  PPetc).  And  an  arrival  time  retrieved  from  travel  time 
tables  is  always  translated  into  a  time  range  with  a  tolerance  for  small  difference. 
The  values  for  these  parameters  may  need  to  be  adjusted  during  te.'^ting  to  ensure  a 
reliable  performance. 

8.3.2  Kiiowlcdgo  application  in  the  procedure  of  seismic  event 
recognition 

The  procedure  of  knowledge  application  is  a  procedure  of  reasoning,  It  is  a  kind  of 
engine  in  an  expert  ,system  -  called  an  ‘‘inference  engine"  in  the  field  of  Artificial 
Intelligence  (Al). 

Many  expert  system  shells  (software  tools  svhich  help  to  construct  expert  systems) 
u.sc  the  architect  nrc  of  a  rule-based  system  (see  e.g.  Ford,  while  rule-based 

.systems  use  decision  trees  to  iiirplenrent  their  inference  engines. 

It  would  be  very  difficult  to  apply  this  kind  of  architecture  to  our  present  appli¬ 
cation  of  seismic  event  interpretation.  Because  at  the  present  time,  we  don’t  have 
a  precise  and  comjjlete  rule  set  and  the  limited  domain  of  knowledge  is  subject  to 
continuous  modifications  and  improvements  due  to  tire  coirtinuous  development  of 
seisrnological  expertise  with  broad-band  data.  The  identification  of  seismic  phasc.s 
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is  by  no  means  a  trivial  exercise.  In  fact,  many  modern-day  scisniologi.sts  have  lii- 
tle  direct  experience  in  tile  routine  'reading''  of  seismograms.  Mveii  experienced 
seismologists  can  rnisidentify  arrivals,  given  the  many  possibilitie.s.  d'hc  ISC  (Intei- 
national  Seismological  Centre)  often  leidentifies  phases  picked  by  station  operatois 
who  do  not  have  accurate  location  cstimalos. 

Since  the  unknown  factors  (origin  time,  distance,  depth)  are  interacting  with  each 
other  through  the  exj^ert  knowledge  in  a  complicated  way  and  the  available  expert 
knowledge  is  limited,  using  a  normal  rule-based  architecture,  the  system  could  be 
bogged  down  in  endless  searching. 

To  solve  the  problent  of  knowledge  application  in  seismic  event  interpretation,  a 
new  tcchniciue  of  “as.sumption  ti'ee”  is  proposed  in  chapter  G.  It  was  ins])ired  from  a 
typical  human  research  strategy  -  making  assumptions  and  then  using  l!;e  inferences 
arising  from  the  assumptions.  It  is  an  attempt  to  make  an  cHirienl  apj^lication  of  a 
limited  domain  of  expert  knowledge  while  there  are  a  number  of  interact  ing  unknown 
features  to  be  constrained, 

In  the  specific  appliciUion  of  seismic  event  interpretation,  a  tlirt'c  hwi’l  assumption 
tree  has  been  used  to  guide  the  apjilication  of  expert  knowlerlge. 

At  the  first  stage,  the  tree  grows  out  nine  brandies  which  repre.scnl  nine  dill'ereiil 
assumiitions  on  the  feature  of  ‘T  S"  [lair.  Then  the  knowlcrlge  liasc  which  onl\' 
depends  on  the  “P-S”  choice  (represented  by  Table  6.1  and  d'ahle  (').‘2)  can  be  applied. 
The  inferred  result  ,  a  constraint  on  some  other  features,  is  used  to  reject  tlie  ciii'ient 
choice  or  to  constrain  choices  for  later  level  unknown  leatures.  lii  other  words,  the 
result  of  the  application  of  knowledge  base  is  used  to  prevoiil  the  tree  from  growing 
I  to  its  full  size. 

Similarly,  at  the  second  stage,  every  surviving  node  from  the  iirsl  lc\el  braiichc.s 
into  four  child  nodes,  each  wit.h  a  different  choice  of  dejith.  h'lien,  know  ledge  whirl i 
depends  on  “P-S”  choice  and  dej^th  choice,  represented  by  differential  time  tables, 
can  be  applied  to  constrain  the  distance. 

When  the  distance  is  known,  the  knoivledge  which  deiiends  on  distance,  repre- 
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sented  by  T?iblc  6.1,  Table  6.2,  and  travel  time  tables,  can  be  ai:>pliccl.  And  tlii 
result  of  the  knowledge-application  is  used  to  measure  the  tilness  of  the  cuneni 
interpretation. 

The  separation  of  domain  knowledge  from  the  inference  engine  makes  it  eas,\  t" 
update  the  expert  knowledge  in  the  future.  The  entities  in  the  seismologicnl  e.xpei  i 
-system  and  their  relationships  are  illustrated  in  Figure  8.4. 

The  domain  knowledge  is  partitioned  into  six  groups;  Table  6-1  --  'Fable  (i-  l. 
travel  time  tables  and  differential  time  tables,  which  reirresent  distinctive  spccialisis 
(or  independent  knowledge  sources).  The  application  of  these  knowledge  source.s  is 
scheduled  by  the  assumption  tree.  Solutions  are  built  up  on  a  global  data  struc- 
Inrc,  the  “information  package”,  which  rejjresents  a  tyi)e  of  ‘'blackboaiil"  (wlncli 
is  a  technique  in  the  field  of  expert  system).  This  particular  type  of  blackboard 
-  the  information  paclcage  -  has  multiple  existences  during  pixicessing.  'Fhey  arc 
dynamically  crisUed  and  evoluteu  wliii  ihu  uesr  l^pnieiil  of  ih.;  a.i.smnpiion  tree.  In 
fact,  the  assumption  tree  method  combines  the  blackboaid  and  ri'.isoning  sclicdulci 
into  an  integral  tree  structure. 

T'l’v,  1  V  mnl;!’!'':'  r  •  and  rix  nn;:;',  kunwled''"  r>"C(VvU,,,,.  the;: 

premises  (or  conditions),  the  assumption  tree  method  has  cnal>led  ellitient  know] 
edge  application  in  the  complicated  ])roblcm  of  seismic  cvcnl  intcipr;'!  ation. 


8.3.3  A  c.oniplcrnoiitary  clesc.rijrtioii  of  the  assuiniJtion  trc('  im't.hod 

'Fhc  assumption  tree  is  a  new  automatic  reasoning  method  which  i.s  presented  in 
this  the.sis.  Although  proposed  for  the  design  of  a  .scisinologiral  expert  .sy.slr:ni,  it  i  ■ 
a  general  ajrjjroach  to  an  efficient  automatic  reasoning  with  limilcd  iloinain  e.\))'  ii 
knowledge, 

I'his  section  gives  a  generalised  description  of  the  assumption  tree  method  ami 
some  important  design  criteria  when  applying  this  mctlrod.  The  material  presented 
may  guide  to  design  a  new  type  of  expert  system  shell  -  “an  assumption  tree  exj-'ci  1 
system  shell”  -  for  a  wide  range  of  expert-system  applications. 


» 


Fig.  ti.4.  Entities  and  their  lelationsliips  in  tlie  seisniological  expert  system. 
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8.3.3. 1  Suitable  applications  of  the  assumption  tree  method 
The  applicution  [)roblenis  which  fall  under  the  umbrella  of  the  assuinpLiou  ti'- 
method  have  the  following  common  description:  given  a  collection  of  expert  knowl¬ 
edge  and  a  partial  description  of  an  unknown  object,  to  constrain  the  choices  !'■: 
the  object,  The  description”  refers  to  a  constraint  on  itart  of  the  feutui' 

which  describe  the  object.  The  expert  knowledge  which  is  available  is  a  collection 
of  jjroccdures  oi  rule.s,  each  of  which  produces  a  con.straint  on  a  .sub-'ioL  of  feattii'  - 
following  from  a  constraint  on  another  subset  of  features. 

'I'o  helj)  later  discussions,  denote  all  the  features  which  describe  the  object  a.s 

d'ho  nuinber  of  choicr;s  for  F,  i.s  denoted  as  n,,  and  all  the  possible  choice.s  lor  /■',  ran 
li','  icpieKcnted  as 

il,  —  {/'i.li/i.'j!i/i.3i''i/i.ti,}' 

8. 3. 3. 3  The.  assumption  tree:  structuix  -uu!  rcusoniny  ulgorithni 
The  .slructure  of  an  a.ssiunption  tree  is  given  in  I'ijpire  (h.'i.  bvery  I'catuie  (or  fact"; 
corre.siionds  to  a  si)ccilic  level  (or  generation)  in  the  nsf;unii)ii()ii  tree.  IJraiichin. 
of  the  tree  ul  ;i  specilic  i(Wel  is  iinule  according,  to  the  choices  ol  tlic  eoi re.spuiidiii 
feature. 

Knowledge  application  is  organi.sed  according;  to  the  lev(!l  in  the  a.ssuniption  ti" 
.Al  each  level,  a  diHereiit.  set  of  ox])ert  knowle(lg,e  is  lelrieveil,  As  ordeu'd  from  neC 
to  leaves  in  Figure  G.3  (/’i,  Fx,  ...,  f'm).  fhc  domain  knowledge  items  is  divided  in'- 
in  gi(juj).s  aecorthng  to  the  order  they  will  be  a|Ji)lierl  : 

(5'])  —  knowledge  items  which  only  lely  on  the  choice  of  /-j, 

{So)  -=  (knowledge  items  which  rely  on  i'j  and  F->)  -  .S'l, 

(bj)  —  (knowledge  items  which  rely  on  /'j,  Fo  and  F^)  -  So  -  b'j, 

(bm)  —  (knowledge  item.s  which  rely  on  7A,  7-2,  ...,  and  F,„]  -  Si  -  Sx  -  -  S,,,  i. 

At  the  ith  level  of  the  assumption  tree,  knowledge  set  S,  is  applied. 


’I 

1 
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Every  node  in  the  assumption  tree  carries  an  information  package  (cienotecl  as  ./ 
wliich  records  tlie  current  constraints  on  all  features: 

I  -  {luh.h,-  ! . A-,},  /.  C 

The  reasoning  algorithm  begins  from  the  root  node.  The  informal  ion  jracku,. 
at  the  root  may  contain  a  constraint  on  some  of  the  features  (which  re]MC,si‘iii> 
individual  evidence),  and  full  ranges  of  choices  for  the  other  features.  The  loci 
node  branches  into  ?ii  child  node.s  according  to  the.  number  of  choices  for  the  fii>i 
level  feature  Fj  (which  is  recorded  in  the  information  package),  and  tlu'  informat im 
package  ol  the  root  is  copied  into  every  child  node.  Each  child  node  has  a  different 
choice  for  l'\  (denoted  as  /i.,),  therefore,  in  its  infonaalion  ]!ack:r'',r  the  constraim 
on  f’l  is  updated  with  i  he  otrly  choice  of  fi  ,.  'j’hcn,  every  node  searches  in  tl'.r 
knowledge  set  .Si  lor  ap]ilicable  items  from  the  knowledge  base.  The  ap]rlicalion  cl 
an  item  of  knowledge  will  produce  new  constraints  on  some  other  feature.s.  I'or  a 
new  eoiistrainl  for  the  feature  F,,  denoted  a.s  when  //  O  /,■  is  an  einiity  set.  a 
contradiction  occurs,  and  the.  growth  of  the  euneiil  node  is  lei'ininated;  otherwise. 
Ill"  old  ci.v.islraiiit  /,  is  replaced  with  / ' /,.  Afou'  the  procedure  of  Ihe  ap))licat imi 
of  knowledge,  a  surviving  node  eontinues  to  hrancli  according  to  tlie  choic('s  fur  llie 
next  level  feature,  ami  the  inoeedure  will  continue  voeursively  unlil  the  current  iind' 
is  rejeeted  or  a  loaf  of  the  assumption  tree  is  reached. 

8.S.‘1.S  How  to  arrange  jmtnrcs  in  an  (issninjilinn  tree 

While  applying  the  assrun];lion  tree  method  lo  a  general  ai)])lication  ofexjrcrL  system 
design,  a  c|uo.stion  arises:  how  to  place  the  in  fealures  on  dilleriuil  Icrels  in  ill' 
assumjition  tree  so  that  it  can  work  most  ellieiently  V  Thai  is,  which  ])ermulatiun 
to  choose  from  m\  possible  choices  ? 

Suppose:  (l)  all  the  knowledge  items  will  be  used  to  terminate  a  current  node. 

(2)  the  order  of  Fi,  /'2,  . ,  F,,,  is  chosen  as  ihe  order  from  root  to  leaves  in  tlie 

assumption  tree. 

For  an  item  of  knowledge,  the  earlier  it  is  applied,  the  more  efficient  it  is  used. 
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For  exainjjk;,  application  of  an  item  of  knowledge  in  the  first  generation  will  rcji'ci 

rin  X  713  ^  ^  .  ^  solutions  (note:  n,  is  the  number  ot  choices  for  the  featiiir 

F,),  while  application  of  an  item  in  the  last  generation  can  only  reject  one  solution. 
Thus,  the  number  of  choices  for  every  feature  shall  be  taken  into  consideration.  The 
feature  with  more  choices  might  bettor  be  ])laced  at  a  later  generation. 

Since  all  knowledge  items  comirete  with  each  other  for  earliei  a]5pli(:ation,  we  shall 
take  the  order  IVom  the  earliest  generation  to  the  latest  while  arranging  the  features 
on  the  assumption  tree  to  achieve  the  Irest  integral  efficiency. 

As  de.srribed  previcu.sly,  the  knowledge  set  which  is  aiiplied  at  the  Uh  level  i.s  .S',. 
Let  the  number  of  knowledge  items  in  5,  be  k,,  and  k,  rejire.scnls  the  upper  limit  o( 
the  number  of  kuuwletlge  items  which  could  be  aiijiliod  at  the  ith  generation. 

As  mentiimed  above,  iirslly  we  shall  choose  a  feature  (from  the  in  features)  loi 
the  first  generation.  The  number  of  solutions  which  would  be  eliininati'd  at  the  first 
j‘eiieratiun  can  be  measured  as 

f?i  “  t/3  y  m  X . >:  u,„  x 

The  irr<'''tr'r  the  value  for  ('ll.  the  I'igher  the  (juality  of  the  assumid ion  tree. 

Fur  the  second  generation,  the  actual  nuinbei  of  knowledge  items  apjilied  at  nm 
liui"  would  depend  on  indivirlual  applications.  'I'liis  is  beeam.e  that,  the  first  gjiinei 
alioii  plicat.ioius  might  have  f.:ut  off  some  branches  so  that  the  second  generation 
would  not  grow  t.o  full  size.  However,  we  can  still  use  tlu'  upjiei  limit  number  k',  to 
evaluate  the  average  efficiency.  The  elliciency  of  the  second  generation  can  still  be 
measui'ed  as 

Q.^  —  7(;|  X  7l.|  X  7?r,  >:  .  X  7t,„  X 

(deiu’rally,  the  quality  of  thi:  rth  generation  is  given  by 

Q,  ~  7i,^  j  X  'ii3  2  Ui+n  X .  X  n„,  >:  A,. 

VVe  can  use  this  cprality  measure  to  choose  apjrropriatc  features  for  all  the  levels  ol 
the  assumption  tree  from  root  to  leaves. 
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8.3.!iJi  The  choices  for  the  ‘'assumption  features” 


Tlie  features  whicli  appear  as  sjrecific  levels  in  an  assumption  tree,  naiiu’d  "assunip 
lion  featAires",  arc  not  limited  only  to  the  apparent  features  wliieh  describe  tli>‘ 
object,,  they  can  also  include  some  redundant  classification  features,  uliich  C'UiM 
sometime  greatly  iupirove  the  performance  of  knowledge  applicalioii  I'oi  examph  , 
in  the  previous  seisinological  interjiretation  system,  the  classification  ol'  ‘d'-S  "  pair 
models  is  a  redundant  classification  which  is  a  comliiuation  of  a  dislaiu  e  const  raiiii 
and  two  constraints  on  idiase-identitj'.  When  some  combinations  of  constraints  oii 
original  features  fre(|uently  ajrpear  in  the  jneinise  part  of  some  knowledge  items, 
t  hey  may  suggest  ;i,  good  choice  of  a  redundant  elajisilicalioii  feat  m  e,  W'ln'tlu  i  oi' 
not  to  include  this  kind  of  choice  into  the  structure  of  the  assumpl.ioii  t.iis'  can  be 
decirled  by  using  the  rjuallty  measure  which  is  described  above. 

S.3, 3. 5  The  'iivw  idea  in  the.  assntnption  tree  method 

Although  the  tissumpliou  tree  me.thod  luis  been  designwl  lur  a  parlienlar  i',i'oU))  ol 
;ipplieatiuiis,  we  can  cuiiiptuc  it  with  ollu'r  gymeral  reasoning  inelhods  in  a  chaiai 
ti’i'i.stic  a.simet, 

'I'he  lutrjor  new  ide.a  in  the  assumption  tree  method  is  to  schedule  icasoning.  ac 
coixliug  to  the  expert  knowledge  which  is  avtiilable. 

hi  the  held  of  expert  sy.stem,  tyjiical  examples  of  the  control  iiieehiiiiism  loi  ma 
soiling  are  “backward  chaiiiiMf;’'  (or  “backward  leasoniiig")  and  “I'oiwaid  clmiiiin;',' 
(or  “forward  rea.soning“ )  (see  e.g.  Jackson,  IDfiO).  Ihu'kwaid  (  liaining  entails  taking, 
a  conclusion  first,  then  seeking  evidence  to  suiiport  if.  Forwiud  chainiii;’;  works  in 
(lie  reverse  (.lirection,  t.lie  system  tries  to  find  a  rule  whose  comlilion.s  are  now  sal 
islied  by  t  he  inforiiiatiuii.  In  the  case  of  iiisullicient  evidence  and  iiisullleieiil  expei  l 
knowledge  while  the  choices  for  the  solution  are  many,  using  the  hackwaril  chain¬ 
ing  method  the  reasoning  procedure  may  never  terminate;  while  using  the  forward 
chaining  the  reasoning  inocedure  might  not  he  able  to  progress  to  any  solntiuii  (due 
to  iusiiflicieiit  evidence). 
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In  contrast,  in  the  assumption  tree  method,  the  reasoning  piorcduie  is  rhaim  ! 


neilliei  from  evidence  to  solution  nor  tire  other  way  around,  il  is  dii'H'ted  in  a  jiii 
decided  order  which  has  taken  into  consideration  tlic  distrilnition  of  knowledge  v  iiii 
the  features  (as  described  in  section  8. 3. 3. 3).  The  new  idea  is  to  follow  from  ‘  whai 
one  knows”  to  decide  “when  tand  linw  to  use  it”.  In  llii.s  way,  the  assumption  iri' 
method  Could  ai)))ly  knowleclge  more  eliicienlly  in  the  c:ase  of  insullici(uit  knowl 
edge  and  evidence,  and  the  rea.soning  |)rncedure  is  guaranteed  to  terminate  in  an\ 
circumstance. 


•  •  •  • 


Bil.)liogr;iphy 


1  '  1  .  \v.>vc  i'i:i  i;i  I  ’.vi':  "'I  j'jiisotro;'!-'  I’V'li:'.,  J.  O' /,'■  ■ 

(lu,  ‘Jon;! 

Aiuli'ii.iia,  1),  l...  mS'.).  '1  lii'oiy  ot’  Lhi?  I'jiirlli,  IMivrk'vcIl  Si'i<;ui.ilu;  rublii  aLimi.-,  Hosl-on,  '5G(i  I'j' 
Aiiili-rhiin,  1),  1..  A-  U/ii'woiiski.  A.  iM.,  1082.  Uppci  -iinuiiU!  aiiisiitnipy:  iioiii  ii  i"' 

(im'illiit ion;;,  (iriipliiis,  ./.  Jloi/nl  Axtroii.  Sor..  00.  .')88-'ll).l, 

Amlii.  M.,  I.sliikaWii  Vaiiiivy.nki,  k’.,  1083.  Shcar-wavi'  polaiizatioii  .•uii-iol  lopy  in  (In'  iipp  : 

iilaliUi'  |i"iii'at.h  llcilislin,  ,la|)ati.  ./.  A'lOp/iiyi.  A'r.'..  08.  OSOO  1 
nimla"',  A.  A’  Mil  I'linovaa,  10(11.  An  appli'Ml.iiiii  ol Atalislical  (iiaiiiiiinalidii  in  clissily 
sni.'iiiiic  vvciiU.  Hull,  sfixtii .  .V()( .  Am.,  fi-1.  901-971 
llnwiii.ui.  .1.  )(.,  A  Amin,  M.,  1087.  flhnar-wavc  !,plii;i:..-  m  llm  upp'M  .man' I"  aluivn  ll; 

'rnlij'a  ,'ail)i lucl.ioli  Ziini',  ,!.  ]i.  Artivii.  8H,  '.!.n  A'‘. 

Ill  I'iiiiaii,  1,,,  kVii'diiiaa,  J.,  Olsliaii,  H.,  A’.  SUinn.  (7.  lO'Sl.  (  'hi.'i.'iijiciilion  v.iiil  (n  i.v 

lli’lliinnl.,  (.lA:  Wadswui  Lli  liiLaniatiimal  j'.iiiiip. 

(Ici'vdiiy,  V.,  Jvlulolknv,  I.  A.  A:  I’.sam'ak.  I..  1077.  Kay  iiirllicxl  in  ;;ni-;innlia>y,  Univi'i  /it.a  Kail 
I’raha, 

('■Ill'll,  (.7  II..  li)82.  Aiililicalinii  (ll  paU,inn  raiaij'iiiLinli  In  Miisiuic;  wav*' iiil  I'l  prcLal.iiiii.  A/>/ilii'i:: 

()/  jialLcvn  IdV  110,  liuc.a  llatnii.  I'la  :  CIUI  I’nis.i. 

C Ilii'isU'iiai'ii,  N'.  I.,  A'  Salisbui’y,  M.  11  ,  1070.  Soiaiiic  aiiisnlropy  in  lln'  upper  iinuilli':  Mvidei.- 
Il'nin  l.lU'  Day  ni  Islands  ()|ilnolil<i  cninplex.  ■).  (Ou/I.'ii/X.  Jtc.i.,  8  1  'KUII  'IGld 
(Halin',  07.1.,  100.3.  'riii;  cuniiiliitLi  oiilnmil  lay  expan.-inn  -  1  C.lalenlalinn  uf  .synl.belii:  snisnU),",!';-. 
Cr'ra/iAy.i.  J.  hit..  115,  421-43-1. 

Dcy,  S.C.,  KennnI.l.  D.L.N.,  Dowinan,  J.ll.  &:  Goody.  A..  1003.  Vaiiatioiis  in  the  upper  inantl' 
Velocity  stnicLurc  under  nui  tbern  Australia.  Geophyx.  J.  hit..  114.  3()1-31(), 

Unity,  T.,  &  Aiidei'suii,  L).  L.,  1980.  Si'isniic  vel()citie.s  in  mantle  minerals  and  the  mineralogy 


Ilibliograpliy 


IT 


till'  vippi'l  lIKUlM'',  J.  Gi'iipbys,  Urs..  94.  ISO-T  1012 
1  A.  M..  i'.',  AiuI'Msdii,  D.  L.  lO'l.  i  rli'i  cm  •  .1|  l  li  iiiuiliT  Uni.',.  Hurt!. 

I’htinl..  Inin..  2,0. 

l',:u  If.  r.  S.  .k',  SIc'aici,  I’.  M..  1001.  i  I'iil  i'i:i  (,1  cloli..  .ini'-  ii'ii,'.  c.ii 

.il  ic-|iiiUiii_c  ;ilL',iH'illilu,  Hull  S’l'i:  Sin.  An-.  8'1.  OOG-.'iTCi. 
l-'oiil,  N..  lOOl.  Expert  .sysL(MU'i  ;iii(l  iiit ill- i.;l  ii.'i'Hivri.'iK  c.  Liluaiy  As-im  i.ii uui  Li.;. 

111.  1\.  S.,  lOV-1  Sy  lit  ail.i'  meil.mU  in  pa’i'iii  i,  •  n.nnitior..  2‘.tj  pp..  Aca.l,  me  I’l-  N'-w  'i  i  ; 
anil  LdiuIihi 

I’luTi'i,  1\  .  lOSO.  Ilceeiitly  liinui’il  clastic  ai.i-'p'.mpy  ;iii‘l  pctiolnci'  al  inmliT-i  Im  tie'  einitiiu'ii'a!. 

siilii  1  u-,;  ill  111  liii.-pliere  ill  sinil  liein  (icnna'.'i.  I'lnj-  l.iirllt  I'tiiiu  l.  Inh  r..  111.  03-118. 
l-'ukai).  y ..  10,S  1.  SeS  t’viileiirc  lin  aiiisut  I'lipy  in  ill'-  i-aitl:  s  mantle.  iSnldti .  300.  C0j-G0,S. 
tlnclze,  C.,  V  Kolil.st eilt,,  1.1.  L.,  1073.  Ealun'atoiy  stmly  of  (lislncatjim  cliiiili  timl  clilliision  in 
olivine,  i;,s.  [ir.i..  78.  3001-3071. 

Cliiinl\,  .A.,  1001.  1  li  uail  liaiul  si  iiilics  ol  lli' upper  inanilc  licncalh  nortlieni  Australia,  honnui-.,, 
t/ie.M.s,  Aiiel  Katl.  Lhiiv.,  CanlifiTa 

( liulniiiiiilssiin.  ().,  Keniu't.t..  Il.E.N.,  .k-  (loo'ly.  A..  1001.  Hi  ii.nl  Ii.iikI  iiliserv.it.iiin-i  ot  upper 
niiiiii.l'-  .sei-.iuie  phases  in  norUiern  AH'tialia  anil  the  at lenmrt imi  slnici.iire  in  the  ni>]>ei 
inaiith'.  yVii/.s.  Hiirlli  I'lnni'i.  Inter..  8-1 

1 1  ay.ei ,  11.  1  i , ,  ,k'  0'(  amiiell.  li . ,  1070.  Kin"mai  i'l  nnnlels  nl  l.iie,e-.--eale  I  lew  in  t  h"  i  I’.rt  h'.s  inai','.  h  . 
./.  (Ir.iplni,.  /I'p-.s,,  84.  1031-1(118. 

Hearn,  'I  .  lOM.  I’li  tnu-el  limes  in  ,Suui.iie;n  ^..■lh.;^llm.>  ./  liT n/i/.-;,--  IS.-.,  8!),  1-^  l,'l-lS.i),T 
1 1  em  Ira  java,  A  .1 ,  .  1081.  A  st  I  lily  of  S-l>oil.'  v.-ave  veil)'  i’ y  si  i  iieiin  ■'  in  lie  mail!  1"  ihiwn  to  1 1  'M 
l.iii.  II -.ii,,  lie'  \ .  all  animi;.,.i  I- ismic  a; ; .  i’li.U.  the.si.-,.  Aii-i.  Hall.  Ini.  Caiiherra, 

Ih'--,  11..  I'.liil.  hi'i.-iiiii- aiiisoi.iopy  111  the  i;p;e-iim>si  m-i:.' h' mi'le;  i.c.  aii-.  .'.r.ti  203,  020 

1  Ini  111'.  ii  r,  S  I,  1077  1 ’eak  I  ei'oi'iiil  imi  ill  V.  e.-.i'ln;  III-- .  .S  i.'i-P/i //I-  p  .  ,'o  1 1  i,,n;i  .  -.  Ajijiiii'nli'n 
31-10,  Sprili;'/'!- T'erlafi  llerliii  lloidelher,.  New  Vii:’.. 

.laeksiiii.  1',.  1000.  liuroiliiittioii  to  o.Kpei  t  ry.-t'.'iiis,  Oml  eih.  A'lihsi'ii  Wi- .h-y 

.hiuriiet.,  11..  .V  .liilii'i  1..  N..  1082.  Varialiiiii  wi’h  a.iv  nl  .'.'.lisi'lriipy  iin'l''i  m  eaii  .  hniii  i';real.-i  i; i ! 

.'iiilfar'i'  waves.  Ctrojilii/.i.  /I'c.c.  I,r.ll..  9.  17‘,l-l,Sl. 
l\'aral.(i,  S,.  S:  Li.  1’,.  1002.  l.lilinsive  rreep  in  pi'i la.-kn e;  Iniphrnt mil  ■  hn  I Im  i  le  "Ini’.y  iil  tie 
liiv.er  iiianLle.  Srirriri  .  25>''i,  123S  12ln 

Eenni'tl,.  M  1,.N,.  Gnilniund  nii,  ()..  S:  Ton'.'.  1001  Ihe  nppe;  niaiiih'S  ami  1’  velorily 

s:rii''tiiie  heiiealh  Nmllierii  Ansli  i-.lia  fii,r.ii  Lroa'l-haiiil  nhs'-i  vat  ion-,  /'/i;;-,.  I.'arth  y'/(iri’,‘ 
/ill..  8(1.  8.3  OS. 

Keiiliett .  B.L.N,.  1991 .  3' he  removal  nl  free  sin  fa-  e  ini'-:  leO. inns  li  nni  1  h:  '■"-e"inpo[i''nl 

seisino^rain.s,  G'eo;i/ii/,s.  J.  Jut.,  10-1.  1-03  103. 

Keimet.t,  D.  L.  K.  &  Eiif’/lahl.  E.  U.  1091.  Travel  times  for  ('.Inhal  cartluiiiake  hjeation  and  plia- 
assoeiatiou,  6'eo/dii/.s,  J,  /lU,  105.  '!20-4G!i. 

Kennell,  D.  L.  N.,  1901,  lASPEI  1991  Sci.imolugical  Tabic.i.  Bibliotoeh,  Canherra. 


Bibliography 


17-'. 


Kciiiiell.  3.L.N..  1903.  The  di.st.nr.ec  clcpcside:.':'' of  rtgioiiKl  phiusc  fliscriiiiiiNiiil.s  Ihilt.  sciiri}. 

Soc.  Am..  S3,  115h-llGG. 

Kojiiu-tt,,  D.  L.  K..  1905.  Event,  locution  aii'l  sj  in.e  diarart.erisatioii.  in  Monilni my  a 

Conijirchcvsn  c  Test  Ban  Trctiiy.  ccl.s  E.  S.  Husehye  fc  A.  M.  Dainty.  Khnvci  .  Doi'diccht. 

Maiupi  icc.  D.  i:  Silver,  P.  G.,  1993.  Intcip'.etr.tio;;  of  SKS  waves  using  .sample.';  from  tli',’ 
sul)CC''t,iiiciiLa!  lilhn.siihere,  Phij.'^.  Earti-  Plerrl.  hit.,  78.  257-2S0. 

Moatagner,  i:  Tanimoto.  T..  1901.  GioO:-.'.  r.fipcr-iiraiitle  tomography  of  .seismic  velociii".' 
and  ani.socropy,  J.  Geophtjs.  Pc.'i..  86.  20337-20351. 

Mykkeltveit,  S.,  Rlngdal,  F.,  Kvaerna.  T..  E  Alewinc.  R.VV.  (1990)  Application  of  regional  ana;.  ■ 
in  .seismic  verification  research.  Dull.  Scis''.  Soc.  Am,  SO.  1777-lSOO. 

Kalaf,  H-C..  Nakauishi,  1..  sE  Anderson,  D.  L..  19SG.  Measurements  of  mantle-wave  velocities  a;, 
inversion  for  lateral  heterogeneities  and  anisotropy,  3,  Inversion.  J.  Geophys.  Pcs.,  91, 
72C1-7307. 

Nicolas,  A.,  Boucher,  J.  L.,  Bondier.  F.  &:  Mercier.  J.  C.,  1971.  Textures,  structures  and  fabrics 
due  to  solid  state  flow  in  some  European  Ihcizolitcs,  Tcclonophysics,  12,  55-S6. 

Nicolas.  A.,  Bondier,  F.  ,k;  Doullier.  A.  M..  1973.  Mechani.sms  of  flow  in  naturally  and 
e.xperimeiitallv  deformed  peridolites.  At:  .J.  Sci.,  273,  S53-S7C. 

Nicohes,  A.,  i;  Poirier.  J.  P.,  197G.  Cry.stalline  Plasticity  and  Solid-State  Flow  in  Metainorphic 
Rocks.  Whiley,  London,  '137  pp. 

Nicohus,  A..  St  Chri,stcn,sen,  N.  I..  1937.  Forn..''.'.:cj;i  of  nni.sotropy  in  upper-manlle  peridotite  -  a 
review.  Rev.  Geopky.s.,  25.  111-123. 

Pavlidis,  T.,  1971.  Linguistic  analysis  of  v.-av-;:;ri:'.s.  Softicun  Engineenny  \'ot.  2.  (,).  Ton.  cdj. 
Academic  Press,  New  York, 

Pioss,  W.  FI,,  Flannery.  B,  P..  Teukofsky.  S.  .A  .  h  \'cttor!ing  W.  T.  19SS.  Numerical  RccIih's, 
Ciambridgr;  Dnlveisity  Pro.ss. 

R.iitt,  R,  \V..  Shot  ,  G.  G,,  Francis,  T.  J.  G..  vk  Moriis.  G.  B..  1909.  Ani-utropy  of  the  Pacific 
upper  mantle.  J.  Gcophys.  lies..  7d.  3u9:--3it'9. 

Regan.  J.  rG  Anderson,  D.  L,.  1934.  Anisotropic  models  ol  the  t\pi>er  nnintle,  Phy%.  Earth  Piin:> 
Inter..  35.  227-2G3. 

Shearer.  P,  M.,  i;  Orcutl,  J..  19SG.  Comprcs.detK'il  and  shear-wave  atiisotropy  in  the  oceanic 
litlios'  hero.  Geophy.^.  J.  Royal  Astron.  S:-...  ST.  907-1003. 

.S''vci.  P.  G..  3;  Chan,  \V.  \V..  1991.  51ioar-v,,-.-.e  splitting  and  subcuntmental  mantle  defnrmali' 

J.  Gco]ihy.<i.  Res..  96.  lC-129-10451. 

Stockman.  G.,  Kanal.  L.  &  Kyle.  M,  C..  1970.  Structural  pattern  iccognition  of  carotiti  pulse 
waves  using  a  general  waveform  junsiuc  Comm.  A.sso':.  Coniji.  Much.,  19,  688-095. 

Tanimoto.  T.  Anderson.  D.  L..  19S4.  Mappir.g  convection  in  the  mantle.  Geophys.  Res.  Lett.. 
li,  2S7-290. 

Ihomson,  C.  J.,  Kendall,  J-M.,  3;  Guest,  W,  S..  1992.  Geometrical  theory  of  shear-wave  splittir.. 
corrections  top  ray  theory  for  interference  in  isotropic/anisotropic  Ir  nsitions,  Geophys.  J. 
Int.,  lOS,  339-363. 


Bibliography 


) 


H 


■ 


ITr. 


Tong.  C.,  Gviclmund.ssnn,  0.  &  Kennctl,  B.  L.  K..  1994.  Shear  wave  .>:plit.t.iiig  in  refractrtl  wavc^ 
rctuincd  from  thv  upper  maiitie  transition  zone  beneath  nortlieni  Au.slralia.  J.  Ccophij-'i. 

Hen,  99,  15,783-15,797. 

Tong,  C.,  1995.  Characterization  of  .seismic  phases  -  an  automatic  analyst.-r  for  seismograms. 
Gi'-opliiis.  J.  hit.,  123,  937-947,  v 

Tong.  C..  i  Kennett,  D.  L.  N.,  1995.  Townrd.s  the  identification  of  later  sei.sini''  (ilia.scs.  Ccoi'l'i,- 
J.  ha.  123,  94S-95S. 

Tong,  C..  S:  Kennett,  D.  L.  K.,  1990.  Automatic  sei.sinic  event  recognition  anrl  later  pha.se 

identification  for  broiul-b.and  .seismograms.  Bull.  Scimn.  Soc.  Ain.,  subniittcd  in  Koveinber 
1995. 

Vinnils.  L.  P..  Kos.arev,  G.  L.,  i,'.  Makeyeva,  L.  1..  19S4.  Anisotropy'  in  the  lithosphere  from  the 
observations  of  SKS  and  SKKS,  Dokl.  Acad.  Nauk  SSSR,  278,  1335-1339. 

Vinnik,  L.  P.,  Makeyeva,  L.  I.,  Milev,  A.,  .t  Usenko,  A.  Y.,  1992.  Global  p.atterns  of  azimuthal 
anisotroijy  and  deformations  in  the  continental  mantle,  Ccvphy.'i.  J.  hit..  Ill,  433-447. 

5’u,  G.  K,,  ic  Mitchell,  B.  J.,  1979.  Rogionalized  shear-velocity  models  of  the  Pacific  upper  mam  1'' 
from  observed  Love-  and  Rayleigh- wave  dispersion,  Gcophya.  ,/.  11,  A.ition.  Snc..  57,  311-311 


I 


I 


I 


I 


> 


» 


V 


